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PREFACE 


The  David  W.  Taylor  Lectures  were  initiated  as  a living  memorial  to 
our  founder  in  recognition  of  his  many  contributions  to  the  science  of 
naval  architecture  and  naval  hydromechanics.  We  invite  eminent  scientists 
in  fields  closely  related  to  the  Center's  work  to  spend  a few  weeks  with  us, 
to  consult  with  and  advise  our  scientific  staff,  and  to  give  lectures  on 
subjects  of  current  interest. 

Admiral  Taylor  was  a member  of  that  very  select  group,  a naval  officer 
who  was  also  a world  renowned  scientist.  His  research  retains  its 
significance  today.  He  was  a pioneer  in  the  use  of  hydrodynamic  theory  for 
solving  naval  problems,  in  the  use  of  mathematics  for  hull  form  development, 
and  in  his  "global"  treatment  of  ship  resistance.  As  a naval  officer,  he 
founded  and  directed  the  Experimental  Model  Basin.  He  personally 
influenced  the  design  of  virtually  all  ships  of  the  U.  S.  Navy  for  many 
years.  Thus,  he  was  more  than  a scientist,  he  was  a manager  and  working 
engineer  who  helped  shape  a great  fleet. 

This  year  our  lecturer  is  Vice  Admiral  Roger  Brard  of  the  French 
Navy.  As  the  director  of  the  Bassin  d'Essais  des  Carenes,  in  Paris, 
for  many  years,  he  contributed  to  the  design  of  many  ships.  And  his 
scientific  work  in  naval  hydrodynamics  has  been  rewarded  by  the  highest 
recognition  our  profession  can  offer.  Thus,  in  every  way,  his  career 
is  a modern  parallel  to  that  of  Admiral  Taylor.  We  are  most  honored  that 
he  has  agreed  to  be  a David  W.  Taylor  Lecturer. 

I present  Roger  Brard,  Vice  Admiral  of  the  French  Navy,  former 
President  of  the  Academie  des  Sciences  of  France,  naval  architect, 
mathematician,  hydrodynamicist,  and  friend. 
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INTRODUCTION 


May  I first  thank  Dr.  Cummins  for  his  kind  words  of  welcome.  Please 
believe  that  I consider  them  much  too  laudatory.  It  is  a great  honor  and 
privilege  for  me  to  deliver  this  Second  Series  of  David  Taylor  Lectures. 

But  it  is  also  a dangerous  task,  in  that  it  is  nearly  impossible  to  give 
lectures  worthy  of  such  a great  man. 

Admiral  Taylor  combined  a genial  intuition,  a scientific  spirit,  and 
an  astonishing  sense  of  the  physical  possibilities.  His  achievements 
deserve  admiration.  Despite  the  great  strides  made  in  naval  architecture 
since  his  day,  it  is  still  difficult  to  design  ship  forms  superior  to  those 
of  the  Taylor  Standard  Series. 

‘ The  subject  selected  for  these  lectures  concerns  ship  maneuverability. 
This  is  a field  where  theory  occupies  a much  less  important  place  than 
empiricism  because  of  the  complexity  of  the  phenomena  involved  in  the 
maneuvering  of  a ship.  No  complete  theory  yet  exists,  and  one  can  presume 
that  this  will  be  so  for  many  years. 

Nevertheless,  it  can  be  fruitful  to  gather  some  facts  revealed  by 
experiments  and  to  connect  them  with  relatively  simple  theoretical 
concepts.  This  does  not  lead  to  what  a theory  should  be  but  rather  to  an 
introduction  to  such  a theory.  This  is  the  purpose  of  the  present 
lectures. 

~>  Two  facts  seem  to  be  of  great  importance: 

One  is  that  a turning  ship  experiences  a system  of  hydrodynamic 
forces,  one  component  of  which  is  orthogonal  to  the  velocity  of  her 
center  of  gravity.  This  component  must  be  considered  as  a lift.  But, 
in  the  classical  theory  of  lifting  surfaces,  a lift  appears  only  if  free 
vortices  are  shed.  It  follows  that  the  vortex  theory  should  be  the 
basic  tool  of  the  maneuverability  theory. 

The  second  important  fact  is  that  the  inertial  forces  of  the  fluid 
may  depend  not  only  on  the  motions  of  the  body  itself  but  also  on  the 
existing  motion  of  the  fluid  when  the  maneuver  occurs.  This  means  that 
the  maneuverability  theory  has  to  take  into  account  not  only  the  classical 
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added  masses”  but  also  "apparent  added  masses. 
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It  follows  from  the  above  remarks  that  in  order  to  obtain  an 
expression  for  the  classical  added  masses  in  terms  of  bound  vortex  dis- 
tributions, we  must  first  examine  how  the  vortex  theory  applies  when  no 
free  vortices  are  shed  by  the  body.  Then  we  have  to  generalize  the 
classical  lifting  surface  theory  to  the  case  of  bodies  with  a finite 
displacement.  This  must  be  done  even  when  there  is  no  singular  line  on 
the  hull  which  could  play  the  role  of  the  trailing  edge  of  a wing. 

Now,  let  me  outline  briefly  the  contents  of  these  lectures.  They 
are  divided  into  seven  chapters. 

Chapter  1 comments  in  a concise  manner  on  the  concept  (which  will 
be  used  throughout  the  lectures)  of  an  "almost"  inviscid  fluid,  on  the 
state  of  the  art,  and  also  on  the  optimizations  of  the  trajectories  of 
a body  moving  in  water  even  though  the  control  of  such  a body  is  beyond 
the  scope  of  these  lectures. 

Chapter  2 summarizes  the  classical  theory  of  the  motion  of  a body 
moving  through  an  inviscid  fluid.  It  recalls  the  concept  of  the  added 
masses  according  to  which  the  inertial  forces  of  the  fluid  are  derived 
from  its  finite  kinetic  energy. 

Chapter  3 treats  the  same  problem  but  with  the  help  of  vortex 
theory.  The  moving  body  is  replaced  by  a vortex  distribution  kinematically 
equivalent  to  the  body  inside  the  domain  Dg  exterior  to  the  hull  and 
actually  occupied  by  the  fluid.  The  domain  interior  to  the  hull  is 
occupied  by  a fictitious  fluid  mass  which  is  at  rest  with  respect  to 
the  hull.  By  stating  that  the  set  of  fluid  particles  belonging  to 
the  vortex  distribution  is  in  dynamical  equilibrium,  we  obtain  a 
generalized  form  of  the  Kutta-Joukowski  theorem  concerning  the 
hydrodynamic  forces  exerted  on  a bound  vortex  sheet.  Furthermore,  we 
can  now  state  the  expression  for  the  hydrodynamic  pressure  on  the  hull. 

This  may  help  to  throw  light  on  the  concept  of  added  masses  and  apparent 
added  masses. 

Chapter  4 begins  with  a review  of  the  theory  of  thin  wings, 
continues  with  the  theory  of  wings  of  finite  thickness,  and  then  deals 
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with  the  case  of  moving  bodies  which  have  no  real  trailing  edge.  In 
our  opinion,  the  lifting  forces  exerted  on  such  bodies  are  due  to  a 
shedding  of  free  vortices;  this  phenomenon  occurs  because  of  the 
intersection  of  trajectories  of  fluid  points  coming  from  the  suction 
side  with  those  coming  from  the  pressure  side.  Two  examples  are  selected. 

One  could  be  that  of  a submarine  which  moves  in  the  vertical  plane 
and  the  other  that  of  a double  model  moving  in  the  horizontal  plane.  The 
problem  would  be  underdetermined  if  the  continuity  of  the  pressure 
through  the  free  vortex  sheets  did  not  entail  the  continuity  through 
the  shedding  line  of  the  pressure  on  the  hull.  Here,  we  encounter  a 
difficulty  because  the  position  of  the  shedding  line  on  the  hull  has 
to  be  estimated  from  experiments. 

Chapter  5 is  divided  into  two  parts.  The  first  one  deals  with  a vertical 
flat  plate  which  can  be  considered  representing  an  infinitely  thin  double 
model  or  an  infinitely  thin  wing  with  a very  small  aspect  ratio.  In  his 
thesis  Casal  had  studied  this  problem,  with  the  view  to  tackle  a limiting 
case  for  the  maneuvering  ship.  Casal' s results  were  interesting.  Here,  we 
have  examined  some  refinements  for  the  boundary  condition  at  the  plate  to  be 
satisfied  everywhere  on  its  surface.  But  the  problem  has  appeared  more 
complicated  than  the  one  studied  in  Chapter  4.  The  second  part  of  the  Chapter 
is  devoted  again  to  the  latter  problem  but  the  free  surface  effects  are  taken 
into  account  (the  boundary  condition  there  being  linearized). 

Chapter  6 is  devoted  to  small  motions  about  a uniform  motion  of 
translation.  Two  kinds  of  small  motions  are  considered  (1)  those 
analogous  to  the  small  motions  of  a wing  in  the  vertical  plane  and  (2) 
those  analogous  to  the  small  motions  of  a double  model  in  the  horizontal 
plane.  The  motions  of  the  second  kind  are  simpler  with  regard  to  the 
pressure  continuity  condition,  but  they  are  more  complicated  in  that 
the  vortex  wake  is  necessarily  three  dimensional. 

Chapter  7 combines  several  subjects  which  could  not  be  attacked 
within  the  framework  of  a very  general  theory:  interaction  of  the 

appendages  with  the  hull,  small  motions  of  a wing  profile,  and  problems 
more  or  less  solved  or  understood  (e.g.,  snaproll  and  pitchup  of 
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submarines,  effect  of  the  history  of  the  motion  or  equivalency  of  the 
apparent  added  masses  on  the  course  stability  condition).  Lastly, 
attention  is  drawn  to  the  help  that  experimentation,  simulation,  and 
theory  should  bring  to  one  another. 
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CHAPTER  1:  NEEDS  AND  STATE  OF  THE  ART 

According  to  the  French  vocabulary,  the  "maneuverability"  of  a 

given  ship  is,  in  the  mathematical  sense  of  this  term,  the  domain  of  all 

her  possible  steady  motions.  The  term  "handiness"  characterizes  the 

shortness  of  the  time  necessary  to  pass  from  a given  steady  motion  to 

a second  steady  motion  close  to  the  first  one.*  Hence,  "handiness" 

increases  with  the  stability  of  the  second  motion.  In  the  past,  the 

opposite  proposition  was  generally  held. 

Important  results  in  the  domain  of  directional  stability  have  been 

2 

obtained  during  the  last  few  years  or  so  by  Contensou,  Davidson  and 
3 4 

Schiff,  Dieudonne,  and  Grim.**  I myself  studied  the  stability  of 
submarines  in  a vertical  plane  in  a paper  written  20  years  ago,  but 
not  published.***  I had  shown  that  if  the  stability  becomes  negative 
beyond  a certain  critical  speed,  the  ship  tends  to  experience  periodic 
motions  of  considerable  amplitude;  these  must  be  avoided  at  all  costs. 

More  recently,  many  other  authors  have  contributed  considerably  to  progress 


*For  the  distinction  between  "maneuverability"  and  "handiness," 
see  Roy.^ 

**Unpublished  reports  of  0.  Grim,  Hamburg  Schiffbau  - Versuchanstalt . 

***Brard,  R. , "La  tenue  de  plongee  des  sous-marlns , " prepared  at  the 
Bassin  d'Essais  des  Car^nes,  Paris,  in  November  1953  but  not  published. 

^Roy,  M.,  "Le  probleme  de  la  stabllite  des  regimes  de  vol," 

Bull.  Ass.  Tech.  Mar.  et  Aeron.,  Paris  (1931). 

2 _ 

Contensou,  P.,  "Mecanique  du  Navire  en  route  et  en  giration," 

Bull.  Ass.  Tech.  Mar.  et  Aeron.,  Paris  (1938). 

3 

Davidson,  K.  S.  M.  and  C.  Schiff,"  Turning  and  Coursekeeping 
Qualities  of  Ships,"  Trans.  Soc.  Nav.  Arch.  Mar.  Eng.  (1946). 

4 

Dieudonne,  J.,  "Note  sur  la  stabilite  des  regimes  de  route  des 
navires,"  Bull.  Ass.  Tech.  Mar.  et  Aeron.,  Paris  (1949). 


5 


in  the  field  of  directional  stability  and  control  of  bodies  moving  in 
water  (see,  for  instance,  the  papers  cited  as  references  as  well  as  the 
Proceedings  of  the  International  Towing  Tank  Conference  (ITTC)  and  those 
presented  at  the  seminar  held  in  London  during  1972).“* 

Optimization  of  the  trajectories  of  a given  ship  should  be  the  goal  of 
research  on  maneuvering  qualities.  But  this  problem  would  be  meaningless 
if  the  domain  of  all  her  possible  study  motions  was  too  small.  On  the 
other  hand,  such  an  optimization  is  practically  impossible  if  the  equations 
of  all  the  possible  motions  are  not  known  accurately. 

Unfortunately  this  is  effectively  the  case.  Because  of  the  short- 
comings of  the  theory,  it  is  necessary  to  resort  to  experiments  on 
captive  and  on  free-running  models.  But  to  obtain  good  agreement  between 
the  measured  trajectories  and  those  derived  from  equations,  it  has 
become  necessary  to  introduce  a considerable  number  of  terms  into  the 
equations.  This  leads  to  expansions  whose  uniqueness  is  questionable 
since  equivalent  predictions  can  be  obtained  by  simultaneously  altering 
several  terms.  Only  the  first  terms  of  such  expansions  have  physical 
meaning. 

Theoretical  research  must  be  pursued  actively  in  order  to  improve 
this  rather  poor  situation.  This  is  not  an  easy  task  for  every  theory 
implies  assumptions. 

In  the  present  lectures,  we  consider  only  the  case  of  liquid 
unbounded  in  all  directions  (except  the  vertical  direction  when  a free 
surface  exists)  and  we  assume  that  the  liquid  is  "almost"  inviscid.  This 
means  that  liquid  adheres  to  the  hull  of  the  ship  but  that  the  other 
effects  of  viscosity  are  ignored.  This  considerably  simplifies  the 
problem  and  allows  us  to  utilize  vortex  theory;  thus,  the  system  of 
hydrodynamic  forces  exerted  on  the  body  can  include  a lift.  But  the 
assumption  prevents  us  from  taking  into  consideration  the  phenomenon 
of  separation  due  to  strong  adverse  pressure  gradients.  Furthermore, 
the  assumption  does  not  provide  us  with  means  to  determine  the  exact 


^Proceedings  of  the  International  Symposium  on  the  Directional  Stability 
and  Control  of  Bodies  Moving  in  Water,  London  (1972). 
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CHAPTER  2:  MOTION  OF  A SOLID  THROUCH  AN  UNBOUNDED  LIQUID 

This  theory  Is  classical;*  It  is  reviewed  here  for  comparison  with 

results  derived  from  the  bound  vortex  theory. 

We  consider  a right-handed  system  S of  axes  0(x^ , x^)  moving 

with  the  body.  Let  u^,  u2»  u^  denote  the  components  on  the  moving  axes 

of  the  velocity  V (0)  of  0 and  u.  , u , u,  those  of  the  angular  velocity 
h 4 j b ^ 

of  S.  Let  I be  the  momentum  of  the  fluid,  J the  moment  of  the  momentum 

about  0,  and  K the  moment  of  momentum  about  the  fixed  point  O'  with  which 

0 coincides  at  the  time  t under  consideration.  Let  (£^,  £2>  £3)  and 

(£. , £c,  {,,)  be  the  components  of  I and  J,  respectively. 

4 0 0 

As  regards  the  force,  the  hydrodynamic  forces  exerted  on  the  body 

are 


(xr  x2. 


x3)  =- 


->■ 

dl 

dt 


and  as  regards  the  moment  about  O',  they  are 


V 


V = - 


dK 

dt 


We  may  express  the  X's  in  terms  of  the  £'s.  For  instance, 

X1  = " (ife+  U5C3  " U6C2>' 

^4 

X4  = <‘~“dt  + U5^6  " U6^5  + U2?3  " U3C2^ 

On  the  other  hand,  we  can  show  that  the  kinetic  energy  T is  given  by 


*After  Kirchhoff,  Kelvin,  and  Lamb  analyses.  For  more  detail, 
see  Lamb.^ 

^Lamb,  Sir  Horace,  "Hydrodynamics,"  Sixth  Edition,  Cambridge  University 
Press  (1932),  Chapter  VI. 
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position  on  the  hull  of  the  line  of  shedding  of  the  free  vortices.  It 
follows  that,  if  a slight  change  in  its  position  can  entail  an  important 
alteration  of  the  calculated  system  of  hydrodynamic  forces  exerted  on  the 
body,  then  it  is  necessary  to  determine  the  position  of  this  line 
experimentally. 


7 


where  $ is  the  velocity  potential  inside  the  domain  D exterior  to  the 

e 

hull.  It  can  also  be  written  in  the  form: 


t . . u , u , 
ij  i J 


where 


ij 


°/L  * i 


dp  . 

- dZ,  <p  = <t>(u,,  u„,  ...  u.)  for  u.  = 1 and 
dn  i 1 i o i 

Uj  = 0 (every  j t i) 


The  t ^ ' s are  the  elements  of  a square  symmetric  matrix  T.  They  depend 
only  on  the  geometry  of  the  hull. 

9T 

Since  we  can  show  that  -r — = then  by  substituting  the  £.'s  into 

3uj  J j 

the  expresssions  for  the  X's,  we  obtain  the  following  Kirchhoff  formulas: 


By  solving  six  Neumann  exterior  problems,  we  can  determine  the 
then  the  t ^ . 

In  the  formula  for  X^,  the  first  term  on  the  right  side. 


and 


namely 


d 3X 

- — -jj—  , represents  the  effect  of  the  inertia  of  the  fluid,  i.e.,  the 

effect  of  the  so-called  "added  masses."  The  other  terms  give  the 
expression  for  the  hydrodynamic  forces  exerted  on  the  body  when  its 
motion  is  uniform  and  therefore  when  the  fluid  motion  is  steady  with 
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respect  to  the  body.  The  resulting  system  of  hydrodynamic  forces 
acting  on  the  body  may  thus  be  written  in  the  form: 


S 


d 


+ S 

q.  s. 


The  part  S.  in  the  contribution  from  the  added  masses  and  the  part  S 

in  q.s. 

would  exist  alone  if  the  motion  of  the  body  were  uniform. 

In  a similar  manner  we  can  derive  the  system  of  forces  from  the 

g 

kinetic  energy  T of  the  solid. 


Ss 


s®  + s 

in  q.s 


which  should  be  exerted  on  it  in  order  to  define  its  motion  by  the  six 
functions  u (t),  u„(t),  ...  u,(t).  In  fact  the  system  of  forces  exerted  on 

1 L d O 

the  solid  is  the  sum  of  S , of  the  forces  due  to  hydrostatic  pressures 
on  the  hull  and  the  system  due  to  external  forces.  By  expressing  this 
sum  as  equivalent  to  the  system  Ss,  we  obtain  a set  of  six  partial 
differential  equations  which  are  the  equations  of  the  motion  of  the 
solid . 

THE  D'ALEMBERT  PARADOX 

If  u.  » uc  * u,  = 0,  that  is  if  U = 0,  and  furthermore  if  u , u„,  u 
4 j o fc.  I z j 

are  constants,  we  obtain 


X1  - X2  = x3  = 0 

« 

but,  in  general,  the  moment  (X.,  Xc,  X,)  is  not  zero.  The  system  of 

4 j b 

hydrodynamic  forces  exerted  on  the  solid  reduces  to  a couple.  This  is 
the  famous  d'Alembert  paradox. 
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STABLE  AND  UNSTABLE  STEADY  MOTIONS 

Let  us  assume  that  the  system  of  external  forces  Is  equivalent  to  zero. 
By  specifying  that  u , ...  u,  are  constants,  we  obtain  the  equations  of 
all  the  steady  motions  of  the  solid.  More  particularly,  we  obtain  the 
equations  of  the  steady  motions  of  translation.  By  doing  this  and 
studying  the  hydrodynamic  couple  exerted  on  the  solid,  Lamb  showed  that 
in  the  case  of  an  ellipsoid  whose  axes  coincide  with  its  axes  of  inertia, 
the  only  stable  translation  is  that  in  the  direction  of  its  smallest 
axis  of  symmetry. 

Actually,  however,  experiments  have  shown  that  other  translations 
can  be  stable.  The  explanation  of  this  apparent  contradiction  is  that 
the  hydrodynamic  forces  on  a body  in  a uniform  motion  of  translation 
do  not  reduce  to  a couple  alone.  This  is  a consequence  of  the  shedding 
of  free  vortices.  We  will  see  in  Chapter  4 that  this  phenomenon  alters 
both  the  quasi-steady  system  of  forces  and  those  due  to  inertial  effects. 
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CHAPTER  3:  MOTION  OF  A SOLID  THROl/CH  AN  UNBOUNDED,  I NV I SC ID  FLUID 

ACCORDING  TO  THE  BOUND  VORTEX  THEORY 

KINEMATICS 

Vortex  Filaments,  Vortex  Tubes,  Vortex  Sheets 

If  the  velocity  V of  a fluid  is  continuous  and  continuously 
differentiable  in  a certain  domain  D,  we  can  define  the  vorticity  w 
at  every  point  M belonging  to  that  domain.  The  vorticity  at  M is 
defined  by 


w(M)  = curl  V (M) 

The  properties  of  a given  distribution  V = (D,co)  of  u>  within  D are 
derived  from  the  Stokes  theorem.  If  C is  a closed  circuit  located 

-f 

within  D,  the  circulation  of  V in  that  circuit  is  equal  to  the  flux  of 
co  through  any  open  surface  S whose  edge  coincides  with  C: 

T (C)  = f V*ds  = / / (0*n  dS 

JC  Js 

— ► 

The  unit  vector  n normal  to  S is  in  the  positive  direction  with  the 
sense  selected  on  C. 

A vortex  filament  L is  tangent  to  co  at  every  point,  and  a vortex 
tube  T is  the  surface  generated  when  the  vortex  filaments  L intersect 
a given  closed  contour  C.  It  follows  from  the  above  theorem  that  if 
reduces  to  C by  a continuous  deformation,  then  the  circulation  r(C^) 
in  any  closed  contour  on  T is  a constant  equal  to  F (C) . F (C ) is 
termed  the  intensity  of  the  tube  T.  A tube  T can  contain  a unique 
filament  L and  its  intensity  can  nevertheless  be  different  from 

-F  _ -F 

zero.  In  that  case,  w is  infinite  on  L and  V is  no  longer  continuously 
differentiable  at  every  point  of  D;  L then  has  to  be  considered  as  an 
infinitely  thin  vortex  tube. 
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A vortex  sheet  can  also  be  defined  by  a limiting  process.  Such  a 
sheet  is  a domain  of  very  small  thickness  e.  The  vorticity  inside  the 

sheet  is  w = — ; T is  finite  when  £ goes  to  zero.  The  sheet  reduces  to 

a surface  E and  T is  tangent  to  E at  every  point  (see  Figure  1). 

Let  P be  a point  belonging  to  £.  Let  £+  and  £ be  the  two  sides  of  £. 

-►  

Let  n denote  the  unit  vector  normal  to  £ in  the  direction  from  £ toward 

-f-  -►->-►  •+ 

£ , and  (n,  0,  T)  a right-handed  system  of  three  unit  vectors,  T being 
in  the  direction  of  T.  We  may  consider  on  £ the  vortex  filaments  £ and 
the  lines  C tangent  to  T and  respectively,  and  define  the  elements  of 
arc  ds  on  £ and  da  on  C so  that  ds  > 0 in  the  direction  of  T and 
da  > 0 in  the  direction  of  S.  The  part  of  the  sheet  located  between  two 
lines  i-(o),  L (a  + da)  is  an  infinitely  flat  tube  of  thickness  £ and  of 
intensity 


dF  = T*T  da 


According  to  the  Stokes  theorem 

datecv*  - v")]p  = dr 


and 


(V  - V )p  = (T  A n) 

“► 

Conversely,  any  surface  through  which  V is  discontinuous  can  be  considered 
as  the  support  of  a vortex  sheet. 

Since  the  intensity  of  a vortex  filament  is  a constant,  a vortex 
filament  cannot  begin  or  end  in  the  fluid.  The  support  of  £ is  a 
closed  contour  or  its  two  ends  are  located  on  the  boundary  of  the  fluid 
domain  (possibly  at  infinity). 
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The  Poincare  and  Biot  and  Savart  Formulas 

Let  denote  the  domain  interior  to  a closed  surface  S and  n the 

unit  vector  normal  to  S in  the  inward  direction.  Inside  D^,  the 
velocity  V is  supposed  to  be  continuously  differentiable.  By  applying 
the  classical  formula 


curl  (curl  A)  = V(div  A)  - AA 


(where  A is  the  Laplace  operator)  to  the  vector 


A(M) 


V(M’) 

MM' 


dD.(M') 


we  can  readily  obtain  the  famous  Poincare  formula: 


curl 


u// 


(nAV)M^ 

MM* 


dS(M')  + 


f 


iu(M'  ) 
MM’ 


dDi(M 


'>] 


+ V 


(n*V)M^ 

MM' 


dS(M') 


_1  fff  div  V(M') 
4 tt  JJJn  ~ MM' 


dDi(M 


■>] 


V (M)  if  M e D. 
0 if  M e D 

e 


(3.1) 


In  this  formula,  D is  the  domain  exterior  to  S,  and  M'  is  defined  by 
e i 1 

M7^  = (0+). 
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In  the  following,  we  deal  with  liquids  and  div  V = 0.  Formula  (3.1) 
gives  the  solution  of  the  partial  differential  equation 

= curl  V,  inside  D. 


when  V is  given  on  the  boundary  S of  D.. 

By  means  of  a limiting  process,  it  is  seen  that  if  S is  at  infinity, 
and  if  contains  a unique  vortex  filament  L of  intensity  I’ , then 


V(M)  = curl 


4TT 


fL 


ds(M') 

MM' 


(3.2) 


This  is  the  Biot-Savart  formula. 


Vortex  Distribution  Kinematically  Equivalent  to  the  Hull  of  a Body 

Let  £ be  the  surface  of  the  body  and  D.,  D the  domains  interior  and 

l e 

exterior  to  £.  Let  £.,  £ denote  the  two  sides  of  £.  The  unit  vector 
- 1 e 
n normal  to  £ is  in  the  inward  direction. 

Our  purpose  is  to  show  that  there  exists  one  vortex  distribution 

inside  and  on  £ which  is  equivalent  to  the  body  inside  and 

which  satisfies  the  condition  that  the  relative  velocity  is  zero  at  £^. 

To  that  end,  we  apply  the  Poincare  formula  inside  Dg.  At  infinity, 

| V | = 0^— R is  the  distance  from  £.  Since  n is  in  the  outward 


direction  with  respect  to  and  w 


0 inside  D 


we  readily  obtain: 


curl 


(-nAV)M' 

MM'  * d£(M')  +V 


(trv)M' 


’ V (M)  if  M e 
, 0 if  M e Di 
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Let  V denote  the  velocity  at  a point  moving  with  the  body.  The 

b, 

■*  — ► 

Poincare  formula  applied  to  V inside  D.  gives: 

E 1 


curl 


_i  rr  <"%)»; 

47T  JJ  MM' 


dE(M')  + V 


(n-VE)M’ 

MM' 


d£(M') 


+ curl 


ML 


2^  (M' ) 


^ JJJV  mm'  dDi(M,)  = 


V if  M £ D . 
E i 


0 if  M £ D 


By  adding  the  two  formulas,  we  have: 


curl 


JL  rr  (~nAVR)Me 
4it  JJ.  MM ' 


dl  (M 


,)+/l 


2&,(M') 


“MM—  dDi(M,) 


{V  (M)  if  M £ D 

. ' 

V£(M)  if  M £ Di 


(3.3) 


-►  -f 

where  V is  the  relative  velocity  V - V of  a fluid  point. 
R E 

This  shows  that  there  exists  a vortex  distribution 


P = (e,  f ) + (Di,  2^e)  with  T(M')  = (-nAVR)M^ 


(3.4) 


which  fulfills  the  requirement.  Since  function  T is  defined  on  E, 
we  necessarily  have 


3_2_[_ 

9o3s 


(3.5) 
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Effective  Determination  of  the  Vortex  Distribution 
-> 

The  vector  T is  determined  by  the  condition  that  the  left  side  of 

— ► 

Equation  3.3  be  equal  to  V when  M e D,  . This  is  satisfied  if  it  is 

- E 1 

equal  to  V (M . ) when  M E L . This  leads  to  the  vectorial  Fredholm 
b 1 1 

equation  of  the  second  kind 


1 

2 


(nYfi 


+ curl 


ISh. 

MP 


— dl(P) 


= V£(M,  t) 


2^e(P,  t) 
MP 


dDi(P) 


(3.6) 


This  equation  is  singular  for  the  vector  T(M,  t)  = X*n^,  where  X is  a 
constant,  cancels  the  left  side. 

The  discussion  of  the  equation  is  rather  tedious  and  it  is  simpler 
to  proceed  in  another  manner.  First  we  replace  V by  the  distribution 


V = with 


i . (e.  + CD,,  2Se),  V2  - (e.  £) 


(3.7) 


The  condition  to  be  satisfied  by  T^  is  that  the  velocity  induced  by 
V is  zero  inside  D^.  If  this  is  possible,  there  exists  within  D^  a 
velocity  potential  such  that 


V'  = V„  - V = V<t ' (within  D.) 
hi  1 


(3.8) 


The  discontinuity  of  through  £ is  tangent  to  E.  Consequently, 
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3$'  _ -*■  ?.  v 

3^  = n'VE  °n  Ei 


(3.9) 


It  follows  that  $ ' is  the  solution  of  a Neumann  Interior  problem. 
This  solution  exists  since  = an^  *-s  unique. 


Conversely,  when  <P ' is  known,  is  determined  inside  by 
Equation  (3.8)  and  we  have 


T1(M)  = nAV  (Mt)  on  Z 


(3.10) 


Let  be  the  velocity  induced  by  V . We  have 


v2  - V 


(3.11) 


inside  D (since  V = 0 in  that  domain).  This  velocity  is  irrotational 
e 1 

also  inside  D . Hence  is  equivalent  to  a normal  doublet  distribution 
on  E : 


V2  'v  (Z,  U'n) 


(3.12) 


The  density  p ' of  this  distribution  must  be  determined  so  that  it  generates 
inside  a velocity  potential  $2  equal  to  this  entails 


V = Vi  + V?  = \ + <VE  - Vx)  = VE  within 


It  follows  that 
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VM>  - -s//--"-*  5^MFdZ(P> 


" \ Vi’(M) 


f^s?dE<p>  ’ 


M C Z 


MM.  = 1^(0+) 


Now  we  have 


T2(M)'  = " nMA[V(Me)  ' V’(Mi)] 


The  boundary  condition  on  Z^,  namely, 


n • V = n • V , on  E 

L 6 


is  satisfied  since  the  normal  component  of  is  continuous  through 

Z and  equal  to  n*V  on  E . 

t 1 

The  above  potential  *s  the  solution  of  a Dirichlet  interior 
problem,  and  (3.13)  has  a unique  solution.  We  have 


4>(M  ) = 4*  (M  ) = <t>  (M.)  + U ' (M)  = 4>'(M.)  + p'(M) 
e i e i i 1 


The  solution  of  the  problem  is 


V - (Z,  | ) + (Dlt  2$e)  with  T - Tx  + 


It  follows  from  the  definition  of  4'  that  at  a point  M located 
inside  Dj  and  moving  with  the  body,  4>'  is  a linear  function  of  the 
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components  u^(t).  We  have 


4>'(M,t)  = 


Uj  (t)  4>  (M) 


(3.17) 


where  the  functions  are  independent  of  time.  Furthermore 
follows  from  (3.13)  that  y'  may  be  written  in  the  form 


it  also 


6 

U'(M,t)  = £ uj(t)  Uj(M)  °-18) 

j = l 


where  the  functions  yj  are  independent 


of  time. 


Generalization 

It  may  be  that  the  body  is  moving  in  an  incident  flow  whose  velocity 

is 


Vq(M,  t)  = V<t*0(M,  t) 

In  this  case,  the  distribution  becomes: 

V2  * (Z,  (y'  + y " )n) 

and  the  potential 


(3.19) 


(3.20) 


^ li (W'  + “">o  ► dE(P) 


P,t  3np  MP 


(3.21) 
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must  satisfy  the  boundary  condition 


$2(M1)  = $’(1^)  - <5>0(M1)  on  St. 


(3.22) 


Let  G be  the  resolvent  kernel  of  Equation  (3.13).  We  have 


P'(M,  t)  = £ ut(t)  U'(M)  = jj  G (M,  P)  *'(P,  t)  dE(P) 

j E 

U"(M,  t)  = //g<m,  P)  [ - * (P,  t)[  dl(P) 


(3.23) 


DYNAMICS 

Forces  Exerted  on  a Vortex  Distribution 

As  in  Chapter  2,  we  consider  two  right-handed  systems  of  axes  S,  S', 
namely  the  system  0(x^,  x2>  x^)  which  moves  with  the  body  and  the  system 
0'(x^,  x2>  x^)  which  is  fixed  in  space.  System  S'  is  selected  so  that 
it  coincides  with  S at  a certain  time  t.  Let  M'  be  the  position  at  t 
of  the  point  M which  moves  with  S,  and  V , Y,.be  respectively  the  velocity 

Ej  Ej 

— ► 

and  the  acceleration  of  M.  Furthermore  V(M,  t)  is  the  velocity  of  the 

fluid  point  P located  at  M at  time  t,  and  V = V - V is  its  relative 

R E 

velocity. 

The  force  F per  unit  mass  of  the  fluid  is  defined  in  system  S'  as 

F =*  F(M\  t) 

According  to  the  famous  Helmholtz  theorem,  the  vorticity  u at  a 
point  M at  time  t is,  in  fact,  a property  of  the  fluid  point  P located 
at  M at  t.  This  follows  from  the  Helmholtz  equation. 
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9 u)  . d u)  , , / u>  „ \ 

U o c>  * dT  p <”•  '>  ■ (p -y 


(3.24) 


where  p is  the  mass  density  of  the  fluid.  For  this  equation  to  hold, 
it  is  necessary  that  the  fluid  be  inviscid  and  that  F be  the  gradient 
of  a certain  function  U of  M'  and  t: 

F(M\  t)  = VU(M\  t)  (3.25) 


If  so,  a vortex  filament  moves  with  the  fluid  and  its  intensity  is 
independent  of  time.  Consequently,  no  external  force  is  exerted  on 
the  fluid  points  belonging  to  the  vortex  filament.  One  can  say  that 
the  vortex  filament  is  free. 

Let  us  come  back  to  the  vortex  distribution 


which  has  been  determined  earlier  so  that  it  generates  inside  Dg  a 
fluid  motion  identical  with  that  due  to  the  motion  of  the  solid.  The 
fluid  points  which  belong  to  the  distribution  V do  not  move  with  the 
fluid.  They  do  not  satisfy  Equation  (3.24)  because  the  adherence  forces 
which  keep  them  at  rest  with  respect  to  the  hull  are  due  to  some 
viscosity  effect  and  do  not  meet  the  requirement  of  Equation  (3.25). 

The  hydrodynamic  forces  exerted  on  the  fluid  points  of  V will  be 
termed  the  hydrodynamic  forces  exerted  on  the  distribution  itself. 

This  misuse  of  language  is  adopted  for  the  sake  of  brevity.  These 
hydrodynamic  forces  are  obviously  connected  with  the  hydrodynamic 
forces  exerted  on  the  moving  body  itself.  To  determine  the  latter, 
we  have  to  begin  by  determining  the  former. 
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Euler  Equation  in  the  Moving  System  of  Axes 

Let  p be  the  pressure,  p the  hydrostatic  pressure,  and  P . the 

s a 

hydrodynamic  pressure: 


Pg  - & , Pd  = P ~ Ps  = P ~ PU  (3.26) 

-►  -► 

Let  V be  the  relative  acceleration  and  U)  the  relative  vorticity.  Then, 
R K 

- d*R  - * 

TR  ' dt“  ' “R  - Curl  VR  ■ <“  - <’-27> 

and  by  the  Coriolis  theorem,  we  have 


Y 


dV 

dt 


+ 2V5r 


where 


The  Euler  equation  gives: 


<’r  + + ’fyfy 


- 

[it 


* "raJr  + 7 


(N)] 


<VE  + 2VVr) 


Let  A denote  the  axis  of  the  helicoidal  motion  tangent  to  the 
motion  of  the  system  of  axes  S,  and  let  A(M,  t)  be  the  vector  orthogonal 
to  A,  with  its  origin  on  A and  its  end  at  M.  The  velocity  of  any  point 
A located  on  A and  moving  with  S is  a vector  Vg(A)  parallel  to  A and 


23 


independent  of  the  position  of  A on  A.  We  obtain 


V (M)  = V (A)  + fy* 


9V, 


Y„(M,  t)  - jf-  (M,  t)  + (Vk(M)-V)  Vf(M) 


and 


(v£*V)  vE  = V^e(a)  + W 


3 a («_M) 

fci  E. 


-af 


V(  - \ nV) 


Hence  the  Euler  equation  becomes: 


- Vp, 

P d 


- 5T  + ,(  f Wt  - v<  i V 


1 „ 2.2, 


1 ..2, 


This  equation  reduces  to 


Ip  = 

P ’ d 


(M*  > t)  + t) 


when 


V (M,  t)  = V4>(M',  t) 
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(3.28) 


(3.29) 


(3.30) 


Hydrodynamic  Forces  on  an  Element  of  Vortex 
Sheet  or  on  an  Arc  of  Vortex  Filament 

Vortex  Sheet.  An  element  of  vortex  sheet  is  the  domain  dD  bounded  by 
two  elements  of  surfaces  dE+  and  dE  ; dE+  is  derived  from  dE  by  a 
translation  n£.  Let  de  denote  the  set  of  fluid  points  located  inside 
dE.  Its  momentum  goes  to  zero  with  £ and  so  does  its  derivative.  Hence, 
the  system  of  forces  exerted  on  dE  is  equivalent  to  zero.  This  system 
consists  of  the  forces  exerted  by  the  adjacent  sets  of  fluid  points 
and  of  the  binding  force  d^.  When  £ goes  to  zero,  we  obtain 


d f*  - (p+  - p“)  n dE  = 0 (3.31) 

The  binding  force  is  thus  normal  to  the  surface  E of  the  bound  vortex 
sheet  and  the  pressure  is  discontinuous  through  the  sheet.  If  the 
vortex  sheet  is  free,  d?  = 0 and 

P+  = P~  (3.32) 

on  the  sides  of  a free  vortex  sheet. 

Vortex  Filament.  Let  be  an  arc  of  length  6 included  in  the 

bi  si 

bound  part  L of  a vortex  filament  L ; see  Figure  3.  Let  A denote 
b 

the  midpoint  of  . The  length  6 is  small  enough  for  L to  be 
bi  si  bi 

regarded  as  a segment  of  the  - axis.  Let  (£,  H,  C)  denote  a 
moving,  right-handed  system  of  axes  and  (r,  0,  O a system  of  polar 
coordinates,  r being  the  distance  from  the  £ - axis.  In  fact,  L 
is  the  limit  of  a segment  of  vortex  tube  with  a very  small  radius 
r^;  its  intensity  is 


f - 


2 it  r. 

JJ 


u)(r,  9)r  dr  d0 


0 0 


(3.33) 
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The  relative  velocity  V may  be  written  in  the  form 

K 


V,  = V„  + (V'  - 5V')  + <$V' 
R Ro 


*■  I I I 

where  is  due  to  all  causes  other  than  t V'  to  1 , and6V'  to  *-  . 

Ro  ^ bi 

The  momentum  i(dE,)  of  the  set  dE,  of  fluid  points  belonging  to  L tends 
11  bi 

to  zero  with  r^  even  when  i remains  finite.  Let  us  put 


V,  = lim  V„  + (V'  - 6v') 
1 r+0  R° 


This  limit  is  the  incident  velocity  on  L,  . Let  D1  be  the  domain: 
bi 

D'  = | r,  0,  C;  0<r<R,  O<0<2it,  ~ f 6s1Si  f 6s1 1 O 

This  domain  does  not  include  the  arc  L.  . Let  I(dE')  be  the  momentum 

b l 

of  the  set  dE'  of  fluid  points  located  within  D'.  We  have: 


k ;<dE'>  *3 


Tt  '(dE'  + dEl>  ’ PJJVR<”'VR 


) dS  + 0(6S  R ) 


where 


S = S + 


• £+  + , SL  = | r = R,  6e  [o.  2tt)  , | C|<  j 6s1  J 

£—  = | 0^r<R  , 0c  J^o,  2tt),  C = + -|  <5S] 


while  n Is  normal  to  S in  the  outward  direction. 
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Let  u^,  v^»  w.  denote  the  components  of  V.  in  the  H,  and  4 
directions,  respectively.  We  have 


lim  §7  I(dE’) 
RiO  dt 


1 im 
R^O 


2tt 

o i f (-  r “in2e  ‘x + r c”s2°  *,) R * 6>1 


(3.37) 


4 p r AV . 6s, 
2 l 1 


Now 


let  - dF  denote  the  force  exerted  by  dE1  on  dE’ . By  the  momentum 


theorem,  we  have 


dt 


I(dE’)  = - dF  + 


li- 


Pn  dS 


(3.38) 


From  (3.28),  we  derive 


1 

P 


^ curl  VrJa(V.  + <$V’) 


V i (V.  + 6v’)Z 


within  D' 


(3.39) 


When  R decreases,  VRq  , V^, 


9Vj 

at 


and  the  first  three  terms  of  Equation 


can  be 
(3.39) 


considered  constant  vectors 
contribute  nothing  to 


— I(dE').  Furthermore, 


6V'  * 


2TTr 


A r 


(3.40) 
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Hence 


d 


(i^AR)  - 


1 

2 


on  Sj 


The  contributions  from  I and  Z cancel  each  other  by  symmetry.  We  thus 
obtain 


lim  7-  I(dE')  = - dF  - 4 p ?AV.  6s. 
RIO  dt  12  11 


Finally,  comparing  with  Equation  (3.37),  we  have 


dFT  = - p ?AV1 


(3.41) 


(3.42) 


This  formula  yields  the  expression  for  the  force  exerted  by  the  incident 
flow  on  an  arc  6s^  of  a bound  vortex  filament.  It  can  be  considered 
as  the  Kutta-Joukowski  theorem. 

System  of  Forces  Associated  with  the  Vortex 
Distribution  V Kinematically  Equivalent  to 
the  Moving  Body. 

The  vortex  distribution  kinematically  equivalent  to  the  moving  body  is 

taken  in  the  form  (3.7).  Its  properties  are  described  by  equations  (3.8)- 

(3.23).  The  fluid  points  belonging  to  V do  not  move  with  the  liquid,  but  with 

the  body.  Consequently,  Helmholtz's  theorem  does  not  apply.  The  vortex 

filaments  whose  union  makes  the  vortex  distribution  and  the  fluid  taking 

part  in  the  general  motion  exert  a mutual  action  on  each  other. 

We  shall  first  assume  that  u. , u„,...,  u,  are  constants  and  therefore 

12  o 

that  the  relative  motion  is  steady. 
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Let  £ > E , be  the  external  and  internal  sides  of  the  hull  surface  £, 
e 1 

respectively.  The  fluid  points  belonging  to  the  vortex  sheet  are  located 

between  £.  and  £ . Let  M be  a point  on  £,  and  M. , M the  points  located  on 
i e r i e 

£.  and  on  £ and  defined  by 
i e 

MMa  = 1^(0+)  , MMe  = -r^(0+).  (3.43) 

The  incident  velocity  on  the  set  dE  of  fluid  points  located  inside  the 
volume  d£  x e is 

VM)  ° ^ (3,44) 

where  V is  the  relative  velocity.  According  to  (3.31)  and  (3.42),  the  hydro- 
dynamic  force  exerted  on  dE  is  given  by 

dFx  = (pd(Me)  - Pd(M1))nM  dj.(M)  = - p T (M)  A ^(M)d£(M) 

= “ 5 pVR  (Me)  dZ(M)  (3.45) 

Let  denote  the  system  of  forces  dF^,.  We  write  symbolically 

ST  = (£,  dFT).  (3.46) 

Similarly,  the  system  of  hydrodynamic  forces  exerted  on  the  body  may  be 
written  as 

* Sd  ■ <V  Pd».>  » 0-c)  »•«) 

Inside  , the  u^'s  being  constant,  we  have: 

piyPd  = ^E=  “e*"  V^Er2)  (3-48) 

whence 

p (M)  - ; p Sip  r“  + constant  (3.48') 

d Z t 

The  system  of  forces  exerted  on  the  vortex  sheet  by  the  liquid  located  inside 
is 

(_Pd(V"M  dV  " (_pi3E  ^(C)dDi)  = ~ 1"“e  = " FC  (3,49) 
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C being  the  center  of  volume  of  D^,  m = pD^  the  mass  of  the  body,  and  the 
vector  with  its  origin  is  on  the  axis  of  helicoidal  motion  of  the  body  and 
its  end  at  C.  The  system  (3.49)  reduces  to  a unique  force,  namely  the  centri- 
petal force  acting  on  the  mass  m. 

One  observes  that 

Pd(Me)  = 3 poj  r2  (M)  - \ V2  (Me),  (3.50) 

34> 

which  is  in  agreement  with  (3.29)  since  the  relative  motion  is  steady  ( = 

dt 

0 at  any  point  fixed  with  respect  to  the  body  and  located  inside  Dg). 

We  have  seen  that  the  velocity  potential  $ inside  Dg  is  identical  with 

the  potential  due  to  a normal  doublet  distribution  of  density  p'  on  Z. 

This  density  is  determined  by  the  Fredholm  equation  expressing  that  the 

interior  determination  of  coincides  with  the  velocity  potential  <J>'  such 

that  V4>'  • V - V.  (equation  (3.8)).  We  have,  in  general, 

Cj  I 

*'  (M,  t)  = Z u (t)  *'  (M)  (3.17) 

i 3 3 

and 


with 


$(M  , t)  = (M  , t)  + p' (M,  t) 
e i 

u'(M,  t)  = Z u (t)  p’(M)  (3.18) 

j J j 


The  vorte:  theory  gives  means  for  calculating  since  the  six  potentials 
4 j ' s follow  from  the  solution  of  the  Neumann  interior  problems  expressing  that 


3 $ ' 

-5S1  <HC  « ■ " • Vi  ™ Ei  <3-51> 


and  from  the  six  Dirichlet  interior  problems  expressing  that 

*j(M1,  t)  = (Mi,  t)  on  Z±  (3.52) 

Let  us  consider  now  the  case  when  the  motion  of  the  body  is  not  uniform. 
Equations  (3.48)  - (3.50)  no  longer  hold.  But  (3.50)  has  to  be  replaced  by 
(3.29)  and  one  has 

(M  , t)  = Z u (t)  (*'  (M  ) + |i ' (M) ) (3.53) 

3t  e j j j i j 
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where  u,  (t)  =*-r— u.(t).  This  gives 
J at  j 


with 


S,  - S,  + S - F 
d in  T C 


(3.54) 


5in  = ("P  Uj(Mi>  + uj^}  dI(M)) 


(~p\  Vj[(M)  ^ dE (M)  ) 


(3.55) 

(3.56) 


1 .2  -*-2 
2 “E  rC 


(3.57) 


is  the  system  of  forces  due  to  the  so-called  "added  masses".  The  system 

ST  ' " “ I VR  (M)  dE  (M)  + “e  ^C-1  (3.58) 

is  the  quasi-steady  system  of  forces  that  is  the  one  obtained  by  neglecting 
the  acceleration  of  the  body. 

The  reason  why  equations  (3.48)  - (3.49)  may  no  longer  hold  is  that  y 

E 

may  be  rotational.  To  maintain  at  rest  the  liquid  located  inside  , it  is 
necessary  to  add  inside  a force  ?'  per  unit  mass,  so  that 


Since 


V A (F*  - y£)  H 0. 


- + 2, 

YE  =*Ti-+  *(2  °E  r >’ 


(3.59) 


condition  (3.59)  implies  the  existence  of  a potential  ip  such  that 


This  gives 


F'  - y = V 

yE  3t 


pd  (Mi * t)  m P |f  t)  + | P r2(M,  t)  on  E^ 


whence 


i*  ,M  . / i*  x 

3t  ^Mi’  t3  * 3t  V , t* 

e 

One  may  select  for  ip  a harmonic  function  and  write 


(3.60) 


(3.61) 


ip  (M,  t)  - E u.  (t)  ip,  (M)  inside  D, 
j j j i 
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The  ij^'s  are  solutions  of  the  Dirichlet  interior  problems 


*J(Mi)  = - *j(Me),  (j  = 1,  ...  6) 


(3.62) 


Equations  (3.59)  - (3.63)  are  not  of  great  interest  and  are  given  for  the 
sake  of  completeness  only.  The  equations  of  importance  are  (3.54)  - (3.57). 

Generalization 

Let  us  suppose  that  the  liquid  motion  is  due  not  only  to  the  motion  of 
the  body  alone  but  also  to  other  causes.  The  velocity  generated  by  these 
other  causes  is  irrotational  inside  D.,  and  the  kinematic  solution  is  obtain- 

i 

ed  as  indicated  in  the  section  on  effective  calculation  of  the  vortex  dis- 
tribution. The  two  systems  g and  are  altered.  The  new  system  of  bound 
vortices  is 


V = 01  + V2  + V'2 


with 


)=  (Z,  p'  n), 


°2  f ’ 

V'2  = ^Z,  = (Z,  p"n) 


The  velocity  in  the  part  of  D close  to  Z is,  in  general,  irrotational.* 

e e 

The  velocity  potential  in  this  region  is 

$(M,  t)  = 4>2  (M,  t)  + <^(M,  t)  + 4>q  (M,  t) 


with 


>;<V  t>  + ♦"(Mj.  t)  - + v>Pi 

- *,(Mi,  t)-  $o(M1,  t) 


From  the  Kutte-Joukowski  theorem  we  have 

d?T  “ " \ p + % + + d5:(M) 

with 


(3.63) 


Vr  ’ "J  + ™2  ‘ VE 


32 


The  system  of  forces  S becomes 

q.  s 

Vs  " [dfT]  - «•<*> 

where  is  given  by  Equation  (3.63).  The  system  due  to  the  true  and 
apparent  added  masses  becomes 


S.  = [S'  + S"  ] 

in  in  in 


(3.65) 


where  S|^  is  the  system  S^  of  the  preceding  section  and  S^n  is 


Sin  = [_P  lr  (M*  t)SM  dT’W] 


(3.66) 


S'^  is  the  system  due  to  the  "apparent  added  masses." 

The  above  equations  for  non-steady  motions  imply  the  assumption  that 
no  streamwise  vortices  are  shed  by  the  body.  In  the  opposite  case,  if  the 
motion  is  unsteady,  the  incident  velocity  induced  by  these  vortices  varies 
with  time  and  new  bound  vortices  appear  on  the  hull  surface.  This  phenomenon 
will  be  studied  in  Chapter  6. 
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CHAPTER  4:  EXTENSION  OF  THE  LIFTING  SURFACE  THEORY 

TO  BODIES  WITH  FINITE  DISPLACEMENTS 

For  the  sake  of  brevity,  the  present  chapter  considers  mainly 
translational  motion  of  bodies.  The  theory  expounded  in  Chapter  3 
provides  the  necessary  tools  when  the  angular  velocity  varies.  The 
problem  of  small  motions  about  a uniform  translation  will  be  examined  in 
Chapter  6 and  the  effect  of  the  angular  velocity  will  not  be  ignored. 

We  begin  with  thin  wings. 

THIN  WINGS 

As  before,  let  S denote  the  right-handed  system  of  axes  moving  with 
the  body.  The  wing  is  generally  symmetric  with  respect  to  the  plane 
y = 0;  it  is  generated  by  profiles  P( y) , and  its  transverse  sections  are 
generated  by  the  planes  y = constant,  with  - Y < y Y.  The  length  t(y) 
of  P(y)  is  the  difference  X(B)  - X(A)  of  the  abscissas  of  the  ends  B and 
A of  the  profile.  The  base  chord  L of  the  wing  is  the  length  of  P(0). 

The  wing  moves  in  the  negative  x-direction.  When  y varies  from  -Y  to  +Y, 

> the  lines  described  by  points  A and  B are  termed  leading  edge  and  trailing 

edge,  respectively.  The  edges  are  not  necessarily  straight  lines  parallel 
to  the  y-axis,  and,  for  y = + Y,  points  A and  B do  not  necessarily  coincide. 
Let  z+(x,  y)  and  z (x,  y)  be  the  third  coordinates  of  a point  respectively 
describing  the  upper  arc  AB  and  the  lower  arc  AB  of  the  profile  P(y). 

The  chord  of  the  profile  P( y)  is  the  length  i of  a straight  line  joining 
A to  B.  The  chords  of  the  profiles  are  not  necessarily  parallel  to  one 
another.  The  curvature  of  a profile  P( y)  is  generally  very  large  at  B, 
and  so  this  point  can  be  considered  an  edge  in  the  geometrical  sense  of 
the  word.  In  contrast,  the  profile  is  rounded  at  A.  In  all  the  cases, 
the  surface  I of  the  wing  is  defined  by  the  equations  z = z+(x,  y)  and 
z = z (x,  y).  The  aspect  ratio  of  the  wing  is  the  ratio  Y/L;  it  is 
rarely  less  than  5 or  6 and  can  be  much  greater. 

The  angle  of  attack  is  the  angle  (AB,  Ox)  and  it  is  generally  small. 

In  this  case,  the  streamlines  of  the  relative  motion  leave  the  wing  from 
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the  trailing  edge  and  there  Is  no  separation  except  whMi  the  curvature 
of  the  leading  edge  is  too  large. 

Although  the  motion  is  three-dimensional,  it  is  close  to  two-dimensional 
in  the  vicinity  of  the  plane  y = 0 when  the  aspect  ratio  is  sufficiently 
large.  Two-dimensional  motions  can  be  realized  in  wind  tunnels.  In  this 
case,  the  wing  is  the  part  of  a cylindrical  surface  that  is  parallel  to  the 
y-axis  and  intersects  the  walls  y = + Y of  the  wind  tunnel.  Experiments 
show  that  the  circulation  F of  the  velocity  in  closed  circuits  surrounding 
the  wing  is  not  zero  and  that  the  relative  streamlines  leave  the  profiles 
at  B.  These  two  facts  are  interconnected.  They  constituted  an  enigma  for 
many  years  because  they  seemed  to  contradict  a consequence  of  the  Helmholtz 
theorem  according  to  which  the  vorticity  in  a fluid  motion  is  zero  when 
the  fluid  starts  from  rest  without  any  shock.  This  question  was  elucidated 
by  Prandtl  who  demonstrated  experimentally  that  a free  vortex  with  an 

dr 

intensity  equal  and  opposite  to  — dt  appears  at  the  trailing  edge  and  is 

carried  by  the  flow  to  infinity  downstream  from  the  wing.* 

When  the  aspect  ratio  of  the  wing  is  finite,  a free  vortex  sheet  is 
shed  from  the  trailing  edge  8.  Let  E denote  the  surface  of  this  free 
vortex  sheet.  When  the  relative  motion  is  steady,  the  free  vortex  sheet 
is  unbounded  in  the  positive  x-direction.  It  is  bounded  ahead  by  6 
and  its  edges  are  the  streamlines  coming  from  the  two  lateral  edges  of  the 

-►  | _ — ► -f 

wing.  The  continuity  of  the  pressure  through  E entails  V - V = nF~T, 

t R R t i 

— ► -4* 

where  n^.  is  the  unit  vector  normal  to  E^  in  the  direction  from  E^. 

toward  Ep  here  E*  and  E^  are  the  upper  and  lower  sides  of  , respectively. 

— ► 

The  velocity  V of  the  fluid  points  with  respect  to  a fixed  system  of 

axes  S is  small,  and  T,  is  approximately  parallel  to  - V . In  fact, 

t E 

the  angle  (Ox,  V ) is  slightly  greater  than  the  angle  - (Ox,  V ),  and 
R R 

at  a large  distance  from  the  wing,  the  free  vortex  filaments  L tend  to 

wind  around  the  free  tip  vortices  coming  from  the  lateral  edges  of  the 


♦Reported  to  the  Third  Congress  of  Applied  Mechanics,  Stockholm 
(1931). 
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wing.  The  vortex  filament  L ^ is  the  free  part  of  vortex  filaments 

L » L + L the  bound  arcs  L,  being  located  on  E. 

f b b + + 

Tb  + 

When  the  wing  is  very  thin,  the  bound  vortices  - — on  E and 


4 b * "f  — 

— on  Z at  points  (x,  y,  z (x,  y))  and  (x,  y,  z (x,  y))  can  be 


replaced  by  a unique  vortex  - — located  on  the  skeleton  E^  of  the 
wing.  In  this  case,  E^  is  considered  as  a vortex  sheet.  The  tangential 
component  of  the  velocity  is  discontinuous  through  E^.  We  can  write 


V(-M)  - V4>(M),  4>(M) 


&JL 


« jjt  hr  hp  dE(p> 

VZf  P 


It  is  generally  assumed  that  E^  is  the  surface  generated  by  half-straight 

lines  whose  origins  are  on  the  trailing  edge  8 and  which  are  parallel  to  the 
x-axis.  The  function  with  which  p coincides  on  E^  is  constant  on 

such  a generatrix.  Inside  E , the  lines  p,  - constant  are  orthogonal 

0 b 

to  the  relative  streamlines.  Let  (n,  T)  be  a right-handed  system  of 

three  unit  vectors;  n is  normal  to  E^,  in  the  direction  from  E* 

- “¥  “►  ■> 
toward  E and  T is  in  the  direction  of  T . Then  V is  in  the  direction 

U D K 


of  5.  In  the  linearized  theory. 


i 6p  » i <p"  - p+) 


9p  ^ 

lvEl  3^  . (with  do  = 3 do) 
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and  the  lift  coefficient  is 


nr/' 

0 E JC 


ds 


is  the  area  of  Eg  and  C the  contour  of  the  wing. 

boundary  conditions  are 
< 

V*n  = V *n  on  E- 
E 0 

< p = constant  on  the  leading  edge 

= 0 on  the  trailing  edge 

The  third  condition  states  that  there  is  no  pressure  jump  on  through 

the  trailing  edge.  This  is  the  expression  for  the  Kutta  condition  in 

three-dimensional  motions.  It  is  known  that  the  above  conditions 

determine  p,  on  Ln.  We  will  verify  this  point  in  Chapter  5. 
b 0 

The  leading  edge  is  difficult  to  determine  because  the  relative 
velocity  is  very  high  in  this  region  when  the  shock-free  condition  is 
not  fulfilled  (that  is,  when  the  relative  velocity  is  not  tangent  to  E^ 
along  the  leading  edge).  The  theory  of  matched  asymptotic  expansions 
is  generally  used  to  circumvent  the  difficulty. 

We  conclude  this  portion  with  a short  discussion  of  the  theory  of 
thin  wings;  note  that  in  this  theory  the  bound  vortex  filaments  on  the 
skeleton  are  all  open.  All  the  other  cases  that  will  be  examined  later 
have  a subclass  of  vortex  filaments  closed  on  the  surface  of  the  body. 


where  Eg 
The 
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WINGS  WITH  A FINITE  THICKNESS 

The  free  vortex  filaments  Lf  are  arcs  of  vortex  filaments  L = LbU  L{, 

L.  being  the  bound  arc  of  L.  Consequently  the  vortex  distribution  V,  which 
b 

is  kinematically  equivalent  to  the  set  consisting  of  the  wing  and  the 
free  vortex  sheet,  can  be  written  in  the  form 


The  vortex  distribution  V ^ generates  the  velocity 


V,  = 0 inside  D 
1 e 

curl  V = 2^„  within  D 

1 Ci  1 


(4.1) 


(4.2) 


As  in  Chapter  3,  we  have 


V^M,  t)  = X uJ(t)V1J(M)  within  n*VL  = 0 on  Z± 


(4.3) 


and  there  exists  a velocity  potential  defined  within  such  that 


-*  -►  ^<J)  * _ 

V„  - V,  - V*'  with  -5-  - n*V  on  Z 

El  dn  E i 


<t> ' (M,  t)  - X u (t)*'(M)  within  Dj 

j J J 


(4.4) 
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On  the  other  hand,  the  vortex  distribution  V ^ is  equivalent  to  a 
normal  doublet  distribution  on  (E  + E^): 


t>2  = (E,  y'n)  + (Ef,  Ufnf) 


Here  n is  normal  to  E and  n^.  to  E^.  Let  us  put 


/ *' 


V »f 


4tt 


JJr<p- 


'>  hr  sf  dE(p> 

P 


i?//“f<p-  hr  hf  d£(p> 


The  total  velocity  induced  by  V is 


/ Vl  + V,  inside  D 
2 f e 


\ V1  + v2  + vf  inside 


Let  us  determine  p'  by  the  condition 


*»«*«_* 
z f 


(4.5) 


(4.6) 


(4.7) 


(4.8) 


within  D^.  We  obtain 


v = vl  + <v*'  - v*f)  + v*f  - vx  + Vg  - vx  - vE 

within  D^. 


(4.9) 
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Since  the  normal  derivatives  of  and  are  continuous  through  Z,  and 

-f  - ► 

as  n*V  = n'^g  on  the  boundary  condition 

n*  V = n*  V on  Z (4. 10) 

t e 


is  satisfied.  To  determine  on  £»  it:  is  necessary  to  know  £f  (M,  t). 

This  implies  that  Zf(t)  and  Mf(P,  t)  are  known  functions  of  t.  Before 

examining  this  problem,  it  is  suitable  to  study  in  some  detail  the 
characteristics  of  the  vortex  sheet  in  the  neighborhood  of  the  trailing 
edge. 


CHARACTERISTICS  OF  THE  FLOW  IN  THE  NEIGHBORHOOD 
OF  THE  TRAILING  EDGE  (See  Figures  4 and  5.) 

The  trailing  edge  8 divides  the  region  of  the  wing  close  to  it 

into  two  parts,  Z^  which  is  included  in  the  pressure  side  and  Z^  which 

is  included  in  the  suction  side.  Let  Z be  the  side  of  Z^_  which 

f 2 f 

corresponds  to  Z^  and  let  Z^  be  the  side  corresponding  to  Z^.  The 

unit  vector  n normal  to  Z is  in  the  inward  direction.  On  Z^,  the 

unit  vector  nf  normal  to  Zc  is  directed  from  Z,  to  Z,  . 

t t f 2 f 1 

Consider  three  vortex  filaments  L L^,  L located  on  Z^,  Z and 

Z^,  respectively,  which  intersect  the  trailing  edge  8 at  the  same  point 

B;  Mf,  , and  M2  are  three  points  close  to  B on  L^,  L^,  L , 

respectively  (see  Figure  4).  Furthermore  consider  three  pairs  of 
points,  namely. 


(M£ 


V;  <*V  V*  «i,'  V 


defined  by 
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defined  by 


MfMf, 

= nf(0+)  ; 

MfMf2 

= - nf(0+) 



“►  -f- 

. 

->  + 

MlMll 

= <V°  > 5 

M,M 

1 ei 

II 

1 

3? 

O 

¥ 

->  + 

M2Mi2 

MM 
2 ez 

II 

1 

aP 

KJ 

O 

Let  C be  the  closed  circuit 

C:  M Mc  M.  M M M M 

e2  i2  ti  ej  ii  I2  e2 

Let  n be  the  order  of  magnitude  of  the  distance  from  B of  the  points 
Mj,  M^,  and  M^.  The  circulation  T (C)  of  V in  C is 

T(C)  •=  [4>(M  ) - <t>  (M . )]  + [<5>(M,  ) - 4>(M,  )]  + 4>(M,  ) - $ (M  ) + 0(p)  = 0 

e2  i2  ti  12  ii  ei 

Hence,  since  p|  and  represent  the  determination  of  p'  on  and 
respectively,  we  have: 

P ' 2 ~ Pf  “ u'j  = 0 or  p = p'2  - p'x  at  B (4.11) 

Let  Lj,  L|,  L'2  be  three  other  vortex  filaments  on  £f,  £^,  and 
intersecting  8 at  the  same  point  B'.  Let  Mj,  M^,  and  Mj  be  three  points 
on  the  relative  streamlines  containing  M^,  M^,  and  M^,  respectively. 

We  introduce  the  following  four  closed  circuits: 
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M M'  Ml  M.  M 


1 

ei 

ei 

fi 

fi 

ei 

C2  : 

M' 

e2 

M 

e2 

Mf2 

Mi2 

M' 

e2 

Ci  = 

Mi, 

Ml 

ii 

m;2 

Mi2 

Mii 

Of  : 

Mft 

Mf, 

Mi2 

Mf2 

Mf, 

• m-2) 

(Mi, 

* ^ 

) 

and 

(m12* 

from  M|. , Mj,  Ml,  as  the  first  three  pairs  are  derived  from  M^.,  M^,  M^. 

The  flux  of  through  is  zero  and  also  through  and  through 

since  no  vortex  filament  intersects  open  surfaces  bounded  by  these 
contours'.  Hence  being  the  circulation  of  the  velocity 


V2  " V2  + Vf  due  t0  V2 


■(E  + Efr)  • 


we  have 


r2(c1)  = o , r2(c2)  = o , r2(ct)  = o 


(4.12) 


Through  C^,  the  flux  of  T^  is 


r2<cf)  ■ 'w  - 


(4.13) 


From  Equations  (4.12)  we  derive: 


v2(m 

v2(m 

v2(m 


e l 


e2 


ii 


)*M^M|  + V2(M 

)-m]m2  + v2(m 

+ V2(M 


fi 


f? 


i2 


)-M^Mf 


)*MfM’ 


)-m^m2 


0(h) 

0(h) 

o(h) 
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By  adding  the  first  two  equations  and  subtracting  the  third  one,  we 
obtain: 


[V2(Me1)  ' W1  M1M1  + tV2(Me2}  " VMi2)]  M2M2  = r2(Cf)  + 0(n) 


Within  D^,  we  have 


v = V'  + V = V'  = V - V 
2 2 f E 1 


The  last  equation  thus  gives: 


v (m  )-m1m;  + v_(m  ) • m 'm  + r (c  ) = r (cf)  + o(n) 

Rej  11  Re2  22  li  2t 


(4.14) 


where  r,(C.)  is  the  circulation  of  V,  in  C..  Since 
1 i 1 l 


i«v  ■ II  2iV5  dsi 


(S^  is  an  open  surface  whose  edge  is  and  v is  the  unit  vector  normal 
to  S^),  we  see  by  taking  for  a surface  close  to  that 


T1(C1)  = O(irBB') 


This  shows  that 


VV’Vi  ‘ ’ W + °<"> 


(A. 15) 


Let  n approach  zero.  This  means  either  that  the  limiting  positions 
Bj,  B^ , Bf  of  , M^,  Mf  coincide  with  B or  that  the  limiting  positions 
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B^,  B^,  B^  of  Mj,  M^,  coincide  with  B'.  In  fact,  B is  the  limiting 

position  of  M^,  , Mf  and  B'  that  of  Mj,  M^,  M^..  To  prove  this,  it 

suffices  to  consider  the  case  when  , M^,  are  at  B.  If  B^  is  not  at  B', 

it  is  somewhere  between  B'  and  M^,  and  BB^  is  orthogonal  to  L ' at 

B^.  Since  the  free  vortices  move  with  the  fluid,  the  fluid  point  P', 

which  is  at  Bj.  at  the  present  time  t,  was  at  B'  at  a certain  time 

t - dt.  Consequently,  the  relative  velocities  at  B^  on  E and  at 

B'  on  E-  , which  are  orthogonal  to  L ' , are  parallel  to  B'B.  Similarly 
12  t2  t 

the  relative  velocity  at  B'  , which  is  orthogonal  to  L is  parallel 

e2 

to  B'B  and  the  relative  velocity  at  B'  , which  is  orthogonal  to  1' 

2 ei  X 

is  parallel  to  B'B.  Since  BB'  is  infinitely  small,  V (B  ) is  parallel 
1 K 6l 

to  BB!  and  V_(B  ) is  parallel  to  BB'  and  both  V (Bf  ) and  V (B.  ) are 
1 K 62  K 1 2 K 1 1 

parallel  to  BB].  But  is  continuous  from  E to  E and  from 
r f R ei  fi 

I to  E . Consequently  V (B  ) = V (B  ) and  V (B  ) = V (Bf  ). 
e2  f 2 R ei  R fi  R e2  R f2 

It  follows  that  BBj  and  BB^  are  parallel  to  BBj  and  that 

BBj^  = BB^  * BB^.  By  comparing  Equations  (4.13)  and  (4.15)  and  taking 

into  account  the  continuity  of  V„  from  E to  E and  from  E to 

J R ei  fi  e2 

E,  , we  could  have  obtained: 

12 

Vd(B  ).(Bll  - BB,')  - Vd(B  )*(BB]  - BB])  = 0 (4.16) 

R ei  r i R e2  t i 

But  the  above  reasoning  proves  more  since  it  gives 

BB^  = BBJ  = BB],  (4.17) 
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Furthermore  we  see  that 


the  trailing  edge  8 is  orthogonal  to  Lp  and  L 2 (4.18) 

Let  us  suppose  that  T (C  ) > 0.  Then  ft  is  in  the  direction  of  BB' 

t K, 

and  ft  , in  the  direction  of  B ' A . Consequently,  the  intensities 
K2 

dr^,  dr  2,  and  of  the  vortex  ribbons  L^,  L^,  bounded  by 
(L^,  L^),  (L2,  L2'),  and  (i.f,  Lp  are  connected  to  one  another  by 

dFf  = dri  + dF  2 (4.19) 

This  equation  evidently  holds  if  T (Cp  < 0.  It  follows  that  the  two 
bound  vortex  filaments  L , 1.  converge  at  B if  Tf  is  in  the  outward 
direction  with  respect  to  the  wing  and  that  they  start  from  B if  T^. 
converges  to  B.  This  also  shows  that  and  L2  cannot  be  considered  as 
included  in  the  same  vortex  ribbon. 


other  and  (see  Figure  5) 

dri  = d^  = ^ drf  , (steady  relative  motion)  (4.20) 


The  properties  indicated  by  Equations  (4.18)  and  (4.19)  do  not  seem  to 
have  been  remarked  on  previously.  They  are  the  extension  to  the  unsteady 
case  for  which  results  have  been  published  by  Roy.^ 

However,  the  above  scheme  implies  a contradiction.  For,  if  the 
relative  velocity  on  Z£  and  Z£  in  the  neighborhood  of  B were  really 

parallel  to  B,  the  fluid  could  not  leave  the  surface  of  the  wing. 


^Roy,  M. , "Theorie  des  alles  sustentatr ices  et  des  helices," 
Gauthier  Villars,  Paris  (1934). 
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The  bound  vortex  filaments  L on  E^  and  L ^ on  E 2 have  probably  such 

very  high  curvatures  in  the  neighborhood  of  6 that  at  a short  distance 

from  B,  the  relative  velocities  at  M and  M are  close  to  the  relative 

ei  e2 

velocities  at  M,  and  Mc  , respectively.  But,  clearly,  the  difficulty 
f 1 f 2 

cannot  be  overcome  within  the  framework  of  an  almost  inviscid  fluid 
theory.  Another  aspect  concerns  the  set  of  fluid  points  belonging  to  the 
free  vortex  filament  Lf.  This  set  increases  continuously  with  time.  This 

means  that  fluid  points  which  were  vortex-free  at  t may  belong  to  at 
t + dt.  This  contradicts  the  Helmholtz  theorem.  By  assuming  that  the 
fluid  is  almost  inviscid,  one  must  necessarily  accept  contradictions 
somewhere.  The  contradiction  here  consists  of  (1)  a certain  discontinuity 
of  the  directions  of  V and  (2)  a discontinuity  of  its  intensity  at  the 
origin  and  the  end  of  each  free  vortex  filament.  In  a real  fluid,  the 
boundary  layer  thickness  is  finite  and  the  flow  inside  it  is  three- 
dimensional;  on  the  outer  side  of  the  layer,  the  direction  of  the  flow 
is  not  orthogonal  to  the  vortex  filament  derived  from  the  almost  inviscid 
fluid  theory  and,  furthermore,  new  fluid  particles  continuously  enter  the 
boundary  layer. 

The  above  contradiction  does  not  prevent  the  almost  inviscid  fluid 
theory  from  being  acceptable  in  regions  other  than  those  where  the  effects 
of  viscosity  are  of  considerable  magnitude  (as  in  the  neighborhood  of  the 
trailing  edge).  In  particular,  and  L are  unquestionably  not  parts  of 

a unique  bound  vortex  filament.  Although  L ^ is  entirely  located  on  the 

suction  side,  L starts  on  the  pressure  side,  reaches  the  leading  edge, 

passes  on  to  the  suction  side,  comes  back  to  the  leading  edge,  and 
rejoins  the  trailing  edge  at  the  same  point  where  L ends.  The  parts 

E^  and  E^  of  E covered  by  and  L ^ are  separated  by  a line  joining  the 

two  lateral  edges  of  the  wing  and  located  on  E^.  If  there  are  two  tip 
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vortices,  this  line  belongs  to  the  vortex  filament  which  includes  the 
tip  vortices.  The  common  ends  of  1 ^ and  L ^ are  symmetrical  to  each 

other  with  respect  to  the  transverse  plane  of  symmetry  of  the  wing.  The 
line  which  starts  from  the  midpoint  of  B and  comes  back  to  it  is  the 
frontier  of  a region  £^  covered  with  bound  vortex  filaments  L entirely 

located  on  I . 

The  wing  surface  £ is  thus  divided  into  three  regions  £^,  E^,  and 
E^  (see  Figure  5).  The  lift  follows  from  the  fact  that  the  number  of 

arcs  of  vortex  filaments  located  on  the  suction  side  E^  includes  all  the 

filaments  i.^,  L and  whereas  on  the  pressure  side,  there  are  only  the 

vortex  filaments  and  L^. 


THE  GENERALIZED  KUTTA  CONDITION 

As  indicated  earlier,  this  condition  states  that  the  relative  velocity 

V_  is  continuous  from  £ to  £,  and  from  £ to  E . Let  us  assume 
R ei  fi  e2  f 2 

that  the  curvature  of  £ along  the  trailing  edge  8 is  large  but  finite, 

i.e.,  that  there  exists  at  every  point  B of  B a plane  tangent  to  E.  Let 

(resp.  v ) denote  the  unit  vector  with  its  origin  at  B,  normal  to  b and 

tangent  to  (resp.  1^) . Let  B^  (resp.  B^)  be  a point  infinitely  close  to 

B and  belonging  to  (resp.  £„  ) . Two  bound  vortex  filaments  1^  on  l , and 


!1 


1 


L 2 on  E^  end  or  begin  at  B.  The  components  in  the  v - direction  of  the 
velocities  induced  by  the  vortex  distribution  V ^ and  by  the  velocity 
potentials  and  are,  respectively, 

0 , (5’<bu).71  * ) + (-  v5i  + ^ ) i Vj-Ij 
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Here  p'  and  p"  are  the  parts  of  p^  such  that 


/ - — i ffp'(P)  - — — i— 
4ti  JJry  *•  ' 9n  M P 
Z pi 

l - ta//“"(p)  hr  ftp 

L pi 


dZ  (P)  = <t>'  (M  ) 


dZ(P)  = - $f(Mi) 


The  component  in  the  v -direction  of  the  relative  velocity  at  B is 

1 ei 


[" v' (b  ) - v 1-V  + /iHl  + 

L ( U>  VE  J vi  \k3vi  + 3vx  ) Bx 


Since  V (B  ) must  be  parallel  to  8,  we  must  have 
el 


.*’<V  VK  + (|^  + |^)  B[  - 


(4.21a) 


Similarly,  we  must  have 


0<V-JE]V  + f£)  #2 


= 0 (4.21b) 


-T  -f 

But  V is  continuous.  Furthermore,  and  are  in  opposite  directions. 

r ■+ 

The  above  two  conditions  are  thus  equivalent.  We  may  take  for  v one  of 

the  two  vectors  and  v^. 
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Let  us  come  back  to  Equation  (4.8).  It  may  be  written  in  the  form 


-sr//“ 2mh  cp  61  m 


(M1) 


^(Mj) 


(4.22) 


Let  G(M,  P)  denote  the  resolvent  kernel  of  this  equation.  We  have 


W^M) 


(M,  P)t|/'(P)  dE(P) 


But  E^  is  unknown.  To  determine  a first  approximation  of  the  solution,  we 
have  to  make  assumptions  concerning  the  motion  of  I.  For  instance,  if 
is  a constant  and  = 0,  we  may  consider  E^  as  generated  by  half-straight 
lines  starting  from  8 in  the  direction  opposite  to  V . Taking  the  x-axis 
in  the  direction  of  V and  with  a a curvilinear  abscissa  on  8,  we  have 


$f(P)  - - 
This  can  be  written: 
$f(P)  " " 


dx(Q) 


dx(Q) 


Hence  at  a point  M belonging  to  £,  we  have 
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U"(M)  = 


- JIG  <“• 


P)*f(P)  dE(P) 


(4.23) 


ff 

r°  2 r-  00  1 

G(M,  P)  dZ(P) 

l 

" a¥  J d°'l  Uf(o')H(o',  P,  Q)  dx(Q)  I 

0 •'x  (o  * ) j 

The  generalized  Kutta  condition  therefore  becomes: 


II 


3G_ 


(B,  P)  dZ(P) 


[-  f do'  f U (o')H(a',  P,  Q)  dx(Q)l 

L o.  Jx(o'y  J 


(4.24) 


[5'<V  - V“>]  V II, 


9 G 


av  (B,  P)<F(Pt)  dZ(P) 
X B 


By  interchanging  the  order  of  integrations,  it  may  be  written 
in  the  form 


- f1 


uc(0')  do'  f dx(Q)  f f 3—-  (B,  P ) H (O  ' , P,  Q)  dZ(P) 


;x(Of) 


'z  9vb 


(4.25) 


- [v-ciy  - ve(b)J  ;b  - J{£  (»,  P)  (Pl)  «(p) 


The  equation  that  determines  on  B is  thus  a Fredholm  equation  of 
the  first  kind.  It  is  certainly  very  difficult  to  prove  that  this 
equation  has  a unique  solution,  but  we  may  admit  it  for  physical  reasons. 
In  practice,  the  use  of  the  resolvent  kernel  of  Equation  (4.8)  would 
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be  difficult  and  if  its  numerical  computation  were  necessary,  other  processes 
would  have  to  be  tried. 

UNIFORM  MOTION  OF  TRANSLATION  OF  A SUBMERGED 
BODY  IN  A VERTICAL  PLANE 

Now  we  consider  a body  with  a longitudinal  plane  of  symmetry  zOx.  This 

plane  is  vertical  and  is  also  symmetric  with  respect  to  the  (x,  y)  plane 

which  is  horizontal  when  the  squat  is  zero.  The  above  symmetries  are  not 

needed  by  the  theory,  but  they  help  to  expound  it.  We  assume  that  the 

— ► 

body  is  in  a uniform  motion  of  translation.  The  velocity  V of  the  origin 

b 

of  the  moving  axes  0(x,  y,  z)  is  in  the  (z,  x)  plane.  The  angle  of  attack 

is  6 = (Ox,  V„) . This  is  positive  when  the  incident  velocity  ( - V ) 

E h 

has  a positive  z-component. 

Experiments  show  that  if  6 0,  the  system  of  hydrodynamic  forces 

-f 

exerted  on  the  body  includes  a lift  orthogonal  to  V , even  when  the  shape 

h 

of  the  body  differs  considerably  from  that  of  a wing.  This  lift  is  due 
to  the  shedding  of  vortices.  Otherwise,  all  the  bound  vortices  would 
be  closed,  and  according  to  the  d'Alembert  paradox,  S1*  would  reduce 
to  a couple.  Various  experiments  (for  instance,  those  carried  out  by  Sears 
with  an  axisymmetric  body)  show  that  separation  occurs  at  the  afterbody 
and  that  the  shedding  of  free  vortices  follows.  (Figure  6) 

More  precisely,  let  us  suppose  that  the  forebody  begins  with  a 
stem  parallel  to  the  z-axis  and  that  the  stern  ends  at  a point  S on 
the  x-axis.  (Figure  7)  When  6=0,  all  of  the  relative  streamlines 
start  from  the  stem  and  end  at  S.  We  have  a leading  edge  consisting  of 
the  stem,  no  trailing  edge,  and  zero  lift.  When  6 > 0,  the  streamlines 
start  not  only  from  the  stem  but  also  from  the  lower  part  of  the 
contour  C formed  by  the  intersection  of  the  hull  £ with  the  (z,  x)-plane. 

The  leading  edge  now  consists  of  the  stem  and  the  arc  C^.  Roughly 
speaking,  the  relative  streamlines  coming  from  the  upper  part  of  the  stem 
seem  to  Intersect  those  coming  from  the  lower  part  of  the  stem  and  from 
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.j 


along  a U-shaped  line  8 consisting  of  an  arc  B^S'lT  located  on  the 

afterbody.  The  ends  Bq  and  B^  of  ® are  slightly  dovmstream  from  the 

maximal  transverse  section  of  the  body  and  symmetric  to  each  other  with 

respect  to  the  (z,  x)-plane.  The  point  S'  is  close  to  S. 

The  hull  £ is  thus  divided  into  three  parts,  £ ^,  £ 5 ^ is 

the  upper  part  of  the  afterbody  and  its  boundary  is  the  line  B and  a 

half-ring  C,  containing  B and  B'  ; £ is  the  lower  part  of  the  afterbody 
z o o 1 

and  its  boundary  is  the  line  8 and  a half-ring  the  part  ^ is 

ahead  of  the  ring  + 'C  (see  Figure  7). 

There  are  three  kinds  of  vortices.  The  vortex  filaments  are 
entirely  bound.  Their  supports  are  rings  on  £^.  The  vortex  filaments 
L ^ come  from  infinity  downstream  and  reach  8 at  points  B located  on 
the  portside  (Figure  8);  on  £^  their  supports  are  half-rings  starting 
from  points  B and  ending  at  points  B*  which  are  symmetrical  to  points  B 
and  located  on  8.  They  leave  B from  B'  and  extend  back  to  infinity 
downstream.  The  vortex  filaments  are  analogous  to  the  but  their 

bound  parts  are  located  on  A free  vortex  is  the  union  of  the 

free  parts  of  one  L ^ and  of  one  L If  d'  ^ and  dT^  are  the  intensities 
of  and  L^,  the  intensity  dT  of  L ^ is 

dFf  = dFl  + dT2  (4.26) 

The  scheme  just  defined  is  analogous  to  that  for  a wing  of  finite 
thickness  except,  of  course,  that  the  curvature  of  the  hull  along  the  line 
8 is  not  necessarily  high.  If  the  fluid  were  really  "almost  inviscid," 
then  not  only  the  bound  arcs  of  and  L but  also  the  free  arcs  L ^ 

should  be  orthogonal  to  8.  This  would  require  a considerable  curvature 

of  the  Lf’s  in  the  vicinity  of  8,  and,  in  fact,  is  impossible.  Similarly, 

-+-*■_  _ 

the  relative  velocity  V„  and  V„  on  £ and  £ cannot  be  in  opposite 

Ri  R2  ei  e2 

directions  along  8.  The  almost-inviscid  fluid  theory  leads  to  contradictions 
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which  are  similar  to  those  encountered  in  the  case  of  the  wings  and 
perhaps  is  still  open  to  criticism. 

However,  the  only  point  of  practical  interest  here  is  the  assumption 
that  free  U-shaped  vortices  really  do  exist. 

One  may  remark  that  the  above  scheme  is  not  basically  different 
from  that  considered  by  the  slender  body  theory.  To  show  this,  consider 
a transverse  section  of  the  afterbody.  (Figure  6).  For  the  vertical  force 
to  be  other  than  zero,  we  have  to  assume  in  the  slender  body  theory  that 
the  flow  is  separated  as  seen  in  Figure  8,  which  implies  the  existence 
of  a pair  of  longitudinal  free  vortices  piercing  the  plane  of  the 
figure  at  T,  T'.  These  free  vortices  necessarily  originate  on  the  body 
and  their  intensity  grows  in  the  negative  x-direction.  The  slender 
body  scheme  thus  seems  to  be  derived  from  that  proposed  here  by  assuming 
that  the  length  of  the  body  is  large  enough  for  the  free  vortices 
originating  ahead  of  the  transverse  section  under  consideration  to  have 
already  wound  up  around  the  two  "tip  vortices." 

My  personal  feeling  is  that  because  the  real  phenomenon  is  so 
complicated,  the  theory  is  able  to  provide  us  only  with  a usable  frame 
to  determine  the  general  form  of  the  equations.  The  number  of  the 
coefficients  found  in  the  equations  should  be  as  small  as  possible  and 
their  numerical  values  should  be  derived  experimentally. 

Another  drawback  is  that  the  theory  does  not  accurately  tell  us 
the  position  of  the  shedding  line  on  the  hull.  But  this  should  not 
entail  the  risk  of  large  errors  in  the  derivation  of  the  system  of 
forces  exerted  on  the  body.  The  generalized  Kutta  condition  is  expressed 
by  the  equation  applicable  to  wings  with  a finite  thickness. 

DOUBLE  MODEL  IN  AN  OBLIQUE,  UNIFORM  MOTION 
OF  TRANSLATION  IN  THE  HORIZONTAL  P1ANE 

Now  consider  a body  whose  shape  is  that  of  a double  model  relating 
to  a surface  ship  (Figures  9 and  10).  The  (x,y)-plane  is  horizontal,  the 
longitudinal  plane  is  vertical.  These  two  planes  are  planes  of  symmetry 
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for  the  double  model.  The  velocity  V„(0)  of  the  origin  of  the  moving 

axes  is  in  the  horizontal  plane  and  the  drift  angle  6 * (V  , Ox)  is 

fc. 

positive.  The  lower  half  of  the  distribution  of  singularities 
kinematically  equivalent  to  the  body  should  be  considered  in  the  case 
of  the  so-called  zero-Froude  number  approximation  r'"r  a s •''face  ship.  It 
is  known  that  this  approximation  very  rough  for  calcuj . _ng  the 
waves  generated  by  a ship.  But  we  will  see  later  (Chapter  5)  how  to 
correct  it.  Thus  the  mathematical  model  described  in  the  present  section 
can  be  considered  as  a first  approach  to  the  maneuvering  theory  of  a 
surface  ship. 

The  z-axis  is  vertical  upward.  Let  E^  denote  the  portside  (y  ^ 0) 
and  E^  the  starboardside  (y  < 0).  Let  C be  the  contour  intersected 
in  the  hull  Z » by  the  (z,  y)-plane,  Cp  its  lower  part, 

its  upper  part,  Cp  the  arc  EE^S^S,  and  the  arc  SS^E^E.  The  arc 
EqSq  is  a horizontal  straight  line  (keel).  The  arc  E^EEp  of  the  double 
model  is  its  stem,  and  the  arc  SqSS^,  its  sternpost.  On  C,  the  selected 
sense  is  that  from  Oz  to  Ox.  The  curvilinear  abscissa  is  A,  and  its 
origin  is  taken  at  E on  the  x-axis.  We  put 

A(Eq)  - XQ,  X(SQ)  - \'Qt  A(S)  = + A 

»<Sj)  - x;  - - xj,  ME,)  - x,  - - x0 

The  drift  6 is  small  enough  for  no  separation  to  occur.  Experiments 
on  surface  ship  models  show  that  the  relative  streamlines  are  more 
inclined  on  the  pressure  side  (starboard)  than  on  the  suction  side 
(port).  Thus,  the  two  relative  streamlines  which  reach  the  same  point 
Bq  on  the  arc  EpSpS  do  not  have  identical  directions  at  Bp.  This 
determines  a shock  and,  as  a consequence,  a free  vortex  filament  is 
shed  from  Bp.  The  origins  on  the  stem  of  the  two  streamlines  are  not 
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at  the  same  level,  the  origin  of  the  streamline  located  on  the  pressure 
side  being  higher  than  that  on  the  suction  side. 

Since  we  are  de-  ' with  a double  model,  the  same  phenomenon  occurs 
on  the  upper  halt  che  body.  The  arc  playini  the  role  of  trailing 
edge  is  the  part  ^S^SSE^  of  the  contour  C,  an<-  the  stem  E^EE^  plays 
the  role  of  leading  edge.  However,  if  the  bottom  is  flat,  then  the 
after  part  of  the  shedding  line  tends  to  move  toward  port  when  6 increases. 
Beyond  a critical  value  of  6,  separation  oc«.  »s.  What  follows  does  not 
apply  in  this  case. 

Let  Lj  be  a free  vortex  filament  which  leaves  arc  SS^E^  from  a point 
B^.  It  goes  to  infinity  downstream  and  comes  back  to  the  body  at  point 
Bg,  the  mirror  image  of  B^  with  respect  to  the  plane  z = 0 and  thus 
located  on  the  arc  SS^E^.  By  studying  unsteady  motions  of  the  double 
model,  we  will  see  later  that  the  bound  part  of  the  vortex  filament  L 
to  which  belongs  consists  of  a unique  arc  B^B^  necessarily  located  on 
the  suction  side.  The  phenomenon  therefore  differs  from  that  for  wings  of 
finite  thickness  since,  in  the  present  case,  L is  not  divided  on  Z into 

two  arcs  L.  , L located  on  Z , Z„,  respectively. 

Di  D2  1 i 

If  6 were  zero,  the  bound  vortex  family  would  consist  of  vortex  rings 
on  transverse  sections  of  the  hull.  The  half-ring  to  port  would  be  in 
the  direction  from  B^  to  B^,  and  the  half-ring  to  starboard  would  be  in 
the  direction  from  B^  to  B^.  When  6 is  positive,  then  to  port,  we  have 
a vortex  filament  of  intensity  dP^  from  B^  to  B^,  and  to  starboard,  we 
have  a vortex  filament  of  intensity  dr£  from  B^  to  B^.  The  intensity 
dT^  is  greater  than  the  intensity  dT  ^ since  L is  located  to  port  and  in 
the  direction  from  B^  to  B^.  The  two  arcs  B^B^  to  port  and  B^B^  to 
starboard  are  no  longer  symmetric  with  each  other.  Let  Cj  denote  the 
arc  BqB^  to  port  and  C^,  the  arc  B^Bq  to  starboard.  These  two  arcs 
depend  on  the  abscissa  A of  B^  on  E^S^S.  Ahead  of  the  arc 
Cj(AQ)  + ^(Aq),  the  bound  vortex  filaments  on  Z are  vortex  rings,  but 
the  parts  of  each  ring  located  to  port  and  starboard  are  not  symmetric 
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to  each  other.  We  can  consider  that  behind  the  ring  C|(A^)  + 
we  are  dealing  with  two  distributions  of  bound  vortex  filaments.  One 
consists  of  vortex  rings  with  intensity  dF^CA);  the  other  consists  of 
bound  vortex  filaments  of  intensity  dF ^ (A ) located  on  the  part  Ej(A^) 
of  the  portside,  which  is  downstream  from  C|(A^). 

The  vortex  rings  generate  a velocity  potential 


$’(M) 


_1 

4 TT 


_1 

MP 


d£  (P) 


(4.27) 


On  the  other  hand,  let  E be  the  free  vortex  sheet  and  E , E , E'  the 

t to  r 1 t 

parts  of  Ej  on  which  we  find  the  free  vortex  filaments  reaching  the 
longitudinal  contour  between  E^  and  S^,  those  which  leave  C from  E^S^, 
and  those  which  leave  C or  reach  C at  points  located  on  S^SS^.  The 
vortex  filament  consisting  of  a free  vortex  filament  L^_  and  of  a bound 
vortex  filament  L"  located  on  the  part  E (A  ) of  the  portside  is 

D 1 U 

equivalent  to  a normal  dipole  distribution  on  the  part  of  E^.  + E^ 
downstream  from  its  contour. 

It  follows  that  the  part  of  vortex  distribution  V ^ due  to  all  the 

vortex  filaments  L"  + L is  equivalent  to  a normal  doublet  distribution 
b t 

on  E^  + E^: 


*£(M) 


JL 

MP 


dE2(P) 


(4.28) 


JL 

47T 


1 

MP 


dEf(P) 


The  sum  may  be  written  in  the  form 
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4>2(M)  + <t> f (M)  , with 

*2<m>  - - JL',<p>  It  sf  dZ<p>  <‘-29> 

Z P 

VM)  - - i 1 “f<p)  lr  mp  dVp> 

Zf  p 


The  function  u'  defined  on  Z has  three  different  determinations 
according  to  whether  P is  located  on  Z^A.^),  on  °r  °n  the  part 

Z^  of  Z ahead  of  C|(Aq)  + C^C^q).  It  is  determined  by  the  condition 
that 

4>2(m)  “ - ^(M)  + Vj.(0)  (x  cos  6 - y sin  6)  within  (4.30) 

This  gives: 

-i  *■'<">-  zj  IJ/'<p)lr  sf  dE(p)- 

l p 

(4.31) 


- (M^)  + Vg(0)  [x  cos  6 - y sin  6]  for  M e Z 

The  three  determinations  of  U'  are  distinct  from  one  another 
because  of  the  antisymmetry  between  the  two  sides  of  Z (due  to  the 
antisymmetry  of  Z^  with  respect  to  the  double  model)  and  also  because 
no  free  vortex  filament  is  shed,  from  the  arc  EqEE^.  On  Z^,  is  a 
constant  along  the  free  vortex  filament  L^.  Hence,  depends  only 
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on  the  function  — — on  the  arc  E.S  SS  E of  C.  Let  us  suppose  that  this 
dA  l)  U I 1 

function  is  known.  Then  p'  is  determined  by  Equation  (4.31).  The  arcs 
C^(A)  and  C9(A)  are  the  lines  p'  = constant  and  the  relative  streamlines  on 
E are  the  curves  C orthogonal  to  these  lines.  To  determine  the  flow  entirely, 
a complementary  condition  is  needed.  This  is  the  generalized  Kutta  condition. 

THE  GENERALIZED  KUTTA  CONDITION  FOR  BODIES 
IN  AN  OBLIQUE,  UNIFORM  TRANSLATION  IN  THE 
HORIZONTAL  PLANE 

Let  Bq  be  a point  on  the  arc  E^S^S  of  the  longitudinal  contour  C of 

the  double  model  and  let  fl^  be  the  plane  tangent  to  £ at  (Figure  11). 

Let  Bq(Vq’  *'»  Zq)  denote  a right-handed  system  of  axes,  B^A'  being  tangent 

to  C in  the  direction  of  increasing  A's.  The  axis  B^y^  is  in  the  positive 

y-direction.  The  axis  Bqzq  is  not  necessarily  vertical  upward  but  its 

direction  coincides  with  the  inward  normal  n at  B^.  We  consider  *~he  right- 

handed  systems  of  unit  vectors  (n,  0^,  x^)  and  (n,  O^,  x^)  for  the  port  and 

starboard  sides,  respectively.  We  recall  that  the  directions  of  V_. , V . 

K1  Rz 
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are  opposite  to  0^  and  0^.  respectively.  In  the  above  two  sketches,  Figure 

11,  we  have  taken  into  account  the  fact  that  the  angle  | (V  , B_A')|  is  less 

K 1 U 

than  the  angle  |(VR2>  BQX  * ) | - 

2 2 

In  the  neighborhood  of  B^,  zQ  = 0(y ^ + A'  ).  Hence,  and  r can 

be  considered  in  that  region  as  functions  of  (y^.  A').  Since  F ^ is 

constant  along  C|(A)  and  ^(A)  is  constant  along  C'^(X) , the  equations 
of  the  tangents  at  B^  to  C|(A)  and  C^(A)  are  as  follows: 


ar. 


3I\ 


yJ  + jxl  »'  ■ ° . yjio 


ar  ar 

a^  y'0  + wr  x'  =°  . y'Q i° 


(4.32) 


The  equations  of  the  tangents  to  the  streamlines  t ^ and  ending  at 


Bo  are 


ar. 


3A ' 'o 
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3A*  *0 


yA  - 


ar 

a^J  A’  " 0 ■ ^o^0) 


ar2 

a^  v " 0 • (yo±0) 


(4.33) 


Note  that 


ar.  ar. 


ar„  ar. 


3A ' 


9yo 


< 0 


3A' 
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> 0 
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Let  Mx,  M2  denote  the  intersections  of  (p  with  C|(A  - dA ) and 
of  with  (A  - dA ) . We  have 


do 


1 


M1B0 


Ri 


This  gives: 


dr. 


V <V 


dA'  cos  (BQy^  , T|) 


1 

2-,  2 


dr 

dA 


, [(5)  ♦(&)  ] 


5 

3A’ 


, for  y0  = 0 


But  dr 


5 

1 = 3A  ' 


dA’  on  EqSqS.  Consequently 


vr1(bo) 


and  similarly 


V <V 


w~-  © 

W-  so 


. 34a) 
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The  generalized  Kutta  condition  states  that  the  pressure  on  the  hull  is 
continuous  through  E^S^S  an<i  therefore  that 


x <v 


(4.35) 


This  gives: 


(4.36) 


dr 

For  — — to  be  other  than  negligibly  small,  the  slope  of  Cl(X)  with  respect 
d a 1 

to  bqVq  and  that  °f  C^X)  with  respect  to  B^X ' must  be  relatively  small. 

Thus  is  seems  that 


!Ll  !!i 

ax'  >>  " 3y^  * 3y^  3X' 


(4.37) 


This  gives 


(4.38) 


It  is  difficult  to  get  more  information  unless  numerical  computation  is 
attempted. 
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FINAL  FORMULAS  FOR  BODIES  IN  AN  OBLIQUE, 
UNIFORM  TRANSLATION  IN  THE  HORIZONTAL  PLANE 


Let  P be  a point  on  Z and  let  P . , P be  the  points  defined  by  PP,  * n 

1 e ip 


(0+)  and  PP  ■ - n (OF),  respectively.  Let  do  be  the  element  of  arc  of  the 
e p 

relative  streamline  which  passes  through  P on  I or  that  of  the  relative 

® ® 3 ij 1 3 r 

streamline  determined  on  I by  the  vortex  distribution.  We  have  *■  — 

3a  3 a 

where  u*  is  the  density  of  the  normal  doublet  distribution  in  (4.31)  and 
3T  -> 

— t is  the  intensity  of  the  vortex  T (P)  on  E.  Thus  we  have 


y ' (P)  - r(P)  - $2(Pe)  - *2(Pi)  (4.39) 

Here  F is  defined  up  to  an  additive  constant;  it  increases  from  the  stem  to 
the  stern  post,  with  r-T^,  T 2 or  F^  according  to  whether  P is  on  Z^,  Z2>  or 
Zy  On  Zj  ' Z^  + E^q  + E^,  we  have  at  any  point: 


"f<V  ' “f<V 


(4.40) 


where  B^  is  the  point  located  on  E^SS^  which  belongs  to  the  free  vortex 
filament  containing  P^. 

The  unit  vector  normal  to  Z^  at  P^.  is  selected  so  that  it  has  a direc- 
tion coherent  with  that  of  n on  Z.  Thus  it  is  in  the  upward  direction  on 
Z,  and  in  the  downward  direction  on  Z,  ; on  Z',  it  is  in  the  direction 

fl  Jo  + f 

from  port  to  starboard.  The  points  P^  and  P^  are  defined  by 

PfP*  - nf (0+)  , PfPj  = -nf(Ob) 

v 

We  easily  see  that 

M f (P)  - *f(P")  - *f(p+)  - Ff(BQ)  (4.41) 

Finally,  we  may  write*. 
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(4.42) 


4>(M)  - 4’2(M)  + *f(M)  with 


• 2(M) 


SJl 


r(p)  s!;® dE(p) 


+ E_  + E, 


(4.43) 


T is  determined  by  the  condition: 


- \ r(M)  - ||  r(P)  dl(P)  - - *f  (M)  + VE(0)  [x  cos  6 - y sin  6] 


MeE, 
(4.44) 


SaW  s(SBo  • x(Bo)e[V 


A] 


Furthermore,  we  have  in  the  first  iteration: 


P f Bq  * [x(BQ)  - x(Pf)]  i(VE)  on  Ef 


(4.45) 


where  the  difference  in  square  brackets  indicates  the  difference  between  the 

abscissas  of  B,.  and  P measured  in  the  direction  of  V . 

0 t, 

The  second  formula.  Equation  (4.44)  is  the  reproduction  of  Equation 
(4.38).  Since  the  latter  is  an  approximate  equation,  used  only  because  at 
the  first  stage  r is  unknown,  it  would  be  convenient  after  this  first  compu- 
tation, to  check  whether  it  is  necessary  to  replace  (4.38)  by  (4.36)  (the 
values  of  and  in  (4.36)  are  those  given  by  the  first  iteration).  This 
would  lead  to  the  second  iteration. 
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CHAPTER  5:  PROBLEMS  RELATING  TO  THIN  SHIPS  AND  TO  FREE  SURFACE  EFFECTS 


It  happens  very  often  in  ship  hydrodynamics  that  drastic  simplifications 
of  a problem  yield  interim  solutions  which  later  become  useful  bases  for  a 
more  refined  treatment  of  the  real  case.  This  has  been  so  for  thin  ship 
theory  in  the  case  of  wave  resistancfe.  For  this  reason,  the  present  chapter 
attacks  the  problem  of  the  maneuverability  of  a very  thin  ship  and  examines 
how  to  take  into  account  the  free  surface  effect  when  the  boundary  condition 
on  the  free  surface  can  be  linearized. 


TENTATIVE  THEORY  FOR  INFINITELY  THIN  SHIPS  IN  AN  OBLIQUE,  UNIFORM 
TRANSLATION  (STEADY  CASE) 

General  Comments 

In  a dissertation  prepared  under  my  direction  about  20  years  ago, 

Casal  outlined  a theory  of  small  aspect  ratio  wings  for  its  applicability 
to  ship  maneuverability.  His  results  are  of  interest  from  many  points  of 
view.  The  picture  he  gave  of  the  hydrodynamic  forces  exerted  on  a maneuver- 
ing ship  is,  in  general,  in  rather  good  qualitative  agreement  with  experi- 
ments. However,  his  assumptions  were  rather  crude,  and  one  wonders  whether 
their  refinement  could  lead  to  numerical  calculations  in  a relatively  simple 
manner. 

As  mentioned  previously,  the  hypothesis  that  the  ship  is  very  thin  has 
not  been  discarded,  but  assumptions  concerning  the  bound  vortex  distribution 
on  the  longitudinal  plane  of  symmetry  of  the  hull  have  been  eliminated. 

This  was  done  in  order  for  the  boundary  condition  on  the  hull  to  be  satisfied 
everywhere  and  not  merely  on  the  longitudinal  axis  of  the  longitudinal  plane. 
Of  course,  no  real  ship  is  thin.  But,  the  addition  of  a source  distribution 
on  the  longitudinal  plane  might  enable  the  treatment  of  bodies  with  shapes 
relatively  close  to  those  of  real  ships.  Indeed,  the  very  thin  wing  gives 
rise  to  difficulties  for  the  leading  edge.  (This  had  been  mentioned  at  the 
beginning  of  Chapter  4.)  But  thanks  to  the  additional  source  distribution, 
it  would  probably  be  possible  to  remove  them. 

It  will  be  seen  below  that  the  greatest  difficulty  originates  in 
combining  the  thinness  of  the  model  with  the  very  small  aspect  ratio  of  a 
real  ship  form.  The  latter  is  sensibly  less  than  unity  but,  in  the  case  of 
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wings,  it  is  usually  much  greater  than  unity.  The  shedding  line  begins  at 
the  junction  of  the  stem  with  a keel;  thus,  the  phenomena  occurring  in  the 
neighborhood  of  the  lateral  edges  of  a normal  wing  are  considerably  magni- 
fied. The  numerical  calculation  of  the  solution  would  probably  be  more 
intricate  than  that  of  the  solution  sketched  in  Chapter  4 for  real  ship 
forms.  For  this  reason,  I was  almost  at  the  point  of  withdrawing  the  sub- 
ject of  the  very  thin  ship  from  the  present  lectures.  I retained  it, 
however,  because  it  may  throw  light  on  the  nature  of  the  problem. 

Notations  and  Assumptions 

The  double  model  of  a thin  ship  is  reduced  to  a vertical  flat  plate  in 
the  (z,  x)-plane.  The  horizontal  plane  (Ox,  Oy)  is  that  of  the  waterline 
of  the  surface  ship.  The  axes  0(x,  y,  z)  move  with  the  body.  Also,  a 
second  system  of  axes  0 (x^,  y^,  z)  is  introduced  for  convenience.  The 
planes  (Ox^,  Oyp  and  (Ox,  Oy)  coincide;  the  x^-axis  is  in  the  direction  of 
the  constant  velocity  V of  the  common  origin  of  the  two  systems  of  axes. 

L 

The  contour  C of  the  plate  consists  of  the  stem  EE^  of  the  surface 
ship,  of  its  keel  EqSq,  of  its  stern  post  S^S,  and  of  the  other  three  arcs 
EEj,  E^S^,  SjS  which  are  the  mirror  images  of  the  former  with  respect  to 
the  (x,y)-plane.  The  x-axis  is  in  the  direction  from  S to  E.  The  z-axis 
is  vertical  upward.  It  intersects  the  keel  in  the  midpoint  of  the  segment 
EqSq.  Let  2Z  be  the  length  of  this  segment.  The  keel  E^S^  is  horizontal. 
The  drift  angle 


(OXj , Ox) 


is  supposed  to  be  positive  and  small  enough  for  no  separation  to  occur  on 
the  real  ship.  No  free  vortex  is  shed  from  the  stem  E^EE^  of  the  double 
model.  The  shedding  line  is  the  arc  E^S^SS^E^.  The  supports  of  the  free 
vortices  are  assumed  to  be  in  the  negative  x ^-direction.  A free  vortex 
filament  L.  consists  of  two  arcs,  L.  and  L,  or  L\,  and  L',..  L,.  is  in 

ft  fti  6o  ft1  ft0  ft1 

the  plane  z **  H;  it  leaves  the  arc  S.E,  at  a certain  point  B, ; L,  is  in  the 

11  1 fto 

plane  z - -H  and  reaches  the  contour  C at  a point  B^,  mirror  image  of  B^. 

1^  generates  the  part  of  the  free  vortex  sheet  1^,  while  L ^ generates 
the  part  The  free  vortex  filament  leaves  the  upper  half  of  the  stern 
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post  at  a certain  point  Bj , and  L.  ends  on  the  lower  half  of  the  stern  post 

l u 0 

at  B',  mirror  image  of  B!  . L\  and  L',  generate  the  part  Z\  of  Z..  A free 

u f 01  00  00 

vortex  filament  L ^ (or  L is  closed  by  a bound  vortex  filament  L ^ which 
joins  Bq  to  Bj  (or  B^  to  Bj). 

The  contour  C is  oriented  in  the  positive  direction  (from  Oz  to  Ox). 

Let  ds  denote  an  element  of  arc  of  and  let  v be  the  unit  vector  normal 

to  C in  the  outward  direction.  Hence,  if  ds>0,  (ds,  v)  = Let  x=-X(z) 

Z T 

-V  V, 

be  the  equation  of  the  stern  post  in  the  plane  y = 0.  Let  w,  = — be  the  bound 

b e 

vortex  at  a point  M of  the  area  Z of  the  longitudinal  plane  of  the  hull. 

-f 

As  e is  a constant  length  equal  to  0+,  T is  divergenceless  on  Z.  Thus, 

b 

there  exists  a function  ip  of  the  variables  (z,  x)  defined  on  Z such  that 


Tfe(M)  = cur  | (i^(M)iy) 


(5.1) 


where  iy  is  the  unit  vector  in  the  positive  y-direction.  The  components  of 
T^  on  the  axes  0(x,  y,  z)  are  denoted 


_3jf) 


. n * 0 , z. 


(5.2) 


The  flux  of  T,  through  an  arc  ds  of  C in  the  outward  direction  is 

D 


dr  = T.  • v ds 

D 


(5.3) 


This  quantity  is  the  common  intensity  of  the  free  vortex  filament  L ^ and  of 
the  bound  vortex  filament  L ^ intersecting  the  arc  eqSqS  at  Bq  or  Bq  and  the 
arc  E^S^S  at  or  The  function  ip  is  an  odd  function  of  z and  dip/dz  is 
even.  According  to  the  above  assumption  concerning  the  stem,  ip  is  constant 
on  EjEEq. 

For  convenience,  we  note 


x = the  abscissa  of  on  E„S. 
o 0 0 0 


-X  (z  ) = the  abscissa  of  B'  on  S„S 

o o 0 0 


(5.4) 
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The  velocity  V induced  at  M by  the  total  vortex  distribution  is  the  sum  of 
the  contributions  V(M/E,  ) from  I,  , V(M/l,  ) from  1.  , V(M/Z,)  from  1 , 

6o  60  fii  oi  6 o 

and  V(M/l)  from  E.  It  will  be  calculated  by  using  the  Biot  and  Savart 

-> 

formula.  Before  performing  this  calculation,  let  us  remark  that  T,  is  not 

b 

tangent  to  C (except  on  E^EE^).  Therefore  L ^ cannot  be  in  the  direction 
parallel  to  V on  E_S  , E.S.,  and  S_SS  . This  contradicts  our  assumption 
concerning  the  geometry  of  the  free  vortex  filaments.  The  assumption  will 
be  slightly  modified  when  dealing  with  the  generalized  Kutta  condition. 


Velocity  Induced  by  the  Free  and  the  Bound 
Vortex  Sheets 

It  is  easier  to  use  the  system  of  axes  0(x^,  y^,  z). 
(xj,  y j , zp  or  (x',  y',  z')  on  We  have 


Let  M'  be  at 


“i(M/E<  + E<,  + Ei> 


(5.5) 


* 1 

-i  -°° 


x„cos  6 


-H-z 


1 


dx’;  (yj=xn  sin  6,  zp- H,  z0=-H) 


MM 


, 1 3 r wi  o 


J !^%  I ° 


x„cos  6 


-SL 


H-z, 


MM 


— - dxp  (yj=xQ  sin  6,  zj=H,  zQ= H) 


(5.6) 


',«/e<>-^  J (%)% ! 0 

-H 


,-X^cos  6 


21~Z1 

I MM'  I 


3 dxi;  (yi=-XQ  cos  6’  zi=zo*  xo=~V 
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■kj 


x„coa  6 . . 

0 xn  sin  6-y 

— o— 1 - dx’ 

Imm'I3  1 


,x„cos  6 


xQ  sin  6-y1 


(5.7) 


-X  cos  6 


-XQ  cos  6-y 
~ dx: 

Imm’I3  1 


Behavior  of  v^(M/Z^).  Let  us  suppose  that 


< » on  S-E  and  ® on  S_S 

3xq  0 0 3zq  0 


(5.8) 


Then  v^  is  bounded  on  I and  on  the  contour  of  E (except  at  and  if 
* 0 at  these  two  points). 


Behavior  of  w^(M/E^).  w^  is  also  bounded  on  Z.  But  on  S^E^  and  S^E^, 

|w|»*  because  of  the  assumption  that  L.  is  parallel  to  V at  its  two  ends.* 

6 E 

Velocity  Induced  by  the  Bound  Vortex  Sheet.  By  using  the  axes  0(x,  y,  z)> 


■I'M  ' iT 


<M/E>  ‘ 47  JJ£  IT  - t»’>  lr]  dE("’> 

“<m/e>--4T  Jj 


(5.9) 


E<M'>  It  s’- dE(M'> 


*See  the  section  dealing  with  boundary  conditions  on  the  hull  and  the  one 
on  the  generalized  Kutta  condition  below. 
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The  components  u and  w may  be  regarded  as  the  normal  derivatives  on  Z 
of  potentials  due  to  sources  with  respective  densities  -£(u')  and  £(M'). 
Let  Z+  and  E be  the  surfaces  derived  from  E by  translation  ± 0 in  the 
direction  of  Oy.  Then 


u+(M/Z)  - u“(M/Z)  = -£ (M) 
w+  (m/ E ) - w~(M/Z)  = +C(M) 


As  C > 0, 


u+(M/Z)  - u“(M/Z)  > 0 


The  absolute  value  of  the  x-component  of  the  relative  velocity  is  larger  on 
the  suction  side  (port)  than  on  the  pressure  side  (starboard).  Furthermore, 
the  z-component  of  V (M/Z),  which  is  negative  for  z = -H,  is  less  in 
absolute  value  to  port  than  to  starboard  if  ( > 0 on  S0E0* 

This  shows  that  our  mathematical  model  agrees  with  the  cause  of 
shedding  along  SqEq  (and  S^E^)  and  that  the  pressure  is  effectively 
smaller  to  port  than  to  starboard.  The  y-component  of  V (M/Z)  is  obviously 
continuous  through  Z.  It  may  be  written  in  the  form 
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that  is. 


v(M/T) 


1_ 

An 


1 

MM' 


A<KM')  dZ(M') 


(5.10) 


In  this  formula, 


The  Boundary  Condition  on  the  Hull 
This  condition  is 

iy  • V (M)  - - C sin  6 for  y = 0 (5.11) 


or 


(M/E ^)  cos  6 + v(M/E)  - - C sin  6 on  £ .5.12) 

According  to  the  previous  section,  this  gives: 


i_  ff_!_ 

An  JJ  MM' 


Ai(/(M')  dE(M)  = c sin  6 + 


5?  J+  F(M-M 


')  |±  (M*)  d4(M') 


h J+  (!!)„,  bp-  <HeE> 


(5.13) 


with 


«■»  "•  EiEEo  <5-“> 

F(M,M')  on  the  right-hand  side  of  Equation  (5.13)  is  a known  function  of  M 
and  M' , which  also  depends  on  the  drift  6. 
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The  Generalized  Kutta  Condition  for  an  Infinitely 
Thin  Double  Model 

For  the  sake  of  brevity,  let  us  put 


ci  ‘ WVi ' c - EiEEo 


(5.15) 


Let  M be  a point  inside  the  bound  vortex  sheet  and  M , M the  points  derived 
from  M by  MM  = iy(0+),  MM  = - iy(0+).  The  Kutta  condition  states  that 


p(M  ) * p(M  ) if  M e C 


1 


(5.16) 


It  is  equivalent  to  the  condition  that 


VrCM)  • [Vr(M+)  - Vr(M  )]  = 0 on  Cl 


(5.17) 


Since  w(M/£)  = 0,  u(M/I)  = 0,  we  have 


Vr(M)  - w(M/5^)  ±z  + [u(M/E ^)  - c cos  6] 


Hence  (5.17)  may  be  written  in  the  form: 


i|/ (M)  w(M/Z^)  + ip(M)  [u(M/E^)  - c cos  6]  - 0 on  (5.18) 


We  have  seen  that  the  assumption  that  the  free  vortices  reach  C or  leave  C 

in  a direction  parallel  to  V entails 

E 


|w(M/l^)|  - « on  S1E1  and  E^ 


Thus,  this  assumption  as  to  be  modified 

We  now  assume  that  if  M^  is  inside  the  free  vortex  sheet,  then  |v(M^)j 
| V (M) | and  VR(M^)  is  parallel  to  V£  (provided  that  |MM^|>  0).  This  gives 


Vr(M^)  - [u(M/Z^)  - c cos  6]  + w(M/2^)  ty 
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1 


Inasmuch  as  this  vector  is  parallel  to  V^,,  its  y ^-component  is  zero  and 

[u(M/Z^>  - c cos  a]  sin  6 + w(M/£^)  cos  6=0,  |MM^|  > 0 (5.19) 

By  substituting  this  into  Equation  (5.18),  we  obtain 


iHM)  tan  6 • 'I'OO  on  SqEq> 

' (5.20) 

~ i^(M)  tan  6 = ~ ^(M)  on  S^E^ 

i 


At  S,  w(M/£  ,)  = 0 for  reason  of  symmetry,  and  w(M/£^)  is  probably  small 
S„S 


on 


JoJr 

(5.18): 


Wc  thus  assume  that  w(M/£ ^) 


0 on  This  gives,  after  Equation 


\HM)  • 0 on  SQSS1 


(5.21) 


and  we  recall  that 


0 on  Ej^EEq 


(5.22) 


Formula  (5.20)  shows  that  the  z-component  of  Tfe  is  much  less  than  its  x- 
component  in  the  vicinity  of  EQS0  and  S^.  Hence  the  curvature  of  the 
vortex  filaments  is  generally  high  on  the  plate. 

The  Integral  Equation  of  the  Problem 

The  general  solution  of  the  equation 


A\J)  = f (M)  inside  £ 


(5.23) 


is  <|»  = ijij  + i|>q  with 
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^JJ 


log  I MM' I f(M')  d£(M'), 


(5.24) 


A^0(M')  = O’  inside  Z 


Our  purpose  is  to  show  that  Equation  (5.13)  allows  iJj  to  be  expressed  in 
terms  of  the  harmonic  function  i^,  and  that  the  generalized  Kutta  condition 
determines  ij/n. 

To  that  end,  let  us  consider  the  operators  F^,  F^  which  define  the 
Fourier  transform  of  a given  function  inside  Z and  its  inverse.  We  have 


/ £(0,k)  = F 


= -i—  exp  [-ik(z’  cos  0 + x'  sin  9)]<?(M')  dZ(H) 

Z 


F^1{^(0,k)} 


if  - r 


(5.25) 


exp  [ik(z  cos  0 + x sin  0)]  ^(0,k)  kdk 


V(M) 

if  MeZ 

, 0 

if  M $ I + C 

Let  us  put 

Y (M)  - c sin  6 if  M e Z,  = 0 if  M ( T + C 


*0(M>  ■ h |+  F("-M,)  (tit) ds‘“'>  + It  J+  sF(ir2)M.<is<“') 


> (5.26) 


4(M)  -h  J+  f(m'm,)(Sl)m,  ds(M’>  J+^(^r)M,  ds(M,) 


In  the  sense  of  the  generalized  function  theory,  we  have 

.it 

d0  I exp  (ik[(z-z')  cos  0 + (x-x')  sin  0]}  dk  (5.27) 


1,1  r 

MM'  ’ 2ir  J 


-IT 
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By  substituting  this  expression  into  (5.13),  we  obtain 


feff 


exp  [ik(z  cos  0 + x sin  0)]  kdk 


-tt  0 


till 


exp  [-ik(z'  cos  0 + y'  sin  0)]  f(M')  d£' 
right-hand  side  of  (5.13). 


Hence : 


whence 


4-  F.  (f  )i 


1 2k  1 


right-hand  side  of  (5.13), 


f(M)  =F  j" 


-1 


{2kF^  {right-hand  side  of  (5.13)}} 
{2ky(M)  + 2kgQ(M)  + 2kg  j (M) } 


From  this  we  derive: 


log  |mm' I 


log  |mm' 


log  |MM' 


f(M')  dZ(M') 

{F"1  [ 2kg  L ] } dl(M') 

{F"1  [2k  (y  + g0(M))]}M,  dZ(M') 


(5.28) 


Equation  (5.28)  is  a Fredholm  equation  of  the  second  kind  which  gives  iji  in 
terms  of  and  of  c sin  6 inside  £.  Let  G denote  the  resolvent  kernel  of 
this  equation.  We  obtain: 
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^(M) 

H-(M) 

Hq(M) 


[H^(M')  + Hq(M')]  dKM')  with 
IMM'I  [F_1{2k-7}]M,  dl(M'), 
|MM'|  [F_1{2kg0)]M,  dl(M’) 


\ 


(5.29) 


The  harmonic  function  ^ is  a linear  functional  of  its  values  on  C.  We  may 
write 


*0(M)  = L{M|^0(C)} 

It  follows  from  this  that 


(5.30) 


Hq( M)  = 1'{M|^0(C)}  , 


(5.31) 


where  L'  is  also  a linear  functional  of  ^(C).  The  generalized  Kutta 
condition  on  C and  the  condition  which  indicates  that  ^ = constant  on  E^EE^ 
may  be  written  in  the  form: 


a(s)  + b(s)  f*  = 0 on_C 


By  taking  into  account  Equations  (5. 29)- (5 . 32) , we  obtain: 
a(s)  1{M|*0(C)}  + b(s)  L{M  1 4<0  (C)  } 
i(s)  jj  G(M.M')  L'(M'!^0(C)}  dL(M') 


+ a( 


+ b (s) 


h jj  c<"-" 


')  L'{M'|*0(C)}  dl 


- [a(s)  — + b (s) 


fejj 


G(M.M')  H^(M')  dl'(M')  on  C 


(5.32) 


(5.33) 
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This  is  an  integral  equation  of  the  second  kind  which  determines  i|/  (C)  . 

When  ij/g  is  determined  on  C,  is  determined  by  Equation  (5.28)  and 

* - *0  + ^ (5.34) 

The  preceding  proves  the  existence  of  the  solution  of  the  problem.*  But  the 
writer  dares  not  claim  that  the  above  method  should  be  used  for  computing 
the  solution. 

The  problem  of  the  infinitely  thin  double  model  does  not  appear  to  be 
simpler  than  that  of  the  double  model  of  a real  ship.  Despite  its  short- 
comings, the  Casal  theory  is  probably  the  best  one  to  use  when  the  purpose 
is  simply  to  give  a qualitative  account  of  the  phenomena  involved  in  ship 
maneuvering. 

A TENTATIVE  THEORY  FOR  NON- INF I NIT ELY  THIN  SURFACE  SHIPS 

We  shall  consider  only  the  case  of  a ship  in  uniform  motion  of  transla- 
tion and  assume  that  the  relative  motion  is  steady.  Furthermore,  we  assume 
that  the  boundary  condition  on  the  free  surface  may  be  linearized. 

Let  Zq  denote  the  wetted  surface  of  the  ship.  Since  the  boundary 
condition  on  the  free  surface  is  linearized,  Eq  is  the  part  of  the  hull 
located  below  the  plane  z = 0 of  the  free  surface  at  rest. 

The  ship  moves  with  speed  c (=V  ) in  the  positive  x -direction.  For 

b 1 

the  sake  of  simplicity,  the  hull  heel  is  neglected.  The  longitudinal  plane 
of  symmet-y  of  ZQ  is  thus  vertical.  Let  (Oz,  Ox)  be  this  plane.  The  drift 


*The  "proof"  given  here  is  incomplete,  however,  since  it  has  not  been  estab- 
lished that  the  determinant  of  Equation  (5.28)  is  other  than  zero. 
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6 = (Oxi , Ox)  is  positive  and  small  enough  for  no  separation  to  occur. 

Vector  Potentials 

To  apply  the  method  used  in  Chapter  4 (Section  entitled  "double  model 
in  an  oblique,  uniform  motion  of  translation"),  we  shall  consider  below  a 
double  model  1 which  lower  half  is  the  real  hull  Zq,  the  upper  half  being 
the  image  Z^  of  Zq  with  respect  to  the  plane  z = 0 of  the  free  surface  at 
rest.  Let  ZfQ  be  the  real  free  vortex  sheet  and  Zf^  its  image.  In  what 
follows,  Z^  denotes  the  union  Zf^UZ^. 

Any  free  vortex  filament  Lf^  located  on  ZfQ  comes  from  infinity  down- 
stream, and  ends  at  a point  Bq  located  on  the  arc  EqSqS  of  the  contour  of 
the  centerplane  of  the  ship.  Its  support  coincides  with  a streamline  of  the 
relative  motion,  but  it  is  in  a direction  opposite  to  the  relative  velocity 
VR.  This  free  vortex  filament  Lf^  is  a part  of  a U-shaped  one  l"  = Lf0+L;-+Lfl. 
The  arc  is  bound  on  Z , the  portside  of  Z;  beginning  at  Bq,  it  intersects 
the  arc  C+  of  the  waterline  C and  ends  at  B^ , the  image  of  Bq  on  the  arc  SS^E^. 
The  third  part  Lf^  is  the  image  of  Lfg  on  Zf^,  image  of  ZfQ.  The  intensity  of 
L"  is  a constant  depending  on  the  abscissa  of  Bq  on  the  contour  of  the  center- 
plane. 

There  exists  on  Z another  family  of  vortex  filaments  L'.  Each  filament 
L'  is,  in  fact,  a vortex  ring.  On  Z , the  starboardside  of  Z,  it  is  in  the 
direction  from  B^eSS^E^E  towards  BqESSqEqE,  B^  being  still  the  image  of  Bq. 

On  Z+,  L'  is  in  the  direction  from  Bq  towards  B^.  The  intensity  of  any  vortex 
ring  is  a constant  depending  only  on  the  position  of  Bq. 

We  have  thus  two  bound  vortex  sheets  on  Z.  One  consists  of  the  union  of 
the  vortex  rings;  it  extends  on  the  whole  surface  Z.  The  second  one  is  located 
on  Z+,  behind  a certain  arc  EqE^.  Let  e be  the  uniform,  infinitely  small 
thickness  of  the  vortex  sheets.  The  vorticity  at  a point  Pq  on  Zq  or  ZfQ  is 

“<po>  ■ 1 V "ith  *0  ■ <!0'  "o'  50>  <5'35> 

At  P^,  image  of  Pq,  we  have 

u(pi)  “ I fj.  With  il  - (-£q,  -n0,  50)  (5.35’) 
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The  two  vectors 


J'(M) 


?^(P) 

PM 


dZ  (P) 


(5.36) 


The  total  velocity  induced  by  the  vortex  distribution  is 

V2  = V’  + V"  (5.39) 

It  is  defined  in  the  whole  space,  above  the  plane  z = 0 as  well  as 

below  this  plane  and  inside  £ as  well  as  outside  E.  But  it  is  only  inside 

D O (z  < 0)  that  has  a physical  meaning,  D being  the  domain  exterior 
ex  e 

to  l. 

Each  of  the  two  vortex  sheets  defined  by  (5.36)  and  (5.37),  respectively, 
is  equivalent  to  a normal  dipole  distribution.  We  may  write: 


V’(M)  - V *'(M),  *’(M)  - — //rw£( P)  ^ d£ (P ) (5. AO) 

V"(M)  - V *"(M),  *"(M)  - ^-//£  u JJ  (P)  3^-  p^-  d£+  (P) 

“ 4tT  3n^  PM  dEf(P)  (5.A1) 

But  we  may  also  write: 

V2  = V'  + V^.,  where  (5.A2) 

h ■ 7V  ‘f(M)  - - K "lf  “£(P)  PS  dEf <p>  <5-43> 
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(5.44) 


% ' '*2-  ‘i(M>  ‘ -4?"l  "(P)^  "SdI  <P> 


U (P)  = y'  (P)  + y^  (P)  on  E+  (port  side) 


U(P)  = < 


(5.45) 


y (P)  *>  y'  (P)  on  E (starboard  side) 


One  has 


y"  = 0 ahead  of  a certain  arc  E„E,  on  E+. 
b 0 1 


As  in  Chapter  4,  let  X be  the  curvilinear  abscissa  along  the  contour  of 
the  centerplane  of  Z,  with  X =*  0 at  E,  X < 0 on  SS^EjE,  X > 0 on  EEqSqS.  We 
have 


dyf 

dr" 

d . + — 

dX 

dX 

" dT" 

= dX  (y  - y ) at  Bo  °n  E0S0S 

(5.46) 


If  the  surface  E^  were  known  and  also,  on  Ef,  the  lines  y^  = constant  (that 
is  the  vortex  filaments  L,)  and  the  function  (XB„),  then  could  be 

I dA  U i 

calculated  on  the  inner  side  E^  of  E.  By  means  of  equation 

($'2  + = c x^(Mj)  + constant  (5.47) 

one  could  determine  y on  E.  The  potential  so  defined  is  the  solution  of 
a Dirichlet  interior  problem.  The  Fredholm  equation  giving  y is  regular  and 
has  a unique  solution.  The  velocity  (equation  (5.42)),  satisfies  the 
boundary  condition  at  the  hull,  namely 

V . n = n . c on  E (5.48) 

2 e 

and  according  to  (5.47),  V . n = n . c on  E.  since  — — *1,  are 

t-  1 dn  l dn  1 

continuous  across  E. 

In  fact,  the  surface  E^  and  the  lines  are  not  known.  But,  having 
regard  to  the  smallness  of  6,  one  may  assume  that  they  are  half-straight 
lines  parallel  to  c (in  the  opposite  direction  on  ZfQ  and  in  the  direction 
of  c on  Efj).  One  may  also  assume  that  they  start  from  the  contour  of  the 
centerplane  (except  from  E^EEj).  For  larger  values  of  6,  the  line  8 would 
be  displaced  towards  portside  on  the  stern  region  and  separation  could  occur, 
particularly  in  the  case  of  a flat  bottom. 
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drf 

The  condition  determining  -jt—  is  the  generalized  Kutta  condition.  It 

dA 

expresses  in  the  present  case,  as  mentioned  in  Chapter  4,  the  equality 

K vi  • iv  <vi on  Vos  (5-19) 

The  exact  treatment  of  this  condition  is  difficult  and  the  approximate 
solution  given  in  Chapter  4 is  rather  crude. 


Wave  field 

Since  we  are  now  dealing  with  a surface  ship,  the  real  velocity  in  the 
liquid  domain  cannot  be  reduced  to  just  in  equation  (5.42).  The  total 
velocity  V must  satisfy  the  free-surface  condition  (here  we  assume  that  it 
can  be  linearized)  and  have  on  the  hull  Zq  and  the  new  surface  of  the  vortex 
sheet  Zf^  properties  analogous  to  those  of 

In  fact,  the  free  vortex  filaments  coincide  with  arcs  of  relative  motion 

streamlines.  Therefore,  they  exert  on  the  flow  a zero  system  of  forces  and 

contribute  nothing  to  the  waves.  (This  would  hold  even  if  the  free-surface 

condition  were  not  linearized).  The  vortices  pearce  the  free  surface  along 

its  intersection  with  the  hull.  To  close  them,  the  free-surface  condition 

being  linearized,  one  may  still  consider  the  images  Z^  of  Z^  and  Zf^  of  Zf^, 

and  the  ^ expression  in  terms  of  y and  also  that  of  in  terms  of  y^  are 

both  unchanged.  But  we  have  to  add  to  <P ^ (or  to  + a velocity 

potential  4>  such  that  the  free-surface  condition  be  fulfilled.  Putting  k = 

2 

g/c  , where  g = the  acceleration  of  gravity,  and  introducing  a vanishing 
friction  force  - e'  V per  unit  mass  of  liquid,  we  must  have: 


(k  — - — — 

o 3z  c 3x, 


3x 


2>  (*2  + 
1 


V 'z  - 0 = 


(5.50) 


Wc  have  already  remarked  that  the  z-component  of  the  velocity  induced  by 
the  vorticity  is  null  at  z = 0 (because  the  free-surface  condition  is 


linearized).  We  note,  furthermore,  that  - — 4>P 

OXl  L 


0 for  z - 0 if  the  L^'s 


are  parallel  to  c.  But  for  reasons  already  mentioned,  <t>j.  contributes  zero 
to  the  x, -component  of  V.  Finally  (5.50)  is  reduced  to 


80 


f,  9 c 3 , 3 . . 

(k  — — + — -)  <t 

o 32  c 3x.  .2  1 

1 3Xj 

.-1 


3x 


A ' 

2 V2 


0. 


(5.51) 


Let  Fj,  F be  the  operators  of  he  Fourier  transform  of  a function 
of  (x,  y)  or  (xj,  yp.  Let  us  put 


IjCM) 


I*  (M) 


i = exp  (-  ik  (x^  cos  9 + y^  sin  9)) 


(5.52) 


and,  Pq  being  a point  on  Z^,  its  image,  we  denote  their  coordinates 


y[,  - z' 


and,  putting  for  the  sake  of  brevity 


= V*0,  II*  = (P0)  = V <V 

we  have 

ipp  **  exp  (ikto),  where  u = (x^  - xp  cos  9 + (y^  - yp  sin  9 (5.53) 

Let,  along  the  waterline  WL,  in  the  plane  z = -0, 

n be  the  unit  vector  normal  to  Z in  the  inward  direction,  (5.54) 

s the  unit  vector  tangent  to  WL  in  the  positive  direction  with 
respect  to  the  z-axis, 

x one  of  the  two  unit  vectors  tangent  to  and  normal  to  WL. 

We  put 


a = cos  (Oxj,  n),  6 - cos  (OXj , s),  y - cos  (Oxj,  i) 


(5.54') 


For  the  sake  of  brevity  we  shall  often  write  f (6, 
Equation  (5.51)  gives 


(9.  k) 


2a 
sec  6 


k(k  - kQ  sec20  + i (0+))  1 


k)  for  Fj  {f  (xL,  y1 , 0) } . 


a2 

).  (5.55) 

3x^  z - 0 
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Let  FS^,  resp.  FSg , be  the  part  of  the  plane  z = -0  interior,  resp.  exterior, 

i 0 

to  the  intersection  WL  of  this  plane  with  the  hull  Zq,  and  , resp.  , 
the  determination  of  ^ within  FS^,  resp.  FS^ . We  have 


‘ *2%'  = " 


\ (any  Mf  6 WL)  (5.56) 


(_jL  — $,e')  = (r  + v 

Udx1  2 3Xj  ?2  'm'  ^ ds  Y 3t;M' 


Furthermore : 


Ff  {^7}  = F1  {^}  = exp  (kz'Hj  (M')/k  for  M £ FS.  + FSe  (5.57) 

(5.57') 


0 


1* 
9 I 


3x^ 


Po  PoM  P1  1 

the  latter  equation  being  due  to 

9 9 


9 1 

9x„  TTFT.for  M £j  FS,  + FS 


3z 


p0 


9zt 


, cos  (0z1  np  ) = - (cos  Oz,  np  ). 


'1 


On  the  other  hand,  because  of  the  discontinuities,  equations  (5.56),  of 

$'  and  across  £,  we  have 
2 9xi 

2 1 

F i { <t> 2 j }=-k2cos2e^+~/  li*((^+Y|^)M.+ik  cosOy  (M' ) )dy  l (M' ) (5.58) 

3Xj  j z=0  WL 


From  (5.57'),  we  derive 

*2  (6’  k)  51  *'l  U(M,)  pTm  + 9^7  PTM) 


Pn  0 


?!  1 z = 0 


dW 


- (after  5.57)  //  - (Pi)  (2  exp(kz')  l!*/k)  dln  (P.)  (5. 59) 

in  4tr  0 9n_  1 0 U 

u PO 

Let  4>j  be  the  contribution  to  from  the  first  term  on  the  right-hand  side 
of  (5.58)  and  42'  that  from  the  second  one.  According  to  (5.55),  we  have 
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p(p0)  9 

*!(e,k)-//  (2  exp(kz')  i;*) 

1 Lo  471  3nP„ 


1 2 

k - k sec  6 + i(0+) 
o 


dEQ(P0)  (5.60) 


S,(e*k)“k;wL"(Bli+Yl^+  lk  cos0u)m\  SBf~~ 


k-k  sec29+i(0+)dnM'Ii*dyl^M'^  (5-61) 


Now,  let  <f>  denote  the  velocity  potential  due  to  a unit  Kelvin  source  located 
at  Po(xl’  yi»  z')»  z'  < 0>  and  moving  with  the  same  speed  c as  the  ship  in 
the  positive  x^-direction.  We  may  write 


♦(M'  V "57  (^M  +i^M  + *i  «•  V> 
57  77m  + 77m  + h CM’  piJJ 


(5.62) 

(5.62’) 


We  have 


♦i  = 


2Ij*  exp(kz') 


k - k sec  6+  i(0+) 
o 


(5.63) 


♦l  = 


2 kQ  sec  01 J*  exp(kz') 
_ 

k(k  - kQ  sec  0 + i(0+)) 


Be  combining  (5.60)  and  (5.63)  we  obtain 


•1  ■ "e„  77  “<po>  77^  *1  dW 


whence 


*2  + *i  " ffZ0  ^ 


p0 


(M,  PQ)  dE  (PQ) 


On  the  other  hand,  by  combining  (5.61)  and  (5.63'),  we  derive: 


(5.63') 


(5.64) 


*1  ■ 77  V HiT  t ItV  + “»’>  7777)  * <"•  "')  <VM'>  (5-65) 


83 


The  total  velocity  potential  (for  z < 0)  is  therefore 


4>(M)  = <tf  + + 4>1  = <t>f  + (*'2  + 4>|)  + t'j1 


4 it  Uf (P)  9np  PM  d5'f 


(5.66) 


+ ffZ0  *(M*  P0)  dW 

u Pq 

+ eV  (-ei7  + lr>M'  +tJ(M')^)  * (M*  M,)  dyi(M,) 

0 1M 

We  recall  that  y^  and  y are  connected  to  each  other  by  the  generalized 
Kutta  condition  and  that  y is  determined  by  the  boundary  condition  at  the 
hull,  namely 


(4>f  + i>'2  + >M  = c Xj  (Mp)  on  Z^ 


(5.67) 


We  also  recall  that  I.  = E.  + Z-  • The  second  term  on  the  right-hand  side 

f fo  fl 

represents  a Kelvin  normal  dipole  distribution  on  the  open  surface  Zq  of 
the  hull.  The  line  integral  represents  the  sum  of  a limiting  Kelvin  source 
distribution  and  of  a limiting  Kelvin  longitudinal  dipole  distribution  along 
the  waterline.  These  last  two  distributions  are  dissymetric. 

We  remind  the  reader  of  the  assumptions: 

- The  intersection  of  I,  with  the  hull  Z.  is  supposed  to  be  in  the 

r 0 U 

centerplane ; 

- The  drift  angle  is  supposed  to  be  small  enough  for  the  flow  to  be 
non  separated; 

- The  heel  is  neglected; 

- The  free-surface  condition  is  linearized. 

There  is  no  need  to  say  that  in'  spite  of  very  important  simplifications, 
solving  equation  (5.67)  would  be  considerably  intricate. 
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CHAPTER  6:  SMALL  MOTIONS  OF  A BODY  ABOUT  A UNIFORM  MOTION  OF  TRANSLATION 


Let  us  first  define  Small  Motions  as  studied  in  the  present  chapter. 

The  fluid  is  unbounded  and  almost  inviscid.  The  speed  V (0)  of  the 

E 

origin  0 of  the  system  S of  axes  0(x^,  x^,  x^)  moving  with  the  body  is, 
in  the  mean,  constant  and  parallel  to  the  x|-axis  of  the  fixed  system  S'  of 
axes  O'  (x! , x',  x' ) . The  components  (u. , u , u ) of  V (0)  on  the  moving 

•Lfcj  ^ b j E 22 

axes  are  such  that,  with  c a positive  constant,  the  quantities  (u  -c)  , u , 

2 2 1 2 
and  u are  negligibly  small  with  respect  to  c . The  angular  velocity  ft  of 
j E 

the  system  S with  respect  to  S'  is  zero  in  the  mean.  Let  L denote  the 

length  of  the  body  in  the  x, -direction  and  let  u . , u , u be  the  components 

-*  * L ^ L^  ®L 

of  ft  on  the  moving  axes.  The  quantities  u.  — , uc  — , u,  — are  small  enough 

4 Z I)  c b ^2 

for  their  squares  to  be  negligibly  small  with  respect  to  c . Free  vortices 

are  shed,  but  the  "shedding  line"  8 is  supposed  to  be  in  a fixed  position 

on  the  hull  l of  the  body. 

The  above  assumptions  define  a small  motion  of  the  body  about  a uniform 

motion  of  translation  of  speed  c in  the  positive  x^-direction . 

Nevertheless  there  are  two  kinds  of  such  small  motions. 

In  motions  of  the  first  kind,  the  line  6 is  approximately  in  a plane 

parallel  to  the  (x  , x )-plane  and  ft  is  approximately  parallel  to  the 
i Z E 

x^-axis.  The  support  of  the  free  vortices  can  be  considered  as  a surface 
generated  by  the  curves  close  to  half-straight  lines  in  the  negative  x|- 
direction  which  will  be  considered  horizontal.  The  small  motions  of  the  first 
kind  are  therefore  analogous  to  small  motions  of  a submarine  in  the  vertical 
plane.  The  rotational  wake  is  confined  on  the  surface  E^. 

In  motions  of  the  second  kind,  the  free  vortices  are  shed  from  a line 
6 analogous  to  the  arc  E^S^SS^E^  of  the  contour  C intersected  in  the  hull 
of  a double  model  by  its  longitudinal  plane  of  symmetry.  Nevertheless,  the 
body  is  not  assumed  to  be  symmetrical  with  respect  to  some  horizontal  plane. 

It  will  be  seen  that  the  motions  of  the  second  kind  differ  from  those  of  the 
first  kind  in  that  the  vortex  wake  generated  by  the  body  is  three-dimensional 
and  thus  not  confined  on  a certain  surface. 
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SMALL  MOTIONS  OF  THE  FIRST  KIND 
Kinematics 

Let  Bg  (to  port)  and  (to  starboard)  be  the  two  ends  of  the  shedding 
line  8.  Let  be  a curvilinear  abscissa  on  8 and  -l,  +S.  the  values  of  o 

at  Bq  and  B^  respectively.  Let  t^  be  the  time  when  the  motion  begins,  let 
8^  be  a line  on  E^  infinitely  close  to  6 (Figure  12),  and  let  P(tg,  o)  be 
the  fluid  point  which  was  at  t^  at  the  point  B^(o). 

At  t this  fluid  point  is  at  a position  B*"  (a)  such  that 

t0 


(P(tQ,  a),  t')  dt' 


When  - with  t fixed  - a varies  from  -i.  to  +£,,  the  line  described  by  the  point 

B^  (o)  is  the  boundary  of  E,(t)  in  the  negative  x' -direction.  Let  8C  denote 
to  f t0 

this  line.  Let  P (t,  o)  be  the  fluid  point  which  was  at  time  t at  B (a) . 

£ 1 
At  t,  this  fluid  point  is  at  a position  B^(a)  such  that 

Bf(o)  B^(a)  = | VR  (P (t , a),  t’)  dt'  (6.2) 

T 

When-with  t and  t fixed-o  varies  from  -i  to  +£,  the  line  described  by  B^(a) 

is  a line  8C  located  on  E,(t)  between  6,  and  B*"  ; 8,.  is  the  line  8 The 

t f f t f t-0 

^ t t 

edges  of  E^(t)  are  the  arcs  described  by  the  points  B^(-£.)  and  B^(+£). 

The  above  definition  of  the  free  vortex  sheet  E^(t)  is  general  and 

involves  no  assumption  concerning  the  nature  of  the  motion  of  the  body 

except  that  for  the  fixity  of  8 on  Z. 

Let  us  consider  now  the  parts  E^  and  Z ^ of  Z which  are  the  supports  at 

t of  the  bound  vortex  filaments  L , L beginning  or  ending  on  8.  The 

bi  b2 

part  Z ^ is  included  in  the  suction  side  and  bounded  by  8 and  by  the  bound 
vortex  filament  beginning  at  B^  and  ending  at  B^.  In  contrast,  E^  is  partly 
on  the  suction  side  and  partly  on  the  pressure  side.  Its  frontier  consists 
of  B on  the  pressure  side,  of  the  bound  vortex  B^  B^  on  the  suction  side. 
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and  of  a certain  vortex  filament  (coming  from  a certain  point  of  B and 
ending  at  the  same  point)  on  both  sides.  The  Z^  part  is  covered  by  closed, 
bound  vortices  only.  The  general  disposition  of  the  bound  vortices  is 
similar  to  that  represented  in  Figure  5 for  wings  of  finite  thickness  in  a 


steady  motion  of  translation.  Any  vortex  filament  whose  bound  part  is 

located  on  Z,  is  continued  on  E,  by  a free  vortex  . The  free  arc  Lr  of 
It  f 1 f l 

L,  ends  at  the  beginning  B of  L and  begins  at  the  end  B'  of  L . The 
1 bx  bi 

point  B is  also  the  end  of  Lc  and  the  beginning  of  L while  B'  is  also 

f2  b2 

the  beginning  of  Lc  and  the  end  of  L . On  I.,  the  supports  of  Lc  and 
1 2 b2  t 1 1 

Lc  coincide  and  the  intensity  of  the  union  Lc  of  1,  and  L,  is  the  sum  of 
f2  f fi  f 2 

the  intensities  of  L,  and  L,  . It  is  suitable  for  the  following  to  consider 

bl  02 

that  L,  and  L are  closed  by  a vortex  filament  located  on  the  arc  B'B  of 
bi  b2 

8 and  that  there  exists  on  B^B^  a filament  of  intensity  doT^(o,  t)  + doT 

(a,  t) , o being  the  abscissa  of  B on  8.  As  B^B^  is  infinitely  close  to  B'B 

and  in  the  direction  opposite  to  BB' , the  effect  of  the  fictitious  vortex 

filament  B'B  is  cancelled  by  that  of  the  vortex  filament  B^B^..  The  advantage 

of  this  procedure  is  that  the  vortex  filament  of  intensity  - do[r^(o,  t)  + 

T^(a,  t)J  located  on  the  arc  B^(o)Bf(o  + do)  is  in  a position  independent 

of  t,  whereas  the  positions  of  L on  I,  and  of  L on  E„  may  vary  with  t. 

Dll  D 2 ^ 

During  the  time  interval  (x,  x + dx),  the  intensity  of  the  vortex 


filament  located  on  the  arc  B^(o)B^(o  + do)  has  an  increment  -d  d [i^Co,  t) 
+ T0(o,  t)]  and  a free  vortex  filament  intensity  +d  d [r,  (o,  t)  + r„(o,  t) 
appears  at  x + dx  on  the  closed  contour  B^(o)  B^X  (o)  BX  T (o+do) 

B^(o  + do)  B^(o).  This  free  vortex  filament  moves  with  the  fluid  and  at  t 
it  becomes 


BX(0K  + dx(0)BxC  + dx(°  + daK<a  + do)B^(o) 


Its  intensity  is  unchanged  and  thus  equal  to  d d [r.(o,  x)  + r.(o,  x)].  The 
t t T O 1 2 

points  B^  + ^(o)  and  B^(o)  are,  of  course,  on  the  line  B^(o)B^(o).  This 

allows  us  to  describe  the  structure  of  the  free  vortex  sheet  Z^(t)  at  time 

t as  follows: 
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Let  (t),  (t)  be  the  two  sides  of  Z (t);  Z (t)  and  Z (t)  are 

fl  f2  4 £1  f2 

the  continuations  at  t of  and  Y.^,  respectively.  Let  n^  denote  the  unit 

vector  normal  to  Zc  in  the  direction  from  Zc  (t)  toward  Zc  (t) . Let  Mf  be 
t f2  fl  1 


the  point  of  Zc(t)  which  coincides  with  B (a)  and  Mr  , 

t T I i 


Mr  the  points 
f2 


defined  by  = -n^  (0+) , M^M^.  “ n^  (0+) . Similarly,  let  Mj.  be  the  point 

which  coincides  with  Bt(o').  When  o'  increases  from  -2  to  a,  M'  generates 

t r 

an  arc  C(o;  r,  t) . Let  C^(a;  z,  t)  and  C^(a;  t,  t)  denote  the  arcs 

described  by  M'  and  M'  , respectively.  Let  S (a;  t,  t)  be  the  open  surface 
f i f2 

generated  by  M'  M'  . It  is  normal  to  Zc(t)  and  its  edge  is  the  closed 
*2  f 1 1 


contour 


M (a)  B f (-2)  bJ  (-2)  M (a)  M (o)  (6.3) 

fl  £t  2t  f2  fl 


The  surface  S(o;  t,  t)  is  intersected  by  all  the  free  vortex  filaments 
which  encounter  C (o ; t,  t).  That  of  these  vortex  filaments  which  passes 
through  has  at  this  point  the  intensity 


V[Ii 


(o\t')  + r2(0,,T,)j 


do'^IVa'  ’ + IV°'  ’ 


The  circulation  of  V in  the  circuit  (6.3)  is  therefore 


a 


t)  + 


r2(av 


:)] 


1^(0,  t)  + r2(a,  t) 


(6.4) 


The  velocity  V generated  by  all  the  bound  and  free  vortex  filaments  is 
irrotational  in  the  domain  [D^  - Z^(t)],  Dg  being  the  exterior  of  E at  t: 

V(M,  t)  = V4>(M,  t)  inside  D - Zr(t)  (6.5) 

e t 
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To  any  point  M^(o)  on  T^(t),  there  corresponds  one,  and  only  one,  point 
B(a)  defined  by 


B(o)  Mf 


T 


VR(P(o,  t),  t')  dt’ 


(6.6) 


Here  P is  the  fluid  point  located  at  Mf(o)  at  t and  t is  the  time  when  this 
fluid  point  was  on  B at  a certain  point  B^.(a)  inside  the  free  vortex  sheet. 
We  have 

4>(M  (a),  t)  - *(M  (a),  t)  = r (a,  t)  + T (o , t)  (6.7) 
f2  *1  1 2 


The  vortex  distribution  0 which  generates  $ may  be  written  in  the  form: 


(6.8) 


(D  + VI)  is  identical  with  the  vortex  distribution  which  is  kinematically 
* ^ -+■ 
equivalent  to  the  given  body  when  no  free  vortices  are  shed.  Let  and 

be  the  velocities  induced  by  V ^ and  V respectively.  We  have 
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<4- 


V.  0 inside  D , curl  V = 2fl  (t)  within  D 
l e 1 E 1 


(VE  - ^ t = V*'(M,  t)  within  D± 


O'*  _ . r \ ^ 

— = (n  • V ) on  Z 

3n  E M,  t e 


2 = V<t^,  with  V4>^(M)  = <V£  - vp  witnin  D. 


= » c)  on  ^ at  t 


♦ ' = Z u.(t)  4>'  (M)  everywhere  within  D. 

j = l J J 


V'2  = (Z,  u'n) 


The  velocity  induced  by  V ^ is  within 


V'  = V - V 
2 VE  1 


The  total  velocity  within  is  thus: 


V - V,  + (V^,  - V,)  = V_  within  D. 

1 E 1 E 1. 


and,  since  its  normal  component  is  continuous  through  Z,  the  boundary 
condition  on  Z 


(v  • ”>« , 1 5 VM-  '> 

e 


is  automatically  satisfied. 
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The  distribution  V " is  equivalent  to  a normal  dipole  distribution  on 
I + If(t).  This  normal  dipole  distribution  generates  the  velocity  potential 


t>2  + with 


^(M,  t) 


«f(M,  t) 


s \\  u,,(m' 

l 

afl 


. t)  T~  uuT  dZ(M’) 


3V  m 


I 


"£(m;(0)'  ° -S^TJS;  d£f(Mf(0>) 


(6.13) 


The  discontinuity  of 


* - ^ + *2  + *f 


(6.14) 


through  I,,  from  (a)  to  M (a)  is  equal  to  p (M'(o),  t) , 

I il  l2  i i 

Equation  (6.8),  we  have 


According  to 


Uf(Mj(o,  t))  = (fj  + r2)o>  t = Tf(Bf(o),  t) 


(6.15) 


Furthermore,  because  of  (6.11)  and  (6.11a),  we  must  have 


4>2(M,  t)  = - $£ (M,  t)  within 


(6.16) 


Hence  p"  is  the  solution  of  the  Fredholm  equation  of  the  second  kind 


- \ u"(M,  t) 


iff  '>^rsFd!:(M> 


*f(M,  t)  (6.17) 


on  E 
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T 

The  vortex  — on  I is  finally  given  by 

^ + fjj  = - n A [-  $1  + (^2 (Mg)  - V^(Mi)) 

- n A [(V*  + V^)M  ^ t - Ve(M,  t)],  (Mel) 
e 

(6.18) 

is  determined  bv 

[V(M  (a),  t)  - V(M,  (a),  t)]  (6.19) 

f2  fj 

Note  that  hitherto  no  assumption  had  been  made  concerning  the  smallness 
of  the  motion,  except  for  the  fixity  of  8 on  E.  As  in  the  steady  case,  the 
above  equations  do  not  entirely  determine  the  fluid  motion.  That  deter- 
mination also  requires  the  generalized  Kutta  condition  to  be  taken  into 
account. 

The  Generalized  Kutta  Condition 

Whether  the  relative  motion  is  steady  or  not,  the  continuity  of  the 
fluid  motion  requires  that  the  three  vortex  filaments  L , L , and  Le  be 

u\  02  * 

orthogonal  to  8 (see  Chapter  4) . We  have  also  seen  in  Chapter  4 that  if 
the  relative  motion  is  steady,  this  entails  the  equivalent  two  formulas  Equa- 
tions (4.21a)  and  4.21b)  and  leads  to  Equation  (4.22). 

Let  us  rewrite  (4.22)  by  taking  into  consideration  the  fact  that  the 
motion  of  the  body  is  no  longer  uniform. 


J(M,  t)  - - (i!  A V ) - f,  + 

e 

+ (V"(M  ) - (M  ) ) ] = 
Z e l 1 


Tf(Mf;  t) 

On  Z^,  the  vortex  

T (M  , t)  = - n A 
M 
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Let,  as  above,  G be  the  resolvent  kernel  of  Equation  (4.22)  in  Chapter 


4 . We  have : 


G (B(o) , M)  MM,  t)  dl(M) 


(6.20) 


Vo)  • '>  - t)]  + 11  <b<0)-  t}  t)d£(M> 


On  the  left-hand  side,  we  have: 


f(M,  t) 


n 


^(t) 


Bf (o ' )Mf (o' ) = j VR(P(o\  t),  t')  dt' 


(6.21) 


The  right-hand  side  is  of  the  form: 


2 u (t)  F (o) 

j-1  2 2 


(6.22) 


where 


U1  ” C’  U2 


v° 


when  the  motion  of  the  body  reduces  to  a uniform  translation  of  speed  c. 

Since  we  are  dealing  with  a small  motion  about  a translation,  we  may 
neglect  the  variation  of  V£  with  t provided  that  we  assume  t » tQ.  We 


therefore  have: 


t - -<*>  B(a')M7(a')  - - c(t  - t) 

O * 
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Let  x - - X(o) , y - Y(a),  ~ - n(a)be  the  characteristics  of  8.  By 
integrating  first  on  E^,  Equation  (6.20)  becomes: 


do'ri(o' ) 


3 

3VB(o*) 


G M>^5^o 


dE(M) 


= c F1(o) 


6 

+ I u (t)  F (o) 
j-1  J J 


(6.23) 


Here  the  six  quantities  u^,  u^»  u^,  u^  xf,  u^  and  u^  — are  small  in 
comparison  with  c. 

Let  us  put 


K(o,  o' ; t - t)  (6.24) 

the  expression  in  square  brackets  in  Equation  (6.23). 

When  u » u = u » u,  ■ u,  >■  u,  M 0,  that  is,  when  the  motion  of  the 
1 2 3 4 5 6 

body  is  a uniform  translation,  (6.24)  reduces  to 


- ~ | da'n(o’)  | K(a,  a';  t - x)  rf(o\  x)  dx  - c F^a) 
-l 


By  changing  x into  t - x'  and  writing  x for  x',  we  obtain: 


+1  r 00 

I do'n(o')  1 K(o , o’;  x)  rf(o',  t - x)  dx  - c F^o) 
-l  0 


and  the  limiting  function  r,(o')  for  t ■ °°  is  solution  of  the  equation 
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r+*  r 

00 

r i 

j ,(.■)[ 

1 K(o,  o' ; t)  dr 

do'  - c Fj(a) 


(6.25) 


Let  us  suppose  that 


Uj(t)  = 0 except  for  j » m 


(6.26) 


For  the  corresponding  function  we  obtain 


K(a,  o';  t)  6T  (o',  t - r)  dr  “ u (t)  F (a) 
r mm 

m 


(6.27) 


Such  an  equation  can  be  studied  by  means  of  either  the  Laplace  trans- 
form or  the  Fourier  transform.  Let  us  use  the  second  method.  We  put: 


u (t)  = I exp  (i2Ttvt)  u (v)  dv 
m I m 


(o',  t)  * I exp  (i2irvt)  61^  (o',  v)  dv 
m » m 


6T,  and  d are  the  Fourier  transforms  of  6T.  and  u . Let 
f m f m 

m m 


(6.28) 


K(o,  o';  v)  - I exp  (-12ttvt)  K(a,  o';  i)  di 


(6.29) 


be  the  Fourier  transform  of  a function  zero  for  t < 0 and  equal  to  K for 
t > 0. 
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By  substituting  (6.29)  into  (6.28),  we  obtain  the  integral  equation: 


-fe  J " 


(a  ) K(o,  a ; v)  6Tr  (a  , v)  do'  = u (v)  F (a) 
t mm 

m 


(6.30) 


This  is  a Fredholm  equation  of  the  first  kind.  For  physical  reasons,  we 
assume  that  there  exists  one,  and  only  one,  solution.  In  this  solution  6F  , 

a m 

the  variables  a and  v cannot  be  separable  because  the  kernel  K depends  on  a 

and  o'  simultaneously.  If  u were  a periodic  function  of  t,  6T  (o,  t) 

m m 

would  be  periodic  too,  but  both  the  amplitude  and  the  phase  would  depend 
on  o. 

The  general  solution  of  (6.20)  is  thus 


6rf(o,  t) 


i f 

j=i  j 


exp  (12-rrvt)  • (o,  v)  dv 

3 


(6.31) 


with  6F.  corresponding  to  the  component  u,  of  (V„,  3 . 

3 Eh 


Hydrodynamic  Forces  Exerted  on  the  Body  in 
a Small  Motion  of  the  First  Kind 

Recall  that  the  velocity  potential  on  is  defined  by 


>(Me,  t)  » *'(Mi,  t)  + u'(M,  t)  + p"(M,  t) 


(6.32) 


*(M  , 
e 


t)  * | ^ Uj(t)  (M±,  t)  + JJ  G (M,  P)  *J(P1,  t)  dZ(P)  | 


'I 


G (M,  P)  *f(P,  t)  dZ(P) 


(6.33) 


In  (6.33),  the  expression  in  braces  would  exist  alone  if  no  free  vortex 
sheet  were  shed  by  the  body. 
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According  to  the  Euler  equation,  the  hydrodynamic  pressure  exerted  on 
the  body  in  the  moving  reference  system  is  given  by 


- P,(M  , t 
p a e 


) - - (i*) 

\ / m' , t 2 

e 

■ ■(*)»■,  t+> 


f (0E  r2>»,  t - 2 VR(V  C) 


1 


* ,n2  2S 

2 (“e  r t 2 '*1 


(T  + f»  + f’^)2 


(6.34) 


(6.35) 


Now  recall  that  on  the  surface  of  the  body  there  are  three  local 
systems  of  axes: 

(n,  3,  t),  (n,  P2.  t’),  (n,  t”,  t") 


T x' t are  in  the  directions  of  f,  respectively.  On  2,  we 

also  consider  the  curves  C !,»  each  of  which  is  tangent  to  the  corresponding 

vector  or  Their  arcs  are  in  the  directions  of  these  vectors. 

We  have 


T = n A V ft  . ini  ?*  ft.  = in!l  > 

T1  n A Vl’  T2  8ol  2’  2 3o"  2 


(6.36) 


The  hydrodynamic  system  of  forces  due  to  the  added  masses  is  the  system 


[pd  ln  n dE]  « p E Uj(t)  [*j(Mi)  + Uj (M)]  ^ dZ(M)j  (6.37) 


in 


The  system  due  to  the  apparent  added  masses  is 


Sln  = [pd  in 


Sd£l'[- 0 {11 


G(M,  P)  *f(P,  t)  dl(P) 


)|  n dl  (M)J 


(6.38) 
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The  system  of  the  quasi-steady  forces  is 


Sq.s  = " Vt}  " {|  P ^(M’  C>  A VV  t)J  dI(M)} 

(6.39) 

- - f.<t)  - [{  p (?!  + + i"  M> 

In  this  formula  for  determining  the  quasi-steady  system  of  hydro- 
dynamic  forces  exerted  on  the  body,  T^(M,  t)  depends  on  the  history  of  the 
motion  since  the  contribution  to  <t>  from  <t>^  is 

l t 

*2  = § | do' x (o' ) J 6rf(o',  x)  KjCM,  o',  t - x)  dx,  with 


K^M,  o',  t - x) 


jj  G<M-  p>  aSJ.sW’W 


(6.40) 


+1  c / 

- I do'x(o')  6rf(a’,  x)  K(o , a',  t - x)  dx  - E u.(t)  F.(a) 
-4.  Lo  j J J 


In  practice,  what  is  called  the  quasi-steady  system  of  forces  is  that 
derived  from 


+SL 

r 

■ fC  -1 

j da'x(a') 

I K1(M,  o' , t - x)  dx 

*2  = ' ^ j 


l t 

4n  J da'x(o')  £ J K (a,  a',  t - x)  dxj  6Tf(a')  = Z (tQ)Fj 


(6.41) 


where  tQ  in  the  last  equation  represents  the  present  time.  The  quasi-steady 
system  of  forces  is  therefore  those  forces  which  would  be  reached  at  t = 00 
with  uj(t)  = Uj(tg)  in  the  interval  [t^,  + °») . 
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1 ->  o 

Let  — denote  the  bound  vortex  corresponding  to  Equation  (6.41)  and 

let  S°  denote  the  quasi-steady  system  which  would  be  derived  from  the 
q.s. 


distribution 


tion  (l,  J . 


We  have 


AS  = Su  - S 
q.s.  q.s.  q.s. 


1 “{(ff  - ^ [h  + h + l(f2  + f2°)] 


■v,  ■“<»>; 

(6.42) 


This  difference  is  the  "deficiency"  due  to  the  difference  Sr^(a,)-6r^(o' , t) . 

For  the  sake  of  simplicity,  we  very  often  accept  the  questionable 

dVE(0) 

assumption  according  to  which  the  accelerations  — and  - — are  small 

dt  dt 

enough  for  ‘ 


+ AS 

q . s 


(6.43) 


to  be  practically  negligible. 

SMALL  MOTIONS  OF  THE  SECOND  KIND 

Example  of  a Double  Model:  Structure 

of  the  Vortex  Distribution 

The  planes  of  symmetry  of  the  body  are  the  (z,  x)-plane  which  is 

"vertical"  and  the  (x,  y)-plane  which  is  "horizontal."  The  contour  C 

intersected  in  the  hull  by  the  (z,  x)-plane  is  described  in  the  direction 

from  Oz  toward  Ox,  Oz  being  vertical  upward.  It  consists  of  the  stem 

EjEEq,  of  the  stern  post  S^SS^  and  of  two  horizontal  segments  of  straight 

lines  EqSq  and  S^Ej  in  the  planes  z = -H  and  z = +H,  respectively.  Let  X 

be  the  curvilinear  abscissa  on  C;  its  origin  is  at  E.  The  abscissas  of 

EE.E  S , S and  S are  as  follows: 

0 1 U 1 

X(S)  - t A,  X(SJ)  - -X(SQ)  - Xjt  X(E1)  - -X(Eq)  = -XQ 
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The  "mean"  velocity  of  the  origin  0 of  the  moving  system  S of  axes  is  in 

the  direction  from  S to  E.  The  cartesian  abscissas  of  and  and  of 

and  Ej  are  -£.  and  +£,  respectively. 

The  contour  C is  thus  identical  with  that  of  the  double  models  con- 

H H 

sidered  in  Chapters  4 and  5.  The  aspect  ratio  is  - or  -.  Here  L is  the 

X L 

length  SE  and  the  aspect  ratio  is  very  small. 

In  the  small  motions  of  the  double  model,  the  components  of  V (0)  on 

E 

the  moving  axes  are  c + u^ , u^,  u^  (c  being  a constant)  and  the  components 
of  are  u^,  u5>  u^.  The  squares  of  u^  u3  and  u^Jt,  u^,  ufeJl  are 
supposed  to  be  negligibly  small  in  comparison  with  c^.  If  these  small 
motions  were  totally  in  the  horizontal  plane,  u^,  u^,  u^  would  be  zero. 

It  is  assumed  that  no  free  vortex  is  shed  from  the  stem.  All  the  free 
vortex  filaments  start  from  or  arrive  at  points  on  EqSqS^E^.  Their  ends 
Bq  and  are  mirror  images  of  each  other  with  respect  to  the  horizontal 
(x,  y) -plane.  The  point  BQ  is  located  on  the  arc  EqSqS  and  the  point  Bj 
on  the  arc  SS^E^. 

A free  vortex  filament  L ^ starting  from  B^  and  arriving  at  B^  at  t has 

at  t + dt  two  short  "horizontal"  branches  B,M,  and  M-  B„  such  that 

1 fl  fo  0 

“B^Mfj  - /C  + dt  fyPjtt),  t')  dtMgS  - f*  + dt  VR(PQ(t),  t')  dt'. 


Pq  and  P^  being  the  fluid  points  which  were  at  B^  and  B^  at  t.  The  fluid 
points  on  were  at  t to  port,  and  therefore,  the  bound  part 

included  in  the  same  vortex  filament  L as  L ^ was  necessarily  located  on  an 
arc  BjBq  located  on  the  port  side  E^  of  E.  On  the  contrary,  if  at  t, 

started  from  B^  and  arrived  at  B^,  the  arc  of  L ^ would  necessarily 

be  in  the  region  y < 0 and  L would  be  closed  at  t + dt  by  a bound  vorted 
filament  located  on  the  starboard  side  of  E. 

This  explains  why,  anticipating  this  circumstance,  we  have  assumed  in 
Chapter  IV,  that  no  vortex  filament  including  a free  vortex  filament  is 
divided  on  E into  two  parts,  one  of  them  located  on  E^  and  the  other  one 
on  E . 
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If  all  the  u ' s were  the  constant  zero,  no  free  vortex  would  be  shed. 

The  two  arcs  of  bound  vortex  filaments  joining  B^  to  on  E^  and  B^  to  B^ 

on  w°uld  be  symmetric  to  each  other  with  respect  to  the  (z,  x)-plane 

and  their  intensities  d^  r^(A)j>  d^  ^(A)  would  be  equal.  This  is  no  longer 

true  in  the  case  of  the  small  motions  under  consideration. 

Let  {t  } denote  the  sequence  of  times  t at  which  d,  F,(A,  t)  and 
n n A 1 

d^  r^A,  are  e<lual  to  each  other.  Suppose  that  the  y-component  of  the 

incident  velocities  at  B.  (-A)  and  B^  (-A)  during  the  interval  (t  , t .,)  is 

1 U n n+1 

positive.  Then  at  any  t e(t  , t , , ) , the  intensity  d.  r,(A,  t)  of  the  bound 

n n+ I A 1 

vortex  Cj(A,  t)  joining  B^  (A)  to  B^  (-A)  is  greater  than  the  intensity 

d,  T (A,  t)  of  the  bound  vortex  C’(A,  t)  joining  B (-A)  to  B (A)  on  E . 

A L L 1 U 2. 

The  difference  d d [r  (A,  t)  - T (A,  t)]  increases  during  the  first  part 
t A I 2. 

(t  , t')  of  t , t ..)  and  decreases  during  the  second  part  (t',  t , . ) . 
n n n n+1  n n+l 

If  t^(t  , t'),  the  positive  variation  d d,[(rt,  t)  - T„(A,  x)]  of 
n n x A 1 l 

d^[Fj(A,  t)  - T2(A,  t)]  determines  the  shedding  of  a free  vortex  filament 

B (-A)  ML  (A,  x;  x+dt)  M (A,  x;  x+dx)  M (A,  x;  x+dx)  B (A)  closed  by  the 
111  I f 0 ^ 

arc  C|(A,  x + dx)  joining  Bq(A)  to  B^  (-A)  at  x + dx  on  E^.  The  fluid 

points  which  were  at  x at  B^  (-A)  and  B^  (A)  are  at  (A,  x;  x+dx)  and 

M,.  (A,  x;  x+dx)  at  x+dx,  the  short  arcs  B,  (-A)  M.  (A^  x;  x+dx)  and  B.  (A) 
f0  1 M 0 

(A,  x;  x+dx)  are  approximately  horizontal.  We  have 

f o 


B.(-A)  M (A,  x;  x+dx)  = V.(-A,x)  dx  x 

1 fl  1 l , 

' -VE(X- 

B„<»)  <A,  IS  T+dT)  - V, (A,  T ) dx  ’ 


x)  dx 


V^(-A,  x)  and  V^(A,  x)  are  the  incident  velocities  at  B^(-A)  and  Bg(A)  at 

time  x.  At  t > x,  these  two  fluid  points  are  at  Mc  (A,  x;  t) , M,.  (A,  x;  t) , 

f 0 f 1 

respectively,  and  the  arc  Mr  (A,  x;  x+dx)  MC(A,  x,  x+dx)  Mr  (A,  t;  x+dx)  is 

f0  f f0 

at  t an  arc  of  free  vortex  filament  whose  position  and  intensity  are 
(A,  x;  t)  Mf (A,  x;  t)  Mf  (A,  x;  t)  and  dT  dA  [F^A,  x)  - T2(A,  x)], 
respectively . 
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If  xc(t',  t the  free  vortex  filament  shed  during  the  time  interval 

n n+i 

(t,  x+dx)  is  analogous  to  the  preceding  one  but  its  intensity  d^  d^ 

[^(A,  t)  - ^(A,  t)]  along  the  arc  B^-A)  Mf  (A,  x;  x+dx)  Mf(A,  x;  x+dx) 

M,  (A,  x;  x+dx)  B (A)  B (-A)  is  now  negative. 

I o U 1 

Whatever  t may  be  inside  the  interval  (t  , t any  point  on  the 

n n+i 

line  B (-A)  M (A,  t ; t)  may  be  considered  the  position  at  t of  a fluid 
i £ i n 

point  which  was  at  B (-A)  at  a certain  time  xe(t  , t) . Thus  it  can  be 

l n 

denoted  (A,  x;  t) . The  interval  t - x is  roughly  proportional  to  the 

distance  from  B^  (-A).  The  intensity  at  (A,  x;  t)  of  the  free  vortex 

filament  whose  support  is  the  line  B (-A)  Mc,  (A,  t ; t)  is  as 

1 £ 1 n 


T 

J df  dx  [^(A,  t«)  _ r2(A,  t')]  = dx  [t^a,  x)  - r2(A,  x)] 


since  dx[T1(A,  tR)  - T2(A,  tR)]  * 0. 


The  points  M (A,  x;  t)  are  defined 
1 0 


in  a similar  manner  on  the  line  B-  Mc  (A,  t ; t) . Along  the  arc  M 

to  tr 


n 


(A,  x;  t)  M (A,  x;  t)  M (A,  x;  t) , the  intensity  d d [r.(A,  x)  - T0(A,  x)] 
r i q x a i i. 

of  the  free  vortex  filament  is  positive  if  xc(t  , t')  and  negative  if 

n n 

,CK‘  Vi>- 


Let  S.(A,  t J 
I n 


t)  denote  the  surface  generated  by  the  arc  M Mc  M 

fl  * f0 

during  the  time  interval  [t  , t) . When 

n 

of  the  free  vortex  Mc  (A,  t - 0;  t)  MC(A,  t - 0;  t)  M..  (A,  t - 0;  t) 

f 1 f f0 

vanishes  and  the  surface  S,  (A,  t ; 

1 n 


t reaches  time  t , , the  intensity 
n+i 


t is  completely  free  for  t > t , , ; it 
n+1  n+1 


moves  with  the  fluid.  For  t'  e(t  t , „)  , there  exists  a surface  S„ 

n+1  n+2  2 

(A,  t t * ) analogous  to  S,(A,  t ; t)  but  starting  from  the  arc  C'(A,  t') 
n+ l in  l 

which  is  the  support  at  t'  of  the  bound  vortex  filament  joining  B^  to  B^  on 

V 

Figure  13  indicates  the  relative  positions  of  the  body  and  of  the 

surface  S (A;  t , t)  at  t e(t  , t ..);  Figure  14  shows  the  relative  positions 
in  n n+ i 


of  the  body  and  of  the  two  surfaces  S^(A,  t^; 
t'  e(tn+l’  tn+2)> 


tn+l)>  S2^'  Si+l’  t ) at 
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If  we  are  dealing  with  a "pure"  swaying  motion,  all  the  incident  velo- 
cities at  Bj(-A)  and  Bq(X)  have  positive  y-components  during  the  even 

intervals  (t„  , t„  , and  negative  y-components  during  the  odd  intervals 
Ip  2p+l 

(t2p+i,  C2p+2^ ’ * * • Therefore,  if,  for  instance,  t e(t2n+i*  c2n+2^ » a11  the 

surfaces  S^X,  t^;  t^)  and  S,,(X,  t^;  t2p+2)  are  free  for  t2p  < t^J 

only  the  surfaces  S-(X,  t„  , , , t.  are  in  contact  with  the  body  along 
Z zn+l  Zn+Z 

the  arcs  C2(X,  t) • 

In  the  case  of  a "pure"  yawing  motion,  V (0)  is  tangent  to  the  (z,  x)- 

£ 

plane  at  any  t.  At  a given  time  t,  the  y-component  of  the  incident  velo- 
cities has  a certain  sign  for  X > X(0)  and  the  opposite  sign  for  X < X(0). 

These  signs  depend  on  that  of  the  z-component  of  ft  . If  the  latter  is 

£ 

positive,  the  surfaces  S2  corresponding  to  the  interval  X > X(0)  and  the 
surface  corresponding  to  the  interval  X < X(0)  are  in  contact  with  the 
body  along  arcs  C^X,  t)  and  C|(X,  t) , respectively. 

In  both  cases,  the  vortex  wake  is  three-dimensional.  The  small  motions 
of  the  second  kind  are  thus  much  more  intricate  than  the  small  motions  of 
the  first  kind. 

In  such  a motion  that  v„  becomes  a constant  different  from  zero  for 

E 

t > tn,  while  w and  & vanish,  then  the  surfaces  S and  S are  carried  by 
U £ £ 1 2. 

the  flow  at  infinity  downstream  from  the  body,  and  in  the  limiting  steady 

motion,  the  vortex  wake  rtduces  to  the  surface  I,  + E,  + El.  considered 

f 1 f0  f 

in  Chapter  4. 


Calculation  of  the  Velocity  Field 

The  total  velocity  V due  to  the  motion  of  the  body  is  defined  inside 
Dg  and  inside  by  a vortex  distribution  V.  It  is  expressed  by  formulas 
analogous  to  those  of  the  preceding  sections. 

We  may  write 


V - p1  + V’2  + P'' 

We  have  first 

vi  ■ (£-  r)*(v  25e) 
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The  velocity  generated  by  V ^ is  zero  inside  D^.  Within  D.,  we  have 

V - V , = V*' 
fc.  1 

We  know  from  discussions  of  the  effective  determination  of  the  vortex 
distribution  (Chapter  3)  how  to  calculate  and  thus  A second  part  of 

V consists  of  vortices  ?'  distributed  on  I.  The  corresponding  vortex 
distribution 

is  equivalent  to  a normal  doublet  distribution  on  L: 

V'2  * (I,  n'  n) 

and  the  density  p'  is  selected  so  that  the  velocity  potential 

‘2-£  If  <M)  d£(M') 

satisfies  the  boundary  condition  *2  = * on  1^. 

The  velocity  thus  satisfies  the  condition 

V,  + 7*1  ■=  V,  + (V  - V,)  = Vr,  within  D 
1 2 1 E 1 E i 

This  shows  that 

n . (V,  + V')=  n . V„  on  I as  well  as  on  Z. 

1 i E e i 

Consequently,  the  vortex  distribution  V ^ + V'2  is  that  with  which  we  should 
deal  if  no  free  vortex  were  shed  by  the  body. 
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In  the  case  of  motions  of  the  first  kind,  the  remaining  vortex  distri- 
bution P"  consisted  of  vortex  filaments  L,  each  L being  located  partly  on 
the  hull  surface  I and  partly  on  the  free  vortex  sheet  Z^.  It  has  been 
shown  that  P " is  equivaler 
1^.  Symbolically,  we  had 


shown  that  P " is  equivalent  to  a normal  doublet  distribution  on  Z and  on 


P'2  ^ (^,  »'  n)  + (Zf,  uf  nf ) . 


The  velocity  (V  - vp  generated  by  P'^  was  expressed  as: 


7*"  + V*  with 


*Jff 


srt'W 


1 


The  density  u"  was  determined  by  the  condition 


*"  (M^  t)  3= 


• f (Mi,  t)  on  Zi 


and  Uj  was  the  solution  of  an  integral  equation  which  states  that  the  pres- 
sure p on  Z is  continuous  through  the  line  B from  which  the  free  vortex 
filaments  are  shed. 

The  same  principle  applies  in  the  case  of  motions  of  the  second  kind 

but  with  some  modifications  because  the  vortex  wake  is  now  three-dimensional. 

Let  us  suppose,  for  instance,  that  t c (t„  , t„ 

2n  2n+ 1 

The  free  vortex  filaments  shed  at  times  t prior  to  t^  are  located 

within  domains  W„  , W„  , , ; each  domain  W_  contains  all  the  free  vortices 
2p  2p+l  2p 

shed  during  the  time  interval  [t„  , t„  and  each  domain  W contains 

2p  2p+l  _ 2p+l 


all  the  free  vortices  shed  during  the  time  interval  [t„  . . , t„ 

2p+l  2p+2 


The 

) . In  a pure 


domain  W„  is  generated  by  the  surfaces  S (X,  t„  . , , t„ 

2p  2p+l  2p+z 

swaying  motion,  all  the  surfaces  S (X,  2p,  2p+l)  are  of  type  Sj  and  all  the 


105 


surfaces  S (X,  t„  , , t ,.)  are  of  type  S . In  a pure  yawing  motion,  each 

zp-rl  Zp+ii  L 

domain  W„  , W„  is  subdivided  into  two  domains;  for  instance: 

2p  2p+l 


w = w'  + w" 

w2p  W2p  + W2p 

where  is  generated  only  by  surfaces  and  W'^  only  by  surfaces  S^. 

The  domains  VJ  and  W„  , . are  disjoint. 

2p  2p+l  J t 

Furthermore,  the  vortex  wake  contains  a domain  W„  whose  structure 

2n 

differs  from  that  of  domain  or  The  explanation  is  that  a surface 

such  as  S.  (X,  t ; t)  includes  an  arc  of  bound  vortex  filament  whose 
1 2n 

intensity  is  d^[(r^(X,  t)  - ^(X,  t)]  and  whose  support  is  C|  (X,  t) ; in  an 
analogous  manner,  if  W^n  also  contains  surfaces  (X',  t^n;  t) , each  of 
these  surfaces  includes  an  arc  of  vortex  filament  of  intensity  t) 

- T^(Xr,  t)]  whose  support  is  the  arc  C^CX',  t) . 

The  velocity  induced  by  the  vorticity  distributed  inside  a domain  W 
is  given  by 


V (M,  t)  = curl 
2p 


if 


uf(M’,  t) 


dW  (M’) 
2p 


with  cu.  = — Tc  if  M'  belongs  to  the  cover  (E  + Zc  + E')  of  W . 
£ £ t £}£()£  Zp  ^p 

Similarly,  we  have 


V2p+1(M*  0 = CUrl 


uf(M',  t) 


Let  us  now  consider  the  velocity  V induced  by  the  vortex  filaments 
t t 

located  inside  W,^.  In  a pure  swaying  motion,  is  adjacent  to  part  Ej^ 

of  the  portside  E of  the  hull  downstream  from  the  arc  C!  (X  , t) . In  a 

1 t t 

pure  yawing  motion,  is  subdivided  into  two  parts  and  W"  respectively 


r J o — I i-utu  tww  ivc»i.ua  n auu  n „ l C3UCL 

2n  2n  2n 

adjacent  to  a part  E|  of  E^  and  to  a part  E^  of  E^.  The  arcs  of  vortex 
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filaments  of  intensity  t)  - r2 ( A , t)]  located  on  Z|  and  those  of 

the  vortex  filament  of  intensity  d^fr^CA,  t)  - Fj(A,  t)]  located  on  form 
two  vortex  sheets.  Consequently  is  discontinuous  through  Zj  and  through 
zp  Its  jump  is  tangent  to  Z'^  or  to  Zj,  an<^  zero  through  Z - (Z^  + Zp . 

The  velocity  is  given  by 


2n 


= curl 


tj((, 

2n 


o)f(M' 


t) 


MM' 


dW^  (M') 
Zn 


with  a),  = — T.  on  the  cover  (Z^  + Z* 
f e f f i f o 


+ r‘  + EJ  + £•)  of  »‘n. 


The  resulting  velocity 


V (M,  t)  = Z [v2r(M,  t)  + V2r+1(M,  t)]  + V^M) 
r<n 


is  irrotational  within  D^.  We  may  write: 


V (M,  t)  = V$  (M,  t)  when  M e D. 
w w 1 

Obviously  V (M  , t)  is  not  tangent  to  Z . Consequently,  the  vortex 
we  e 

distribution  V 2 contains  closed  bound  vortex  filaments  located  on  Z.  They 
induce  inside  a velocity  V’^  which  is  irrotational  everywhere  except  on 


Z and  which  meets  the  requirement 


V'2  (M,  t) 


- V (M,  t)  at  every  point  M e D 
w i 


This  complementary  vortex  distribution  is  equivalent  to  a normal 
let  distribution  on  Z.  Let  4>2  c 
doublet  distribution.  We  may  write 


doublet  distribution  on  Z.  Let  4>”  denote  the  potential  generated  by  this 


V'2  = V4>2>  with 


tJ  Jjf 
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and  p"  is  determined  by  the  regular  Fredholm  equation  which  states  that 


= - $ at  every  point  MJ  on  I, 
2 w J i i 


The  corresponding  bound  vortex  — T'^  is 


i.  £n 
e 


T'2 (M,  t)  - - ^ A[v$^(Me,  t)  - t)] 


The  total  vortex  distribution  P”  may  symbolically  be  written  in  the  form: 


(*•  1 *0 


+ I C(w2r*  ^ + (W2r+l* 
r<n 


«f)] 


(W2n’  “f> 


The  total  bound  vortex  on  Z is 


l 

T' 

on 

Z! 

c 

*1 

1 

1 

T' 

on 

Zl 

e 

f2 

2 

where  and  belong  to  and  respectively.  T is  tangent  to  Cj 
(A,  t)  on  Z^  and  to  C^A,  t)  on  Z^«  and  the  total  intensity  of  the  arc  of 
bound  vortex  filament  which  is  the  support  of  ? is  d^T^CA,  t)  on  and 

^A^2^’  ^ 0n  ^2 " 

The  Generalized  Kutta  Condition  in  Small 
Motions  of  the  Second  Kind 

Let  Bq(A)  be  a point  on  E^S^S,  and  let  P,  Q denote  two  points  very 
close  to  Bq  on  C * ( A , t)  and  C^CA,  t)  respectively,  (see  Figure  15).  We 
assume  that 

d^a,  t)  > dxT2(A,  t)  (6.44) 
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and,  therefore,  that  is  the  end  of  a free  vortex  filament  L of  intensity 

0 f o 

d.[r!(A,  t)  - T'(A,  t)].  Let  P , Q be  two  points  on  I with  P~ ? = n_(CH-), 
Ai  z e e eer 

Q?  = n^  (0+) . The  generalized  Kutta  condition  means  that 


lim  £ [p  (P  , t)  - p (Q  , t)]  = 0 
p d e d e 


(6.45) 


when  P,  Q tend  to  (Figure  15).  To  express  this  condition  more  explicitly, 

consider  a point  b very  close  to  P , b P being  in  the  direction  of  . 

t e e f0 
This  point  thus  belongs  to  W ..  . . , . . , 

. . We  also  consider  a point  P on  the  lower 
zn 

side  of  Zc  , with  b P'  orthogonal  to  I.  . The  flow  is  irrotational  at  P' 

*0  f0 
and  at  Q^.  The  condition  (6.45)  is  thus  equivalent  to 


lim—  [p  (b,  t)  - p,(P',  t)]  = 0 when  P,  Q tend  to  B„ 
p d d 0 


(6.46) 


W’an  A decreases  so  that  B moves  in  the  positive  x-direction,  d 

U A 

Lr! (A,  t)  - T'(A,  t)]  decreases  and  vanishes  at  a certain  point  B'(A'). 
i 2 ^ 0 

Simultaneously,  P*  describes  an  arc  P'Pq  and  b an  arc  b b^.  Consider  the 
path  P consisting  of  the  arcs  P'Pq,  PQbQ>  b^b.  According  to  the  Euler 
equation,  we  have: 


| [Pd(b,  t)  - pd(P',  t)]  - - ■§£  J fas  - J (u  A VR)  ds  - | [VR(b)  - vJ(P')] 

P P 

(6.47) 


At  b,  u is  tangent  to  the  arc  M , (A,  t - 0;  t)  M,  (A,  t - 0;  t)  and  thus  is 

tl  + f0 

nearly  parallel  to  Cj(A,  t).  Hence  u A VR  is  normal  to  Z,  and  when  P tends 
t°  Bq,  it  is  normal  to  the  arc  b^b.  As  u = 0 on  P'Pq,  and  P^b^  is 
infinitely  small,  we  have 


J (2  A VR)  . 


ds  = 0 
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Un  the  other  hand,  I V 3s 
Jp+bP' 


>.  f 

JP 


becomes  equal  to  the  flux  of  T through  the 

t n 


closed  contour  P'P^b'bP',  and  we  have 


A 

j v ds  =J  djr^x,  t)  - r2(x,  t)]  = (x,  t)  - r2(A,  t) 


Ao 


Finally,  since.  V and  V_  are  the  relative  velocities  at  B.  on  and  on 
K j K.2  0 1 

£2>  respectively,  Equation  (6.46)  reduces  to 


[^(A,  t)  - r2(A,  t)]  - j [v*  (Bq,  t)  - V*  (B0,  t)]  = 0 (6.48) 


This  is  the  generalized  Kutta  condition.  It  should  be  compared  with 
the  result  for  the  general  Kutta  condition  given  in  Chapter  4 for  bodies  in 
an  oblique,  uniform  translation  in  the  horizontal  plane. 
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CHAPTER  7:  REMARKS  ON  SOME  APPLICATIONS  TO  SHIPS 


INTERACTION  OF  APPENDAGES  WITH  THE  HULL 

These  lecture  notes  would  be  very  incomplete  if  nothing  were  said  about 
the  appendages  witn  which  a ship  is  fitted,  and,  more  particularly,  those 
whose  role  is  to  allow  changes  of  heading  or  depth  or  to  give  to  the  ship 
good  coursekeeping  characteristics. 

In  fact,  the  first  systematic  researches  on  maneuverability  began  a 
century  or  so  ago  with  the  testing  of  isolated  rudders.  At  that  time, 
nothing  was  known  about  lifting  surface  theory  and  the  experiments  were 
therefore  very  useful.  In  France,  experiments  were  carried  out  by  Joessel 
at  Indret  on  the  River  Loire.  His  empirical  formulas  were  used  for  several 
decades  by  naval  architects  as  guides  for  predicting  the  horizontal  components 
of  the  hydrodynamic  forces  exerted  on  the  rudders  and  for  determining  the 
size  of  the  axle,  the  position  and  dimensions  of  the  bearings,  and  the 
structural  design  of  the  stern  post. 

Experiments  on  isolated  rudders  are  now  obsolete  inasmuch  as  the 
properties  of  a given  rudder  depend  considerably  on  the  interaction  between 
rudder  and  hull.  For  example,  the  flow  about  a rudder  located  behind  a skeg 
is  very  different  from  the  flow  about  the  same  rudder  when  there  is  a large 
gap  between  the  hull  and  the  leading  edge  of  the  rudder  (see  Figures  16a 
and  16b).  The  effect  of  the  rudder  is  delayed  in  the  first  case,  but  is 
quasi-immediate  in  the  second.  In  the  first  case,  the  rudder  and  the  skeg 
form  a whole;  in  the  second  case,  the  circulation  of  the  velocity  in  a 
circuit  surrounding  the  rudder  is  established  almost  instantaneously.  The 
time  to  establish  the  circulation  around  an  obstacle  is  the  time  it  takes 
the  flow  to  cover  the  length  of  the  obstacle. 

It  is  known  that  the  moment  about  the  ship  CG  of  the  hydrodynamic  force 
on  the  hull  plus  rudder  is  generally  greater  than  the  moment  due  to  the  con- 
tribution from  the  rudder  alone.  This  effect  is  greater  when  the  vertical 
gap  between  the  hull  and  the  upper  edge  of  the  rudder  is  small  and  the 
transverse  sections  of  the  hull  in  the  region  of  the  rudder  are  flat.  The 
vicinity  of  the  hull  prevents  the  upper  edge  tip  vortex  from  shedding  and 
increases  the  apparent  aspect  ratio  of  the  rudder.  It  may  even  happen  that 
the  bound  vortices  on  the  rudder  are  continued  on  the  hull  until  the  free 
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surface.  In  this  case,  the  amplification  factor  may  be  greater  than  2 (Figure 
17). 

A rudder  may  give  rise  to  separation.  The  smaller  the  aspect  ratio  is, 
the  greater  the  angle  at  which  separation  occurs.  Therefore,  the  best 
rudder  is  not  necessarily  one  with  the  highest  aspect  ratio.  Furthermore, 
course  stability  often  requires  the  addition  of  a fin  ahead  of  the  rudder. 

If  the  rudder  is  joined  to  a fixed  strut  whose  length  is  equal  to  one- 
quarter  or  one-third  of  the  chord  of  the  rudder,  the  resulting  rudder  works 
somewhat  like  a cambered  wing  at  rudder  angles  other  than  zero.  This 
decreases  the  risk  of  separation  and  increases  the  lift  coefficient  of  the 
rudder  (Figure  18). 

On  modern  submarines,  the  rudder  consists  of  two  parts  that  are 
approximately  symmetric  to  each  other  with  respect  to  the  longitudinal  axis 
of  the  vessel.  Clearly  the  bound  vortices  on  the  two  parts  are  connected  by 
bound  vortices  on  the  afterbody  (Figure  19).  In  this  region,  the  two  parts 
of  the  rudder  must  be  regarded  as  included  in  the  longitudinal  contour  of 
the  ship.  The  same  is  true  of  the  sail  of  a submarine. 

In  consideration  of  the  preceding  chapters,  I do  not  believe  that  it  is 
possible  to  set  forth  a general  rule  which  can  take  into  account  the  effect 
of  the  apoendages  on  the  maneuvering  qualities  of  a given  ship.  Each  case 
should  be  considered  separately  before  deciding  whether  it  has  to  be  incor- 
porated into  the  hull  or  be  regarded  as  separated  from  it. 

UNSTEADY  MOTIONS  ABOUT  A WING  PROFILE 

It  has  been  shown  in  Chapter  4 that  by  virtue  of  the  generalized 
Kutta  condition,  the  relative  velocity  on  both  sides  of  a wing  with  a finite 
aspect  ratio  and  a finite  thickness  is  tangent  to  the  trailing  edge  at  any 
point  from  which  a free  vortex  is  shed.  From  this  it  follows  that  the 
assumption  of  an  almost  inviscid  fluid  necessarily  contradicts  either  the 
continuity  of  the  motion  or  the  Lagrange  theoren.  The  Lagrange  theorem 
itself  is  a consequence  of  the  Helmholtz  theorem  which  implies  continuity 
of  the  motion. 

Since  the  Kutta  condition  requires  that  the  relative  velocity  be  zero 
at  the  trailing  edge,  there  is  an  analogous  contradiction  in  the  case  of  a 


112 


wing  profile  when  the  suction  side  and  the  pressure  side  are  not  tangent 
to  each  other  at  the  trailing  edge.  We  will  now  show  that  the  Kutta  con- 
dition  entails  another  unexpected  contradiction  when  the  relative  motion 
about  the  wing  profile  is  unsteady. 

Let  P denote  the  wing  profile,  M a point  describing  P,  n^  the  unit 
vector  normal  to  P at  M in  the  inward  direction,  and  Pg  the  profile 
described  by  the  point  defined  by  MM^  = - 11^(0+).  Let  B,  be  the 
trailing  edges  of  P and  P , respectively  (Figure  20). 

Let  (Ox,  Oy)  be  a system  S of  axes  moving  with  P,  and  let  (O'x,  O'y) 
be  a system  S'  of  fixed  axes.  We  select  S'  so  that  it  coincides  with  S at 
the  time  t under  consideration.  Let  M'  be  the  point  of  S'  with  which  a 
point  M of  S coincides  at  t. 

We  suppose  that  the  velocity  V of  0 has  on  the  X'-axis  a component  in 

h 

the  negative  x'-direction  and  that  the  fluid  is  at  rest  at  infinity  ahead 
of  profile  P.  Let  V(M',  t)  = V$(M' , t)  be  the  velocity  of  the  fluid  point 
located  at  M at  time  t. 

We  have 

\ 

V(M',  t)  = Ve(M,  t)  + Vr(M,  t)  with  $e(M,  t)  = V£(0,  t)  + n£(t)  A 0M 

where  V is  the  relative  velocity  of  the  fluid. 

K 

Let  T(t)  be  the  circulation  of  V in  a closed  circuit  surrounding  the 

dr  dr 

profile.  If  — ^ 0,  a free  vortex  of  intensity  — — dt  is  shed  from  B 
at  dt  e 

during  the  time  interval  [t,  t + dt).  This  free  vortex  is  carried  by  the 

fluid,  and  a fluid  point  P(B  (x))  which  was  at  B at  x < t is  at  t in  a 

e e 

point  such  that 

t 

Vf.  ' J VP(V  ’>•  '’>  dt’ 

T 

V 

Let  t denote  the  time  at  the  beginning  of  the  motion.  The  line  L described 
^ t 

by  the  points  since  time  t^  makes  a "cut"  in  the  plane  (Ox,  Oy)  in  the 
sense  of  the  theory  of  functions  of  a complex  variable  (x  + iy) . 
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Let  L ^ be  the  lower  and  upper  edges  of  the  cut,  respectively.  Let 

U be  a geometric  point  on  L,  and  let  v be  the  unit  vector  normal  to  L at 

MO 

Pq  in  the  direction  from  L ^ toward  L^.  Let  p^,  P2  denote  the  intersections 

of  v with  1 and  L , respectively. 

M 0 1 ^ £ 

Let  £ denote  a path  from  p^  to  p^  following  first  from  p^  to 

and  then  L„  from  M*i  to  p_;  tn  is  the  time  of  the  beginning  of  the  motion, 
r tg  / U 

Let  Pq  be  the  fluid  point  located  at  p^  at  t.  This  fluid  point  was  at  B^  at 
a certain  time  t and 


j V 2s  = J - ^ dt  = - r(x) 


(7.1) 


Let  4>  denote  the  velocity  potential  in  the  domain  exterior  to  the  closed 

contour  B B M B . The  points  B and  B are  infinitely  close  to  B 

ej  ftQ  e2  e ej  ej  e2  e 

on  the  pressure  and  suction  sides,  respectively,  and  the  arc  B M B is 
. e2  e ej 

the  profile  itself.  The  velocity  potential  <t  is  uniform  in  that  domain 
since 


J 

B M B 


ei  e e 


V 2s  = r(t) 
2 


Mft0B2 


\/  2s 


but  it  is  discontinuous  across  L since 


(7.2) 


4>(p^,  t)  - 4>(p',  t)  = J V 2s  = - T (t) 

i 

Let  Pq,  Pj,  P2  be  the  fluid  points  respectively  located  at  t at  p^, 
p^,  p^.  The  trajectory  of  P^  from  t to  t is  the  arc  B^p^  of  1,  and  the 
velocities  of  the  geometric  points  p^  and  p^  are  equal  to  the  velocity  of 
Pq.  Let  V(pj,  t) , V(p2,  t)  be  the  velocities  of  P^  and  P^  at  t.  We  have 
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This  shows  that  at  every  t,  the  pressure  is  continuous  through  L.  This 
result  is  in  agreement  with  the  general  theory  since  the  vortices  located 
on  L are  free. 

Let  us  consider  now  two  points  M and  M respectively  located  on  the 

e!  e2 

pressure  and  suction  sides  of  P . For  any  point  Me»  we  have 


i >d‘V  c>  ‘ 


H <M\  t)  - | V2(M  , t)  + C(t) 
dt  e l.  e 


- ~ (M^,  t)  + V£(M,  t)  . V(Me,  t)  - |[Ve(M,  t)  + VR(Me,  t)]2+C(t) 

- ^ (M;,  t)  + j VR (M,  t ) - j VR(Me,  t)  + C(t) 


Consequently, 


i *>  - »d<V  t>] 


--[3?  c>] 


(7.8) 


* 7 '>  - VE(V  t>]  - 7 [VV  C)  - VV  t)] 


Thus  when  M -*•  B and  M -+•  B , we  obtain 
e i e i e2  e2 


7 ^d'V  - Pd <B.2 • ‘>1  ‘ - 57  - 7 tvR(Bei-  ‘>  * VR(V  t)]  (7'9) 


But  because  of  the  Kutta  condition,  we  have 


IVV  C)l  “ IVBe2‘  C) 


(7.10) 


The  common  value  of  these  two  relative  velocities  is  zero  if  both  sides  of 

P are  not  tangent  to  each  other  at  B . 
e e 

Consequently 


p [pd(V  ° ■ pd(V  t)] " 


(7.11) 
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It  can  easily  be  shown  that  on  and  ^ on  L0  are  infinitely 


e2 


close  to  B and  to  B , respectively, 

e 1 e2 


7 [pd(V  r)  - pd(uei-  t)]  ' - 2 ST 
7 tpd<V  c) 


(7.12) 


p,(u  , t)] 
d 62 


+ - — 
2 dt 


I am  unaware  as  to  whether  the  rather  surprising  results  expressed  by 
Equations  (7.11)  and  (7.12)  have  previously  been  mentioned  by  other  investi- 
gators. The  jump  of  the  pressure  from  the  pressure  side  to  the  suction  side 
at  the  trailing  edge  of  a wing  profile  in  unsteady  relative  motion  is  an 
effect  of  the  apparent  added  masses  due  to  the  shedding  of  a free  vortex 
sheet.  At  the  instant  t^  of  the  beginning  of  a finite  perturbation,  this 


jump  is  infinite  for 


It  probably  decreases  with 


the  aspect  ratio.  We  have  shown  [8]  that  in  the  case  of  wing  profiles,  the 
apparent  added  masses  considerably  affect  the  pressure  distribution  on  the 
profiles.  We  may  wonder  whether  it  is  not  a source  of  vibration  and  noise 
in  the  case  of  wings  or  of  propeller  blades  with  large  aspect  ratio. 


SNAPROLL 

When  a submarine  is  in  a turning  motion  about  a vertical  axis,  the 
system  of  its  inertial  forces  and  of  the  hydrodynamic  forces  exerted  on  it 
contains  a couple  whose  axis  is  longitudinal.  This  couple  is  due  to  the 
centrifugal  force  which  is  horizontal  and  passes  through  the  center  of 
gravity  G and  to  a centripetal  force  which  is  equal  and  opposite  in  direction 
to  the  former  and  passes  through  the  center  of  volume  C of  the  submarine. 

As  G is  located  below  C,  this  couple  creates  a heeling  moment  to  the  interior 
side.  Since  the  motion  of  gyration  is  delayed  with  respect  to  the  time  at 
which  the  rudder  has  reached  the  ordered  angle,  one  would  expect  the  heel 
to  increase  monotonically  until  the  existing  couple  is  cancelled  by  the 
couple  due  to  the  weight  and  to  the  hydrostatic  force. 

However,  this  is  not  what  happens  with  modern  high-speed  submarines. 

In  general,  these  ships  are  fitted  with  two  rudders  approximately  symmetric 
with  each  other  in  terms  of  a horizontal  plane  and  with  two  diving  planes 
symmetric  in  terms  of  the  longitudinal  plane.  The  vessel  heels  (to  port  if 
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che  rudders  go  to  port)  almost  immediately  after  the  rudders  reach  the  de- 
sired angle,  and  a strong  overshoot  occurs.  The  maximum  of  the  roll  angle 
may  be  twice  or  three  times  the  final  heel  angle,  which  is  reached  after 
some  roll  oscillations.  Figure  21  roughly  sketches  the  phenomenon. 

Of  course,  when  the  rudders  are  at  the  angle  a,  the  incident  flow  on 
the  diving  planes  is  no  longer  symmetrical  with  respect  to  the  longitudinal 
plane  of  symmetry  of  the  vessel.  Thus  we  can  ask  whether  the  rudders  may 
generate  a sheltering  effect  on  the  diving  planes,  consisting  of  a lift  in 
the  upward  direction  on  the  starboard  diving  plane  and  in  the  downward 
direction  on  the  port  diving  plane.  A simple  experiment  on  an  isolated  set 
of  two  rudders  and  two  diving  planes  could  probably  determine  whether  the 
above  explanation  is  correct.  Personally,  I feel  that  another  explanation 

is  also  possible.  The  horizontal  lift  of  the  rudders  is  generated  quasi- 

2 

instantaneously  and  it  varies  like  V . In  the  first  phase  of  the  motion, 

h 

the  vessel  experiences  an  acceleration  to  starboard  (if  the  rudders  are  to 

port)  and  after  a time  t very  close  to  the  beginning  of  the  motion,  the 

sail  has  an  angle  of  attack  proportional  to  V x (Figure  22).  The  quasi- 

3 E 

steady  lift  on  the  sail  is  thus  0(V  x).  The  real  lift  on  the  sail  is 

L 

greater  because  of  the  effect  of  the  apparent  added  masses  due  to  the  free 
vortices  shed  by  the  sail  as  mentioned  in  the  previous  section.  Possibly 
this  may  generate  a heeling  moment  of  importance  and  thus  causes  the 
overshoot  shown  by  experiments. 

"PITCH  UP"  OF  SUBMARINES 

When  a submarine  is  turning  in  a horizontal  plane,  the  system  of  hydro- 
dynamic  forces  exerted  on  it  contains  a moment  about  a transverse  axis.  The 
effect  of  this  moment  seems  to  be  an  increase  in  the  apparent  weight  of  the 
stern  or  a decrease  in  that  of  the  bow. 

\ 

In  the  past,  the  asymmetry  of  the  submarine  with  respect  to  the 

horizontal  plane  of  the  center  of  volume  was  important.  There  was  one 

rudder  only,  and  the  pitch-up  condition  could  be  explained  by  a suction  on 

the  flat  part  of  the  hull  above  the  rudder  due  to  the  circulation  around  the 
9 

rudder. 

__ 

Isabelle,  M. , "Remarques  sur  une  cause  perturbatrice  de  l'4quilibre 
dynamique  d'un  sous-raarin  en  plongde,"  Bull.  Ass.  Tech.  Mar.  et  A4ron.,  Paris 
(1938). 
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Today's  high-speed  submarines  are  much  more  symmetrical  and  the  inter- 
actions with  the  hull  of  each  of  the  two  rudders,  which  are  in  opposite 
directions,  could  cancel  each  other.  However,  a pitch-up  moment  still 
exists.  In  fact,  the  heel  caused  by  the  gyration  modifies  the  structure  of 
the  bound  vortex  half-rings  on  the  forebody  since  the  bow  is  inside  the 
circle  of  gyration.  The  suction  side  on  the  forebody  is  starboard  if  the 
gyration  is  to  port  and  because  of  the  heel  a lift  in  the  upward  direction 
may  appear  on  the  forebody. 

COURSE  STABILITY 

The  course  stability  of  the  motions  in  a horizontal  plane  and  also 
in  the  vertical  plane  are  very  important  properties  of  a ship.  The 
phenomenon  of  instability  of  certain  surface  ships  has  been  known  for  many 
years,  but  it  would  be  impossible  to  understand  the  phenomenon  without  the 
help  of  an  approximate  theory.  This  theory  has  suggested  means  for  experi- 
mentally determining  the  forces  acting  on  the  hull,  and  now  the  condition 
to  be  satisfied  is  rather  well  known.  Although  the  control  systems  now 
available  permit  the  ship  to  keep  a route  closer  to  the  required  one,  it  is 
still  desirable  that  ships  be  naturally  stable.  Stability  decreases  the 
power  necessary  to  sustain  speed.  In  the  case  of  submarines,  the  instability 
in  the  vertical  plane  is  dangerous  and  the  instability  in  the  horizontal 
plane  induces  sway  and  yaw  which  generate  rolling  motions. 

The  forces  exerted  on  a ship  model  are  generally  measured  in  conditions 
which  do  not  permit  the  history  of  the  motion  to  be  taken  into  consideration. 
Thus  stability  is  studied  by  assuming  that  the  forces  are  quasi-steady.  The 
equations  are  differential  rather  than  integral-differential  as  they  should 
be. 

In  his  thesis,  which  concerns  surface  ships  only,  Casal  attempted  to 
take  into  account  the  terms  which  depend  on  the  history  of  the  motion.  He 
concluded  that  the  effects  of  these  terms  are  rather  small;  in  general,  the 
stability  - or  the  instability  - would  be  reinforced  if,  according  to  the 
differential  equations,  the  ship  is  stable  or  unstable.  Nevertheless,  he 
showed  that  in  certain  circumstances,  the  integral-differential  terms  may 
decrease  the  damping  coefficients  so  that  the  ship  has  self-yawing  and 
swaying  motions. 
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The  Casal  conclusions  are  based  on  the  concept  that  during  small  motions, 
all  the  free  vortex  filaments  are  gathered  in  a unique  free  vortex  filament 
shed  from  the  stern  post.  For  reasons  explained  in  Chapter  6,  it  seems  to 
me  that  this  scheme  is  too  simple  and  that  the  effect  of  the  integral- 
differential  terms  needs  new  refinements. 

EXPERIMENTAL  STUDIES  ON  MANEUVERABILITY 

All  that  we  know  on  the  subject  of  maneuverability  comes  from  experience 
acquired  either  by  observing  ship  motions  in  a seaway  or  by  carrying  out 
tests  on  ship  models. 

At  the  beginning,  the  main  problem  was  to  determine  the  characteristics 

of  the  steady  gyration  of  surface  ships.  The  hydrodynamic  forces  exerted  on 

the  models  were  not  measured.  Then  course  stability  was  considered.  Recall 

2 

the  mention  in  Chapter  1 that  criteria  were  formulated  by  Contensou  and  by 
3 

Davidson  and  Schiff  . Unfortunately  these  criteria  were  difficult  to  apply 
because  they  were  expressed  in  terms  of  the  forces  exerted  on  the  ship  and 
such  forces  were  unknown.  A finding  by  Dieudonne  is  very  important,  namely, 
that  a ship  is  stable  only  if  the  resulting  force  in  a motion  of  translation 
at  a very  small  angle  of  attack  pierces  the  longitudinal  plane  of  symmetry 
at  a point  located  downstream  from  the  intersection  with  that  plane  of  the 
resultant  of  the  centrifugal  force  and  of  the  hydrodynamic  transversal  force 
in  a gyration  at  zero  drift  angle.  Dieudonne  also  showed  that  this  condi- 
tion means  that  the  ratio  wL/U  in  steady  gyrations  is  a uniform  function  of 
the  rudder  angle  6.  It  has  been  mentioned  in  the  previous  section  that  this 
condition  does  not  take  into  account  the  history  of  the  motion  and  that  in 
some  "pathological"  cases,  self-yawing  motions  may  happen  at  zero  rudder 
angle.  Very  simple  tests  in  a rectilinear  tank  or  (and)  in  a steering  tank 
will  suffice  to  determine  whether  a ship  model  is  stable  or  not. 

The  planar  motion  mechanism  originated  at  DTNSRDC  more  than  15  years  ago 
was  designed  to  measure  the  hydrodynamic  forces  exerted  on  submarine  models 
either  in  oblique  translation  in  the  horizontal  plane  and  in  the  vertical 
one  or  in  forced  periodic  swaying,  yawing,  heaving,  and  pitching  motions. 

The  torque  about  the  longitudinal  axis  of  the  model  can  also  be  measured 
simultaneously.  The  planar  motion  mechanisms  permit  determining  the  response 
of  a free-running  model  to  slight  variations  of  the  rudder  or  of  the  diving 
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planes  under  the  assumption  that  the  hydrodynamic  system  of  forces  reduces 
to  the  quasi-steady  one.  Brard  has  shown^  that  this  technique  can  also 
provide  informations  on  the  effect  of  the  history  of  small  motions  about  a 
uniform  motion  of  translation  provided  that  the  range  of  the  w values  is 
sufficiently  large.  If  the  reduced  angular  velocity  wL/U  acts  on  the 
dimensionless  system  of  forces,  the  part  due  to  the  added  masses  has  to  be 
derived  from  the  measurements  at  high  wL/U  values  and  the  out-of-phase 
forces  from  the  measurements  at  small  values  of  this  parameter. 

Planar  motion  mechanisms  for  surface  ships  are  now  available  at 
several  establishments;  generally  the  amplitudes  of  the  horizontal  periodic 
motions  are  sensibly  greater  than  those  of  the  motions  permitted  by  the 
initial  variants  of  this  device. 

Tests  on  models  in  forced  motions  do  not  fill  all  the  needs  for 
several  reasons.  In  the  first  place,  the  present  state  of  the  theory  does 
not  permit  predicting  all  the  phenomena  to  which  the  maneuvering  ship  can 
give  rise.  This  is  clearly  the  case  when  the  problem  is  to  determine  the 
behavior  of  the  ship  in  waves,  in  restricted  waters,  or  during  fierce 
breaking.  Perhaps  this  is  still  more  evident  when  six  degrees  of  freedom 
have  to  be  considered  simultaneously,  as,  for  example,  in  the  case  of  sub- 
marines or  of  nonconventional  ships  such  as  hydrofoil  craft  or  surface 
effect  vehicles.  A review^  is  available  of  the  various  methods  used  in 
the  Paris  Naval  Tank  together  with  a description  of  the  free  submarine  models 
that  we  are  testing  on  the  bay  of  St.  Tropez  on  the  French  Mediterranean 
coast  for  experimental  determination  of  the  maneuvering  qualities  of  sub- 
marines. 

The  purpose  of  such  tests  of  free-running  submarine  models  is  not  only 
to  compare  the  real  motions  to  those  computed  by  using  analog  or  (and)  hybrid 
computers;  we  also  wish  to  discover  all  the  unexpected  phenomena  which  may 
occur,  for  example,  in  maneuvers  to  be  attempted  in  cases  of  emergency. 


^Brard,  R. , "A  Vortex  Theory  for  the  Manoeuvering  Ship",  Proc. 

Fifth  Symposium  on  Naval  Hydrodynamics,  Bergen  (1964). 

Brard,  R.,  et  al.,  "Le  module  libre  de  sous-marin  du  Bassin  d’Essais  des 
Carfenes,"  Bull.  Ass.  Tech.  Mar.  et  A6ron.,  Paris  (1968). 
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FINAL  REMARKS  AND  CONCLUSION 


During  the  last  two  or  three  decades  considerable  progresses  have  been 
made  in  the  field  of  ship  maneuverability.  They  have  resulted,  essentially, 
from  the  pressure  of  needs.  Today  the  final  purpose  is  very  ambitious  since 
it  is  to  control  any  maneuver  in  a complete  manner.  Desired  changes  of  course 
(or  of  depth  in  the  case  of  submarines)  have  to  be  performed  according  to 
laws  programmed  in  advance.  For  instance,  the  time  allowed  for  reaching  the 
new  heading,  the  maximum  angular  velocity  of  the  rudder,  the  magnitude  of 
the  possible  overshoot,  are  assigned,  so  that  the  ship  must  maneuver  as  she 
were  on  a rigid  trajectory.  This  implies  that  various  types  of  trajectories 
have  been  completed,  that  one  of  them  is  selected  on  the  occasion  of  each 
maneuver  and  that  the  control  system  is  completely  automated. 

For  designing  such  a control  system  and  computing  the  roads  of  the  ship, 
it  is  necessary  to  know  how  the  ship  responds  to  any  given  maneuver.  This 
can  be  done  by  integrating  the  equations  of  the  motion,  which  implies  that 
the  hydrodynamic  forces  exerted  on  the  hull  and  on  the  moving  control  surfaces 
have  been  determined  in  a number  of  conditions. 

In  the  last  Section  of  Chapter  7 mention  has  been  made  of  several 
experimental  methods  which  can  be  used  to  that  end.  Measurements  of  forces 
on  captive  models  in  forced  motions  give  an  important  part  of  the  necessary 
information.  But  the  estimates  so  obtained  obviously  differ  from  the  forces 
which  would  be  exerted  on  free-running  models.  Furthermore,  certain  forces 
are  not  attainable  experimentally.  Resort  to  theory  is  useful  for  obtaining 
at  least  the  order  of  magnitude  of  terms  missing  in  the  equations.  Theory 
is  also  of  interest  for  yielding  the  structure  of  the  equations.  In  case  of 
discrepancy  between  measured  and  calculated  forces  or  between  observed  and 
computed  roads,  one  might  be  tempted  to  include  into  the  expression  of 
forces  terms  the  theoretical  examination  of  which  would  reveal  that  they  are 
physically  meaningless.  There  is  a risk  that  the  equations  after  modification 
lead  to  worse  results  in  some  other  conditions. 

Vortex  theory  is  the  only  one  which  explains  how  force  orthogonal  to 
the  direction  of  the  motion  can  be  created  in  an  almost  inviscid  liquid.  It 
therefore  seemed  to  me  that  vortex  theory  should  be  the  primary  tool  for 
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building  up  a mathematical  maneuverability  theory.  This  is  the  problem 
investigated  in  the  present  lectures. 

Before  attempting  to  extend  the  vortex  theory  to  ships,  I have  believed 
it  important  to  show  that  it  permits  to  evaluate  the  system  of  forces  due  to 
the  inertia  of  the  liquid,  that  is  the  forces  due  to  the  so-called  added 
masses.  The  classical  theory  of  added  masses  is  briefly  recalled  in  chapter 
2,  while  chapter  3 treates  in  detail  the  corresponding  vortex  theory  from 
both  kinematic  and  dynamic  points  of  view.  Attention  is  to  be  paid  to  the 
splitting  up  of  the  vortex  distribution  into  two  parts.  One  induces  zero 
velocity  field  outside  the  body,  whereas  the  second  generates  outside  the 
velocity  potential  due  to  its  motion. 

Chapter  4 extends  the  wing  theory  to  wings  with  a finite  thickness 
and  small  aspect  ratio  and  to  bodies  in  a motion  of  translation.  In  the  case 
of  a relatively  thick  wing  the  Kutta  condition  is  interpreted  as  due  to  the 
shock  between  streamlines  coming  from  the  suction  side  and  those  coming  from 
the  pressure  side.  When  the  flow  is  steady,  the  relative  velocities  on  both 
sides  are  equal  and  opposite  at  the  trailing  edge.  This  paradoxical  result 
is  due  to  the  fact  that  the  thickness  of  the  boundary  layer  is  infinitely 
small  in  an  almost  inviscid  liquid. 

The  problems  relating  to  hips  in  a uniform,  steady,  oblique  motion  of 
translation  give  rise  to  an  analysis  more  intricate  than  that  of  wings 
because  there  generally  exists  on  the  hull  of  a ship  no  line  which  can  be 
considered  as  the  trailing  edge  of  a wing,  at  least  from  a geometrical  stand- 
point. However,  in  the  case  of  a submarine  moving  in  its  vertical  plane  of 
symmetry,  the  visualization  of  the  flow  seems  to  show  that  the  relative 
streamlines  passing  on  the  upper  side  and  those  passing  on  the  bottom  inter- 
sect along  a U-shaped  line  open  in  the  direction  of  the  motion.  This  turns 
out  to  generate  shocks  and  thus  the  shedding  of  a free  vortex  sheet.  This 
phenomenon  is  analogous  to  the  one  just  been  mentioned  in  the  case  of  thick 
wings  (but,  in  the  present  case,  for  calculating  the  flow  around  the  ship, 
the  place  on  the  hull  of  its  "trailing"  edge  should  be  determined  after 
visualizing  the  flow). 

I hope  that  the  reader  will  easily  accept  my  explanation  of  the  existence 
of  a vortex  sheet  downstream  from  the  body.  As  in  the  case  of  a wing,  the 
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problem  would  remain  unsolved  If  no  additional  condition  were  prescribed. 

This  is  what  we  call  the  generalized  Kutta  condition.  It  states  that  the 
pressure  on  the  hull  is  continuous  through  the  shedding  line.  It  takes  the 
form  of  an  integral-differential  equation  similar  to  that  encountered  in 
the  case  of  wings  with  a finite  aspect  ratio,  but  it  is  sensibly  more 
complicated  and  it  would  probably  be  very  difficult  to  solve  mathematically. 
In  practice  this  should  not  be  of  very  great  importance  since  all  the 
necessary  information  could  be  derived  from  measurements  and  observations 
carried  out  on  models  in  uniform  translation  (provided  that  several  angles 
of  heel  are  used). 

In  the  case  of  oblique  translation  in  the  horizontal  plane,  the  stream- 
lines on  the  port  and  starboard  sides  also  seem  to  intersect  so  that  free 
vortices  are  shed  from  the  rear  part  of  the  keel  and  from  the  stern  post. 

The  same  scheme  also  applies  to  submarines  moving  in  a horizontal  plane.  A 
generalized  Kutta  condition  determines  the  problem  completely.  Although 
apparently  simpler  than  the  condition  for  motions  in  the  vertical  plane,  it 
would  probably  lead  to  many  computational  difficulties.  But  pure  forced 
swaying  motions  should  also  yield  all  the  information  needed  for  practical 
purposes  (provided  that  the  tests  include  several  values  of  the  trim  and 
heel  angles). 

Chapter  5 was  initially  planned  to  deal  with  surface  ships.  The  arcs  of 
vortex  filaments,  the  union  of  which  forms  the  vortex  sheet  located  on  the 
hull,  pearce  the  free  surface.  To  determine  the  velocity  induced  by  the 
vorticity  a vector  potential  is  needed.  If  the  free-surface  condition  is 
linearized,  it  is  possible  to  continue  the  real  vortex  field  above  the  free 
surface  by  taking  its  mirror  image  with  respect  to  the  plane  of  the  free- 
surface  at  rest.  In  his  thesis,  Casal  hfld  reduced  the  part  played  by  the 
free  surface  to  that  of  mirror.  Furthermore,  he  had  assumed  that  the  ship 
behaved  as  an  infinitely  thin  flat  wing,  the  bound  vortex  filaments  being 
vertical  segments  of  straight  lines  in  the  centerplane.  By  means  of  a 
complementary  assumption  concerning  the  directions  of  the  free  vortices  near 
the  hull,  he  obtained  an  approximate  solution  with  many  interesting  features. 
But  the  boundary  condition  at  the  hull  was  satisfied  only  along  the  inter- 
section of  the  centerplane  with  the  plane  of  the  free  surface  at  rest. 
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Casal's  scheme  therefore  implied  that  the  draft  was  very  small.  We  have 
shown  here  that  it  is  possible  to  satisfy  the  boundary  condition  on  the 
whole  surface  of  the  centerplane,  but  this  requires  a Kutta  condition. 
Because  of  the  smallness  of  the  beam/ length  ratio,  the  procedure  sketched 
at  the  end  of  Chapter  4 for  double  models  is  no  longer  usable.  Finally, 
we  have  arrived  to  the  conclusion  that  the  approximation  of  infinitely  thin 
hulls  is  much  too  crude.  For  this  reason,  we  have  come  back  to  normal  ship 
forms  in  the  second  part  of  Chapter  5.  It  has  been  shown  that  the  assumption 
that  the  free  surface  condition  can  be  linearized  leads  to  a solution 
involving  a Kelvin  singularity  distribution  on  the  contour  of  the  waterline. 

Chapter  6 tackled  the  problem  of  small  motions  first  in  a nearly 
vertical  plane  and  then  in  a nearly  horizontal  plane,  i.e.,  the  so-called 
motions  of  the  first  and  the  second  kinds,  respectively.  Of  course,  the 
generalized  Kutta  condition  gives  for  dT/dt  along  the  shedding  line  an 
integral  equation  of  mixed  Volterra  and  Fredholm  types.  Despite  the  more 
complicated  structure  of  the  equations,  small  motions  of  the  first  kind  are 
on  the  whole  similar  to  small  motions  about  a wing  with  a finite  aspect 
ratio.  The  vortex  wake  consists  of  a vortex  sheet. 

For  an  analogous  phenomenon  to  happen  in  the  case  of  small  motions  of 
the  second  kind,  the  free  vortices  should  be  shed  from  the  sternpost  alone. 
But  it  is  not  so,  in  general.  For  instance,  in  a pure  swaying  motion,  they 
are  shed  from  a narrow  region  close  to  the  keel  line,  alternately  to  port 
or  to  starboard,  according  as  the  component  u^  on  the  transverse  axis  of 
the  velocity  of  the  center  of  gravity  is  in  the  negative  or  in  the  positive 
direction.  Consequently  the  vortex  wake  does  no  longer  consist  of  a surface, 
but  of  two  sequences  of  domains,  each  of  finite  volume,  one  created  to  port 
and  the  other  to  starboard.  In  a pure  yawing  motion  u^  is  zero  at  any  time, 
but  the  angular  velocity  u^  is  alternately  positive  and  negative.  When  it  is 
positive,  the  incident  velocity  comes  from  the  starboard  side  at  the  fore- 
body and  from  the  port  side  at  the  afterbody.  During  each  time  interval 

corresponding  to  u,  >0,  two  vortex  wakes  of  finite  volume  are  shed  to 
b 

starboard  from  the  forebody  and  to  port  from  the  afterbody.  When  u,  < 0, 

D 

the  volumic  wake  shed  from  the  forebody  is  to  port  and  that  shed  from  the 
afterbody  is  to  starboard. 


125 


Chapter  7 deals  with  the  interaction  between  hull  and  appendages, 
particularly  those  playing  an  important  part  in  maneuver:  rudder,  diving 
planes,  fixed  fins.  Vortex  theory  explains  why  certain  fittings  may  in- 
crease their  efficiency  and  make  the  response  of  the  ship  quicker.  Tentative 
explanations  have  also  been  proposed  as  for  unexpected  phenomena  (snap- 
roll,  pitch  up  of  submarines).  The  last  Section  presents  a brief  survey  of 
experimental  methods  used  in  the  field  of  ship  maneuverability. 

Let  us  come  back  now  to  the  fundamentals  of  maneuverability  theory 
considered  as  an  application  of  the  vortex  cheory. 

In  two-dimensional  flows  ground  wing  profiles,  the  key  to  the  problem 
is  yielded  by  the  well-known  Kutta-condition.  Or  j might  justify  it  by 
putting  forward  that,  if  this  condition  were  not  satisfied,  the  velocity 
at  the  trailing  edge  would  be  infinite.  However  such  a reason  has  not 
always  seemed  to  be  decisive.  I remember  to  have  heard  Prandtl  speaking 
about  this  question  at  the  3rd  Congress  of  Applied  Mechancis  (Stockholm, 

1931).  What  has  entailed  the  firm  belief  of  the  attendents,  has  been  a 
moving  picture  showing  how  a vortex  appears  at  the  trailing  edge,  how  it 
grows  and  then  travels  with  the  flow,  followed  by  a vortex  wake  which 
progressively  vanishes  far  downstream  from  the  profile  when  the  circulation 
around  the  latter  becomes  a constant.  If  the  fluid  were  inviscid,  no  vortex 
could  raise,  except  in  case  of  shocks  (Helmholtz  theorem).  In  a real  fluid, 
however,  particles  coming  from  the  pressure  side  tend  to  wind  around  the 
trailing  edge  instead  of  flowing  back  along  the  suction  side  in  the  direction 
of  the  rear  stagnation  point.  So  a shock  happens.  The  phenomenon  is  connected 
with  viscosity  effects.  Lagrange's  theorem  according  to  which  the  circulation 
in  any  closed  fluid  circuit  is  a constant  applies  in  the  regions  where  the 
viscosity  effects  are  negligibly  small.  So  the  circulation  in  a fluid 
circuit  surrounding  the  profile  and  the  vortex  wake  born  at  the  trailing 
edge  remains  zero  indefinitely. 

In  the  case  of  a non  infinitely  thin  wing  with  a finite  aspect  ratio, 
it  has  been  seen  in  Chapter  4 that,  because  of  Stokes'  theorem,  the  velocities 
on  the  two  sides  of  the  wing,  near  the  trailing  edge,  are  in  opposite 
directions.  When  the  aspect  ratio  tends  to  infinity,  both  tend  to  zero,  which 
gives  again  the  Kutta  condition.  But  before  the  limiting  situation  is  reach- 
ed, shocks  have  appeared  at  the  trailing  edge. 
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The  assumption  underlying  the  theoretical  considerations  developed  in 
the  present  lectures  is  that  the  existence  of  a lifting  force  on  a ship  hull 
implies  that  of  free  vortex  filaments  shed  from  a line  acting  like  the 
trailing  edge  of  a lifting  surface.  But  the  theory  does  not  yield  the  exact 
position  of  this  line.  The  requirement  to  be  met  is  the  continuity  of  the 
pressure  across  the  free  vortex  sheet,  close  to  the  hull.  The  equation 
expressing  this  condition  depends  on  the  hull  geometry  and  also  on  the 
motion.  There  is  no  need  to  draw  the  attention  to  the  formidable  intricacy 
of  the  computational  problem. 

Nevertheless,  vortex  theory  knows  some  successes  of  importance  in  the 
field  of  ship  maneuverability.  For  instance  the  course  stability  means  that 
the  hull  in  forced  gyration  at  zero  drift  angle  experiences  a lift  which  is 
not  much  smaller  than  the  centrifugal  force.  Otherwise  the  resultant  of  the 
two  could  not  pearce  the  centerplane  ahead  of  the  intersection  with  the 
latter  of  the  lift  in  a translation  at  a very  small  drift  angle.  The  problem 
related  to  course  stability  could  not  be  tackled  by  the  older  methods  re- 
called at  the  end  of  Chapter  2. 

To  conclude  these  lectures,  I should  like  to  say  that  my  firm  opinion 
is  that  neither  pure  empiricism  alone,  nor  pure  theory  alone  can  provide  us 
with  means  sufficient  for  resolving  all  the  problems  relevant  to  ship 
maneuverability.  But  I am  convinced  that  by  making  a judicious  use  of 
experiments  and  of  theoretical  schemes,  many  points  still  rather  obscure 
will  be  gradually  removed. 
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Figure  5 - Vortex  Distribution  on  a Rectangular  Wing  with  a 

Finite  Aspect  Ratio 


Lj  and  are  two  bound  vortex  ribbons  whose  edges  are  the  vortex  filaments 
Lj,  Lj  and  L^,  L ^ respectively.  L^.  is  a free  vortex  ribbon  whose  edges  are 
the  vortex  filaments  and  L^.  L ^ and  L ^ are  closed  at  infinity  downstream 
from  the  wing.  The  intensities  of  1^  and  L ^ sum  up  to  the  intensity  of  L 
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Figure  6 - Free  Vortex  Distribution  for  an  Axisymmetr ic  Body 


Figure  7 - Shocks  between  Upper  and  Lower  Streamlines  on  the 

Hull  of  a Submarine 
(Case  of  a small  angle  of  attack) 


+1 
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Figure  7 a - Families  of  Vortex  Filaments  L^,  L L 
on  the  Hull  of  a Submarine 
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Figure  9 - Lower  Half  of  a Double  Model  of  Surface  Ship  In  Oblique 
Translation  With  a Positive  Drift  Angle  6 
(Bound  vortex  filaments  and  relative  streamlines  on  the  hull) 
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8 : Shedding  line  on  the  hull 


6 : Line  on  E^  infinitely  close  to  8 

Bt : Position  at  t of  the  fluid  line  which  was  on  8r  at  x. 

x f 

E,  E' : Lateral  edges  of  E^.. 

a : Curvilinear  abscissa  on  8 and  6^.. 

B^(o):  Position  at  t of  the  fluid  point  which  was  on  Bj.(a)  at  x. 


Bl(-i),  Mj.  (o'),  M (o) : Points  located  on  the  lower  side  Ef  of  E 

1 *1  1 1 1 

Mf  (o'),  M (o):  Points  located  on  the  upper  side  E,  of  E 

z 1 2 *2  1 2 

Figure  12  - Small  Motion  of  the  First  Kind.  Growth  of  the  Vortex 
Sheet  Behind  a Deeply  Submerged  Body 
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Figure  13  - Double  Model  in  Swaying  Motion  at  te  (t  , t .)  S,  (X,  t ; t) 
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Figure  14  - Double  Model  in  Swaying  Motion  at  te  (t  t 2>  Tht 

Surfaces  S,  ( , t , t + .)  are  now  Quite  Free.  S„(  , t + t) 
1 n n 1 2 n 1 

is  Generated  by  the  Free  Vortex  Arc  Shed  to  Starboard  from 

Bo(  ),  B^(-  ) during  the  Time  Interval  (tn  + j,  t) 
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Figure  15  - Geometry  for  Derivation  of  Generalized  Kutta  Condition 
in  the  Case  of  Small  Motions  of  the  Second  Kind 
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Figure  20  - Small  Unsteady  Motion  of  a Wing  Profile 


Figure  21  - Heeling  Motion  of  a Submarine 
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